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SUMMARY 

PART  A - OPTICAL  TEMPERATURE  MEASUREMENT 

This  study  develops  the  theoretical  background  necessary  for  selecting  the  correct 
spectral  Intervals  for  one  color  and  two  color  optical  pyrometers  used  to  determine  the 
temperature  history  of  explosively  generated  alrblast  environments. 

To  Implement  this  theory  a five  channel  optical  pyrometer  was  designed,  calibrated 
and  tested.  The  design  considerations  separated  Into  three  different  disciplines: 
mechanical  Installation,  electronic  signal  amplification  and  optical  design.  The  system 
consisted  of  two  ultraviolet  channels  at  0.26y  and  O.Sby  wavelengths  to  measure  the  high 
temperature  alrblast  and  three  Infrared  channels  at  1.9)j,  2.9y,  and  4.5u  wavelengths  to 
measure  the  later  arriving  lower  temperature  combustion  products. 

In  addition  to  the  primary  temperature  history  results,  ancillary  data  on  mean  shock 
front  velocity,  estimate  for  depth  of  optical  penetration  of  the  radiometers  and  verifi- 
cation that  the  shock  wave  had  an  Irregular  front  boundary  were  obtained.  The  primary 
results  show  good  agreement  of  measured  temperatures  with  expected  temperatures.  Ihe  two 
Hugonlot  temperatures,  computed  from  the  predicted  peak  pressures  of ~83  atmospheres  (1225 
psla)  and  from  the  "smoothed"  value  of  recorded  peak  pressure  of~67  atmospheres  extra- 
polated from  downstream  data,  completely  bound  the  measured  radiation  temperature  his- 
tories. These  measured  temperatures  suggest  that  the  ratio  of  specific  heats  (y)  strongly 
effects  the  calculation  of  alrblast  temperatures.  Hugonlot  calculation  at  the  above  two 
pressures  is  sensitive  to  variations  or  y»  thus  the  choice  of  y will  Impact  the  Inter- 
pretation of  test  results.  The  results  show  that  the  often  used  y ■ 1.3  mean  specific 
heat  ratio  for  pressure  calculations  at  the  shock  front  also  Is  in  reasonably  good 
agreement  with  temperature  test  results  for  shock  front  pressures  In  the  -1300  psla  to 
-900  psla  range. 

Combustion  product  temperature  Information  was  not  obtained  due  to  electrical  noise 
Interference  on  the  three  Infrared  channels.  The  only  data  that  could  be  verified  was 
the  0.8  msec  predicted  arrival  time  for  combustion  products.  The  ultraviolet  channels 
showed  a steep  drop  In  signal  In  the  0.8  to  1.0  msec  range  Indicating  the  arrival  of 
combustion  products. 

Conclusions  from  the  DABS  test  can  be  summarized  as  follows: 

• Successful  operation  of  the  ultraviolet  and  infrared  channels  was  demonstrated 
In  the  field. 

• Air  shock  temperature  data  was  accurately  measured  by  ultraviolet  pyrometers 
and  correlated  with  theoretical  predictions 

• The  mechanical  design  of  the  open  port  Installation  provided  protection  against 
damage  from  the  shock  wave. 
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The  main  recommendation  for  improved  performance  is  to  redesign  the  electronic 
packaging  to  preclude  any  electronic  noise  interference  with  the  infrared  channels.  The 
second  recommendation  would  be  for  an  overall  package  size  reduction.  This  would  be 
accomplished  by  providing  for  an  optically  transparent  window,  instead  of  the  existing 
open  port  design,  which  would  allow  for  a considerable  size  reduction  and  simplification 
of  installation,  calibration  and  operation. 

PART  B - NOVEL  MEASUREMENT  TECHNIQUES 

This  part  of  the  study  effort  discusses  the  concept  of  deploying  special  techniques 
to  measure  the  density  and  temperature  histories  within  the  Interior  of  a nuclear  driven 
shock  tunnel.  The  extremely  hostile  environment  makes  conventional  Instrumentation 
Ineffective.  The  implementation  of  any  novel  technique  must  first  be  based  on  sound 
physical  principles,  and  secondly  must  lend  itself  to  practical  hardware  design  solu- 
tions. This  study  addresses  the  former  question  by  presenting  a feasibility  analysis 
of  the  proposed  two  techniques.  The  latter  question  while  not  specifically  studied,  is 
addressed  by  listing  a series  of  suggested  areas  of  further  study.  In  this  later  section, 
special  design  problems  are  identified. 

. The  first  technique  discussed  is  the  gas  bag  and  pyrometer  experiment  for  measuring 

f the  temperature  inside  the  shock  wave  as  it  passes  the  observation  station.  The  gas 

bag  is  a device  which  provides  an  optical  "window"  into  the  measured  media  so  that 
passive  optical  pyrometer  systems  can  be  used  to  obtain  temperature  histories. 

i The  second  technique  uses  a radioisotope  source  in  an  active  X-ray  system  in  con- 

ij  junction  with  a detector  system  to  provide  a measure  of  the  areal  density  in  the  path  of 

j observation.  Here,  because  areal  density  measurement  will  depend  on  the  distance  between 

I source  and  detector,  multi-frequency  X-ray  or  y-ray  sources  must  be  used  which  also  will 

i measure  radial  displacement  of  the  tunnel  walls. 

[!  The  feasibility  analysis  presented  in  this  study  shows  that  the  physical  basis  for 

!i  each  of  these  two  concepts  is  sound  and  should  provide  the  desired  Information.  The 

I I study  shows  that  the  gas  bag  technique  can  work  with  the  candidate  cooled  neon  gas 

;i  (-220  °K) . The  radioisotope  source  method  should  give  good  results  using  two  radioactive 

j sources,  one  of  Cobalt  - 60  and  the  other  of  Ytterbium  - 169.  Analyses  of  the  data  pro- 

vided by  these  two  sources  can  yield  histories  of  both  areal  density  and  radial  displace- 
f ment  of  the  tunnel  walls.  There  also  is  hope  of  using  an  additional  Yb  radiosotope 

t X-ray  emission  frequency  to  separate  areal  densities  owing  to  lofted  debris  from  that  of 

the  shock  compressed  air. 
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Future  study  areas  needed  to  evaluate  both  techniques  are  presented  so  that  hard- 
ware design  solutions  can  be  evaluated  and  establish  a high  degree  of  confidence  In  the 
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PART  A 

OPTICAL  TEMPERATURE  MEASUREMENT  j 

1 . INTRODUCTION  i 

1.1  Objective 

The  objective  of  this  program  was  to  demonstrate  the  practical  application  of  using 
optical  radiometers  for  measuring  alrblast  temperatures.  The  dynamic  airblast  simulator 
(DABS)  test  series  uses  a high  explosive  charge  to  generate  the  shockwave.  A airblast 
produced  in  this  manner  has  three  fundamental  structural  components  of  interest; 

• Shock  compressed  air 

• Combustion  product  gases 

• Combustion  product  gases  plus  soot  generated  by  explosives 
Temperature  measurements  from  all  three  regions  are  of  prime  Interest  in  this  program. 

The  theoretical  background  necessary  for  determining  the  proper  wavelength  selection 
is  presented.  A system  required  to  measure  these  three  regions  must  include  both  ultra- 
violet and  Infrared  wavelengths.  In  addition  the  design,  development  and  calibration  as 
it  relates  to  the  fielding  of  the  portable  system  is  outlined.  The  three  disciplines 
which  encompass  this  phase  of  the  program  were  mechanical  installation,  optical  design, 
and  electronic  amplification  of  detector  signal.  Problems  associated  with  each  phase  of 
these  design  parameters  is  discussed  in  the  body  of  this  report. 

1 . 2 Background 

During  discussions  with  DNA/SPSS  and  AFWL/DED  it  was  decided  that  new  instrumentation 
systems  should  be  developed  to  support  the  upcoming  DABS  II/HAVE  HOST  field  tests.  These 
two  test  series  are  in  support  of  the  MX  full  scale  development  program.  HAVE  HOST 
involves  the  interior  and  exterior  air  blast  loading  of  scale  model  trench/shelter 
structures  for  MX  concept  validation,  while  the  DABS  II  test  series  developed  the  phen- 
omenology needed  to  simulate  the  dynamic  air  blast  nuclear  environment  for  MX  related 
external  structural  loading.  In  particular  for  the  DABS  II  test,  there  is  a need  to 
measure  flow  properties  associated  with  the  measured  dynamic  pressure  wavefront  anomolles. 

Thus,  air  blast  patterns  over  the  target  area  need  to  be  determined  along  with  temperature 
history  at  various  ranges. 

During  the  course  of  this  program  the  cross-sectional  dimension  of  the  facility  had 
been  changed  from  a 0.191m. (3  ft.)  diameter  tube  to  the  present  configuration  discribed 
below.  Since  much  of  the  theoretical  development  was  based  on  the  smaller  diameter  test 
[i  facility,  and  extrapolated  to  the  larger  configuration,  this  report  will  make  a clear 

distinction  between  the  two  DABS  II-D  configurations  by  adding  a -1  or  -2  designation  to 
the  primary  callout. 
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The  final  AFWL  DABS  II-D  [1]  (TRW  DABS  II-D-2)  configuration  consisted  of  an  open  ended 
trench  with  parallel  concrete  walls  1.83m  (6  ft.)  high  and  compacted  soil  floor  4.27m. 

(14  ft.)  wide.  The  facility  was  inclosed  with  a prismatic  roof  1.07m.  (3.5  ft.)  high 
from  top  of  concrete  wall  to  roof  peak.  The  whole  roof  structure  was  covered  with  a sand 
overburden.  The  facility  measured  36.2m.  (119  ft.)  long,  with  a steel  reinforced  concrete 
driver  section  at  the  closed  end.  This  section  contained  the  explosive  charge  which  was 
distributed  over  71%  of  the  rear  wall  cross-sectional  area.  A TNT  sensitized  ammonium 
nitrate  castable  slurry  explosive  called  IRETOL  30T22-C  was  used  to  generate  the  shockwave. 
The  total  weight  of  explosive  was  823.7  kg.  (1816  lbs.).  This  charge,  when  detonated, 
produced  a planer  airblast  wave  front  which  proprogates  down  the  length  of  the  facility 
at  a peak  initial  pressure  of~163  atmospheres  (2400  psia). 

The  temperature  at  a given  station  of  the  airblast  shock  front  can  be  calculated  with 
good  accuracy,  based  on  either  shock  front  velocity  or  shock  front  pressure  measurements 
using  the  Rankine-Hugoniot  relationship.  The  temperature  in  the  airblast  after  passage 
geometry,  and  the  multl-dlmenslonal  effects  of  the  surrounding  wall  boundaries  interacting 
with  the  shock  front  and  the  following  high  speed  airflow.  Recent  advancements  in  compu- 
tational modeling  Include  the  Influence  of  mass  addition  and  wall  expansion  effects  on 
calculated  temperature  history. 

Combustion  products  temperature  history  can  be  calculated  based  upon  isentropic 
expansion  from  the  explosive  driver's  Chapman-Jouguet  detonation  point.  Results  calculated 
from  this  relationship  can  vary  since  some  of  the  explosives  driving  the  air  shock  do  not 
completely  detonate  and  continue  burning  downstream  of  the  detonation  point.  In  addition, 
some  of  the  gaseous  chemical  species  are  not  in  equilibrium  at  each  pressure  and  tempera- 
ture state  during  the  available  residence  time  within  the  shock  tunnel.  Further  complica- 
tions are  present  because  the  DABS  II-D  explosive  had  have  additives  that  result  in  fine 
particles  of  soot  mixing  with  the  combustion  products.  These  fine  particles  must  be 
considered  in  designing  any  temperature  measurement  system. 

All  these  factors  introduce  uncertainty  into  the  temperature  calculation.  This 
calculation  is  important,  since  a significant  of  the  blast  environment  around  a target  is 
obtained  by  Interaction  with  combustion  product  gases. 

1.3  Measurement  Environment 

The  measurement  environment  consists  of  air,  chemical  species  and  soot  at  continu- 
ously varying  pressure  and  temperature  states.  The  airblast  peak  pressure  is  established 
by  the  weight  of  explosives,  explosive  type,  distance  from  explosive,  and  the  length  to 
diameter  ratio  of  the  shock  tunnel  as  defined  in  section  1.2.  For  the  DABS/HAVE  HOST  test 
series  the  airblast  pressures  of  Interest  range  from  ~163  to  ~41  standard  atmospheres 
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(2400  to  500  psla) . The  corresponding  time  to  half-pressure  can  be  as  small  as 
1 millisecond  with  times  to  zero-pressure  from  5 to  10  milliseconds. 

Any  ideal  alrblast  shock  front  can  be  represented  by  the  Ranklne-Hugonlot  relation- 
ship 


(1.1) 


The  dashed  line  in  Figure  1-1  represents  this  equation  for  a range  of  peak  pressure  (P^) 
alrblast  shock  fronts.  The  solid  line  curves  branching  off  from  the  Rankine-Hugonlot 
curve  represent  the  Isentroplc  expansion  of  the  gases  as  defined  by 


T 

s 


(1.2) 


For  large  pressure  changes  the  isentroplc  expansion  results  in  corresponding  small 
changes  in  temperature.  For  example,  note  that  for  200  atmosphere  shock  front  pressures, 
the  temperature  is  almost  1 ev  (one  electron  volt  = 11,604.9  °K) . After  severe  isen- 
troplc expansion  to  a 75  percent  reduction  in  shock  pressure,  the  temperature  only  drops 
by  33  percent. 

The  combustion  process  is  controlled  by  the  Chapman-Jouguet  detonation  point 
(Pj,_j  “ 24,000  ATM,  Tj,_j  = 2500  "K).  These  combustion  gases  Initially  arrive  with  the 
contact  surface  after  the  pressure  has  fallen  from  its  shock  front  pressure.  For  the 
case  where  the  pressure  at  the  shock  front  is  200  atmospheres,  the  combustion  products 
may  arrive  at  100  atmospheres.  This  represents  an  isentroplc  expansion  of  the  combustion 
gases  from  the  Chapman-Jouguet  state  of  24,000  atmospheres  to  a contact  surface  pressure 
of  100  atmospheres.  Moreover,  the  value  of  y is  unknown,  but  is  in  the  range  of  1.20  to 
1.35  [2],  and  strongly  Influences  the  isentroplc  expansion.  Figure  1-2  summarizes  contact 
surface  pressure  and  temperature  combinations  for  isentroplc  expansion  from  the  Chapman- 
Jouguet  point  for  different  values  of  y.  If  the  combustion  products  consist  mostly  of 
water  vapor  and  nitrogen  gases  (y  • 1.25),  the  predicted  lower  bound  contact  surface 
temperatures,  neglecting  afterburning,  stay  below  ~900  °K  for  pressures  below  ~150 
atmospheres. 
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with  the  pressure  and  temperature  limits  of  Interest  defined,  additional  information 
relating  to  the  combustion  constituents  must  be  defined  for  radiometer  spectral  windows 
corresponding  to  black  body  radiation.  Because  of  the  time  scale  Involved  it  is  assumed 
that  the  chemical  components  created  at  detonation  will  be  the  same  constituents  propagated 
down  the  shock  tube,  or  more  precisely  a frozen  equilibrium  state  exists.  The  explosive 
manufacturer,  IRETOL,  states  in  their  specification  sheet  that  the  mole-fraction  of  gaseous 
combustion  products  generated  at  the  2500  °K  Chapman-Jouguet  detonation  point  will  be: 

• Water  Vapor  0.64 

• Nitrogen  0.22 

• Carbon  Dioxide  0.13 

At  temperatures  less  than  2500  "K,  the  above  gases  do  not  dissociate. 

In  addition  to  the  gaseous  combustion  products,  explosives  laced  with  TNT,  which  are 
fuel-rich,  can  contain  as  much  as  0.03  mole-fraction  of  gas  products  in  the  form  of 
carbon  soot.  These  soot  particles  will  arrive  behind  the  air-combustion  contact  surface 
and  will  maintain  the  same  temperature  of  the  surrounding  combustion  gases. 

1.3.1  Predicted  Station  Temperature  and  Pressure  Values 

During  the  life  of  this  study  the  predicted  temperature  and  pressure  values  had  con- 
tinually been  updated  by  Mr.  Joe  Renick  of  AFWL  as  a result  of  data  obtained  from  previous 
DABS  tests  and  computer  models.  In  addition  the  original  DABS  II-D  shock  tunnel  configura- 
tion was  updated.  Since  much  of  the  theoretical  development  is  based  on  the  0.91m. 
diameter  version  of  the  DABS  II-D  shock  tunnel  configuration,  while  all  the  design  and  test 
results  refer  to  the  final  rectangular  configuration,  this  report  will  make  a clear 
distention  between  the  two  DABS  II-D  configurations. 

The  DABS  II-D-1  designation  will  refer  to  the  original  AFWL  defined  February  1977 
configuration  which  had  a three-foot  diameter  cross  section  and  two  measurement  stations, 
one  close-in  at  a ~122  atmospheres  (1800  psla)  station  and  a second  downstream  ~81  atmos- 
pheres (1200  psla)  station.  This  shot  was  cancelled  by  AFWL  in  May  1977,  and  later 
replaced  by  the  DABS  II-D-2  configuration.  This  latter  shot  had  a single  measurement 
system  at  the  ~83  atmospheres  (1225  psla)  station  (range  of  8.2  meters).  All  of  the 
theoretical  work  presented  in  this  report  was  initially  developed  for  the  DABS  II-D-1 
configuration,  but  is  readily  transferable  to  DABS  II-D-2.  Since  the  same  explosive  was 
planned  for  both  shots,  the  temperature  and  pressure  domain  of  the  latter  is  completely 
within  that  predicted  from  the  former.  Originally  two  optical  temperature  measurement 
stations  were  planned  for  this  DABS  II-D-1  shot,  corresponding  to  initial  shock  front 
arrival  pressures  of  ~122  and  81  atmospheres  (1800  and  1200  psia).  Owing  to  the  short 
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procurement  time,  existing  DABS  Il-D-1  hardware  had  to  be  utilized  for  the  DABS  II-D-2 
shot.  This  compromise  situation  only  allowed  for  the  fielding  of  a single  station  cor- 
responding to  a shock  front  arrival  pressure  of  ~S3  atmospheres  (1225  psla). 


For  the  original  DABS  II-D-1  shot  the  close-in  station  1 was  primarily  for  measure- 
ing  combustion  products  while  the  further  downstream  station  2,  which  provided  a longer 
duration  between  shock  front  arrival  and  contact  surface  arrival,  would  provide  a 
clearer  picture  of  airblast  temperatures.  This  can  be  clearly  seen  in  Figure  1-3  where 
shock  front  and  contact  surface  arrival  times  are  plotted  with  respect  to  range.  The 
DABS  II-D-2  location  is  situated  near  the  original  DABS  II-D-1  station  2 location,  thus 
providing  a time  duration  of  sufficient  length  to  clearly  record  the  airblast  temperature, 
and  with  the  remaining  recording  time  to  measure  the  combustion  product  temperatures. 

Figure  1-4  presents  TRW-predicted  temperatures  for  the  shocked  air  front  at  range  ' 

dependent  initial  pressures,  , and  after  a isentropic  pressure  decay  to  0.5  P^.  It 
is  assumed  that  the  airblast  combustion  gases  contact  surface  will  pass  the  stations  ^ 


approximately  when  the  pressure  drops  to  0.5  P^;  there  will  be  a pressure  continuity, 
but  with  a sharp  drop  in  temperature.  Owing  to  afterburning,  it  is  futile  to  attempt 
to  predict  temperatures  of  combustion  products  in  the  close  range  portion,  thus  the 
discontinuation  of  the  combustion  product  curves  at  15  meters  in  Figure  1-4.  The  only 
statement  that  can  be  made  is  to  say  that  the  expansion  of  the  gases  should  ensure 
temperatures  lower  than  the  2500  “K  Chapman-Jouguet  detonation  point  temperature.  At 
the  longer  range,  approximately  15  meters,  the  lower  bound  temperature  computed  for  the 
assumption  of  no  afterburning  is  shown  for  the  combustion  gases. 

Table  1.1  gives  intervals  of  pressure,  temperature  and  time  for  stations  1 and  2 of 
DABS  II-D-1  and  the  DABS  II-D-2  station.  The  6t  parameter  is  the  predicted  time  interval 
between  station  arrivals  of  the  shocked  air  front  and  the  combustion  gases  contact  sur- 
face. The  equipment  was  designed  to  survive  and  record  useful  data  for  a total  of  5 
to  6 milliseconds  from  the  onset  of  the  shocked  air  front. 


SHOCK  FRONT  PRESSURE,  Pj 
SHOCK  FRONT  ARRIVAL  TIME 
CONTACT  SURFACE  ARRIVAL  TIME 
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Figure  1-3.  Predicted  Shock  Front  Pressure  and  Arrival  Times  as  a Function  of 
Distance  from  the  Explosive  for  DABS  II-D 
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Figure  1-4.  Predicted  Temperatures  as  a Function  of  Distance  From  The  Explosive  for  DABS  II-D 


Table  1-1.  Predicted  Station  Parameters 


•FOR  7-1.4  (UPPER  BOUND  ESTIMATE) 


2 . THEORY 


The  technical  approach  Is  to  use  the  radiation  emitted  from  the  shock  tunnel  gases 
to  determine  its  instantaneous  temperature.  Namely,  by  measuring  the  radiant  energy 
transfer  rate  per  unit  projected  area  per  unit  sol-*  i angle  and  per  unit  spectral  band 
width  (directed  spectral  radiance,  I),  one  can  determine  the  temperature  of  the  source. 
For  the  particular  case  of  the  thermodynamic  equilibrium  radiation,  the  spectral  radi- 
ance, I,  is  called  Black  Body  spectral  radiance,  1°. 

A hypothetical  black  body  is  completely  opaque  to  all  wavelengths,  X,  of  optical 
radiation  at  all  temperatures,  T.  It  has  an  emlsslvity  e (radiant  energy  emission 
efficiency)  equal  to  unity  for  all  T and  all  X.  The  equations  for  black  body  spectral 
radiance,  1°  (X;  T) , are  presented  in  Appendix  A. 

It  is  Impractical  to  approximate  a true  black  body  except  by  artificial 
construction.  Some  solids  and  liquids  are  said  to  be  gray  bodies.  These  exhibit 
an  emissivity,  e , which  while  less  than  unity,  shows  relatively  slow  variation 

with  temperature  or  wavelength  over  fairly  broad  AX  and  AT  intervals.  Most  real 
bodies,  however,  and  particularly  the  gaseous  volumes  of  current  Interest,  exhibit 
large  fluctuations  of  spectral  absorptivity  and  emissivity  both  with  changing  wavelength 
and  changing  temperature.  Kirchoff's  Law  states  that  the  spectral  emissivity  e(X;  T) 
and  the  spectral  absorptivity  a(X;  T)  for  any  gases  are  equal  for  the  same  X and  T.  In 
this  report  the  curved  parenthesis  will  always  denote  functional  relationships,  l.e., 
£(X;  T),  e is  a function  of  X and  T.  The  spectral  radiance  I(X;  T)  of  the  body  is 
related  to  that  of  a black  body  at  the  same  temperature  by  the  expression: 


I(X;  T)  = e(X; 

T)  I“(X;  T)  (2.1) 

Black  Body 

e = 1 

for  all  X and  for  all  T 

Gray  Body 

£ = e <1 

gray 

£ approximately  constant  over 

gray  ^ 

large  AX  and  AT  regimes 

Gases 

e(X;  T) 

large  variation  with  X and  T,  but 

£(X;  T)  may  = 1 for  selected  AX 

and  AT  regimes 

For  gas  volumes,  e(X;  T)  is  strongly  dependent  upon  the  identity  and  relative  proportions 
of  the  constituent  gaseous  species,  and  also  upon  the  temperature  and  pressure  in  the 
volume  and  the  length  of  optical  path,  L,  through  the  volume.  In  a long  path  length  of 
sufficiently  high  density  and  high  temperature  gas,  spectral  regions  for  a wavelength 
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Interval,  + 1/2  AX,  can  usually  be  found  for  which  e(X,AX  ; T)  =»  1,  and  for  which 
I(X;  T)  « I*  (X;  T). 

2.1  Spectral  Region  Selection  Criteria  for  Optically  Thick  Gas  Volumes 

In  general,  it  is  considerably  easier  to  determine  wavelength  (X),  bandwidth  (AX), 
temperature  (T)  and  temperature  intervals  (AT)  regimes  for  which  a specified  gaseous 
layer  of  known  constituent  distribution,  should  exhibit  a spectral  emissivity,  e(X,  AX;  T) , 
closely  equal  to  unity  than  to  determine  accurate  values  of  spectral  emissivity  which 
depart  appreciably  from  unity.  The  former  approach,  therefore,  was  selected  for 
implementation. 

The  most  general  form  of  the  simple  exponential  decay  law  for  the  transmission  of 
optical  radiation  through  a substance  whose  "optical  thickness"  t increases  from  zero 
with  increasing  depth  of  penetration  L,  is, 

I(X)  = I exp  [-t(X;  L)]  (2.2) 

o 

I(X)  = I (X)  when  t = 0 (i.e.  when  L = 0) 
o 


In  various  forms  and  applications,  the  preceding  equation  is  known  as  Beer's  Law, 

Lambert's  Law,  etc.  While  fairly  valid  for  propagation  through  most  types  of  semi- 
transparent solids  and  liquids,  it  is  in  general  not  directly  applicable  for  propagation 
through  gaseous  layers  for  other  than  very  closely  monochromatic  radiation.  In  the  problem 
of  current  interest,  where  high  pressures  and  high  temperatures  tend  to  blur  the  individual 
spectral  lines  of  the  gaseous  constituents  into  pseudo  continue  (See  Figure  C-3b),  an 
exponential  decay  law  is  again  quite  valid  for  application  to  relatively  narrow  but  finite 
wavelength  intervals  X + 1/2  AX  rather  than  to  strictly  monochromatic  radiation. 

An  overbar,  i.e.,  t,  a,  etc.,  will  be  used  throughout  this  report  to  distinguish  spectral 
parameters  "smoothed"  or  averaged  over  a X+  1/2AX  interval  from  the  corresponding  mono- 
chromatic T,  a,  etc.,  parameters.  Using  this  overbar  convention,  the  mean  spectral 
absorptivity  of  a gaseous  layer  of  spectrally  averaged  optical  thickness  t(X,  AX;  P, 

T;  L)  is  seen  from  the  preceding  equation  to  be  represented  by  the  expression; 

a(X,  AX;  P,  T;  L)  = 1 - exp[-7(X,  AX;  P,  T;  L)]  (2.3) 

In  writing  Equation  (2. 3)  it  was  assumed  that  t and  t are  absorption  rather  than  total 
"extinction"  parameters.  In  any  case,  for  gas  layers,  the  scattering  (or  diffuse 
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reflection)  contributions  to  "extinction"  are  generally  of  importance  only  in  spectral 
regions  of  very  weak  absorption.  If  this  is  true,  then  Equation(2.3)is  valid  as  an 
excellent  approximation  even  for  total  "extinction"  (absorption  plus  scattering)  para- 
meters. 

Then  applying  the  earlier  mentioned  Klrchoffs’  Law  to  Equation  (2. 3)  one  can  write: 


e(X,  AX;  P,  T;  L)  = 1 - exp  [-t(X,  AX  ; P,  T;  L)] 
7(X,  AX;  P,  T;  L) 


Tj^j(x,ax;  P,  T;  L) 


(2.4) 

(2.5) 


l.j 


,th 


with  the  summation  extending  over  each  j spectral  band  in  the  X + 1/2  AX  region  of 
interest  for  each  i*"^  species  of  gas  in  the  gaseous  mixture.  The  individual  "optical 
thickness"  contributions,  turn  be  exhibited  in  the  form: 


(X,  AX;  P,  T;  L)  = k^^  (X,  AX;  T)  - L • [N^(P,T)/N^]  (2.6) 

Here  N (P,T)  is  the  number  density  of  the  i^*^  species  of  radiating  (and  absorbing) 

^ -3  19  3 

molecules,  in  units  of  cm  , while  is  equal  to  2.67  • 10  molecules  of  gas  per  cm 

at  standard  temperature  and  pressure  (STP).  The  optical  path  length  is  L.  With  elimin- 
ation of  the  lower  temperature  boundary  layers  near  the  tube  walls,  L reduces  to  ~81  cm 
for  DABS  II-D-1  and  to  -410  cm  for  DABS  II-D-2.  The  k^^  (X,  AX;  T)  are  spectral  mean 
absorption  coefficients  specified  in  the  units  of  "Inverse  reduced  atmosphere  centi- 
meters" [(atm  cm)  or  often  more  precisely,  (atm  cm  - STP)  most  commonly  found  in 
the  literature.  This  does  not  mean,  as  is  sometimes  mistakenly  assumed,  that  these 
absorption  coefficients  are  evaluated  under  STP  conditions.  The  rationale  behind  the 
preceding  designation  is  that  the  multiplier,  L[N^(P,  T)/N^],  occurring  in  equation(2.6) 
can  be  said  to  represent  an  equivalent  path  length  that  would  result  if  the  1^^  species 
gas  were  purified  and  reduced  to  STP  conditions.  This  readily  can  be  seen  by  first 
introducing  the  mole  fraction,  T),  where; 

(P,  T)  = (P,  T)/N(P,T)  (2.7) 

Here,  N(P,  T)  represents  the  total  number  density  of  all  gas  particles  per  cm^  at 
absolute  temperature  T and  pressure  P.  Then,  for  the  ideal  gas  assumption  which  is 
completely  Justified  in  the  (P,  T)  domain  of  current  interest. 
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P • T 

N^/N  • L i (N/N^)  • L - Xj  • p ' L (2.8) 

o 

The  explicit  (P,  T)  dependence  of  Xj^»  and  N have  been  omitted  In  the  preceding  equation 
for  convenience  and  also,  when  unambiguous,  will  often  be  omitted  from  later  equations 
for  the  same  reason.  The  subscript,  o,  on  P,  T or  N,  etc.,  will  always  Indicate  STP 
conditions.  Thus,  since  the  (P/P^)  ratio  gives  the  pressure  In  units  of  standard  atmos- 
pheres, the  terms  on  both  sides  of  Equation  (2.8)  are  said  to  be  In  units  of  effective 
reduced  (atm  cm  - STP). 

A substitution  of  Equation  (2.8)  Into  Equation  (2.6)  yields: 

^^j(X,  AX;  P,  T;  L)  = Xi(P.  T)  k.j(X,  AX;  T)  • P/P^ * T^/T  • L (2.6a) 

The  preceding  equation  Is  used  when  the  mole-fractions,  Xj^»  rather  than  the  species 
number  densities,  N^,  are  the  directly  available  data.  The  spectrally  smoothed  absorption 
coefficients  for  the  molecular  and  Ionic  species  of  Interest  generally  were  either 
directly  available  In  the  literature  or  could  be  calculated  from  data  presented  In  the 
literature. 

The  fundamental  absorption  parameters  are  the  monochromatic  "per  molecule  cross 
sections",  a(X),  rather  than  the  somewhat  arbitrarily  defined  monochromatic  absorption 
coefficient  counterparts,  k^j(X;  T),  or  the  k(X,  AX;  T)  parameters  introduced  in 
equation  (2.6).  The  fundamental  a-cross  sections  are  In  general  only  weakly  temperature 
dependent,  but  the  k-coef ficlents  can  exhibit  from  one  to  several  orders  of  magnitude 
variation  with  temperature  within  the  temperature  domain  of  current  Interest.  The 
reason  for  this  variation  Is  related  to  the  dependence  of  the  k-(and  k-)  parameters 
upon  "excited  state  populations".  A brief  discussion  of  the  relation  between  the  k-  and 
a-parameters,  and  of  the  temperature  dependence  of  the  former  Is  given  In  Appendix  C. 

To  recap,  a large  fraction  of  the  Initial  effort  was  spent  In  determining  spectral 
regions  for  which  either  air  or  the  combustion  products  could  be  expected  to  exhibit 
optical  thickness,  t ^ 3.  From  equation  (2.4)  it  can  be  seen  that  e ^ 0.95  for  t ^ 3. 

Equation  (2.6)  or  equation  (2.6a)  was  used  to  predict  (X,  AX;  T)  values  from  N^,  or 
X^»  and  k (X,  AX;  T)  values  that  were  either  available  in  the  literature  or  could  be 
computed  from  the  In-house  computer  program  TRW  ACE. 
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2.2  Spectral  Region  Selection 


The  spectral  intervals  X + 1/2  AX  investigated  were  all  those  between  ~0 .25u  and 
~5y  wavelength  for  which  the  respective  constituent  gases  (and  carbon  soot  solids  for 
the  combustion  products)  were  known  to  exhibit  elevated  temperature  absorption  bands 
that  gave  promise  of  possibly  being  useful  and  were  not  excessively  expensive  to  imple- 
ment. Optical  instrumentation  to  view  the  very  strong  ISp  CO2  band,  for  example,  is 
approximately  a factor  of  10  more  expensive  than  for  wavelengths  below  ~5)j.  Also  these 
detectors  must  use  liquid  helium  cryostats.  The  longer  wavelength  apparatus  also  is  very 
time-consuming  and  difficult  to  adjust  in  the  field.  The  same  applies  to  ultraviolet 
radiation  at  wavelengths  much  below  0.25y  where  the  hot  air  layer  viewed  would  certainly 
be  "optically  thick".  Moreover,  normal  temperature  ambient  air  is  completely  opaque 
below  ~0.21y  wavelength  for  even  a few  centimeter  path  lengths,  and  exhibits  very  strong 
absorption  near  15y  wavelength  for  moderate  path  lengths.  The  high  temperature  radia- 
tion being  investigated,  therefore,  would  be  at  least  partially  absorbed  in  the  far 
ultraviolet  or  in  the  ISy  wavelength  CO2  region. 

2 . 3 One  and  Two  Color  Pyrometer  Systems 

The  output  voltage  V(X,  AX;  T)  of  a radiometer  filtered  to  respond  to  a X + 1/2  AX 
wavelength  increment  is  directly  proportional  to  the  spectral  radiance  I(X,  AX;  T) 
being  viewed.  The  proportionality  factor,  which  is  dependent  upon  the  optical  system, 
detector  cell  and  amplifier  parameters,  is  subject  to  either  computation  or  preferably 
to  empirical  absolute  calibration.  If  the  spectral  emissivlty  e(X,  AX;  T)  is  known 
to  be  wl  [so  that  I(X,  AX;  T)wi°(X,  AX;  T)  ] throughout  the  temperature  regime  of 
Interest,  then  there  is  a unique  one-to-one  relationship  between  I°(X,  AX;  T)  and  T, 
and  hence  between  V(X,  AX;  T)  and  T.  Equation  CA.2b)  in  Appendix  A gives  the  Black 
Body  law  relation  between  I°(X,  AX;  T)  and  T.  This  equation  could  be  solved  for  a 
somewhat  awkward  relation  expressing  T as  a function  of  1°,  and  hence  a function  of 
the  directly  measured  quantity  V. 

It  is  more  convenient,  however,  to  make  calibration  graphs  representing  both 
black  body  radiance  I°(X,  AX;  T)  and  signal  voltage  V^(X,  AX;  T)  as  unique  functions  of 
T as  the  Independent  variable.  Figure  A-3  in  Appendix  A presents  1°  vs.  T curves  for 
five  selected  wavelength  Intervals  X^  + 1/2  AX^.  Figure  6-4,  in  Section  6.0  presents 
some  equivalent  vs.  T curves.  Any  one  of  these  latter  type  of  vs.  T calibration 
curves,  together  with  the  corresponding  X^  + AX^  radiometer  system  is  said  to  con- 
stitute a "one  color"  pyrometer. 


A pair  of  radiometers  filtered  to  respond  respectively  to  + 1/2  AX^  and 
X + 1/2  AX  together  with  a V (X  , AX.;  T)/V  (X,,  AX.;  T)  vs.  T calibration  curve  is 
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said  to  constitute  a "two  color"  pyrometer  system.  The  V ratio  vs.  T curves  are 

s 

similar  to  the  1 ratio  vs.  T curves  shown  in  Figure  A-4  Appendix  A. 

2.3.1  One  Color  Pyrometer  Systems 

It  is  possible  to  make  fairly  reliable  temperature  determinations  at  least  over 
a limited  temperature  interval,  with  a one  color  pyrometer  system  such  as  was  described 
in  the  preceding  discussion.  This  type  of  "one  color"  pyrometer,  however,  suffers  from 
the  disadvantages  of  requiring  accurate  absolute  calibration,  a relatively  narrow 
spectral  region  for  accuracy  and  of  falling  off  fairly  rapidly  in  relative  accuracy 
outside  of  a favorable  temperature  Interval  which  is  a function  of  the  selected  X . As  j 

will  be  shown  in  Appendix  A,  the  ultraviolet  spectral  region  is  most  favorable  for 
the  measurement  of  very  high  temperatures,  while  the  near  Infrared  is  most  favorable 
for  lower  temperatures.  This  latter  weakness  possibly  could  have  been  overcome  by 
the  use  of  several  "one  color"  pyrometer  systems  at  different  wavelengths  X^  + 1/2  AX^ 
each  selected  for  t ^ 3. 

2.3.2  Two  Color  Pyrometer  Systems 

All  three  of  the  preceding  objections  can  be  more  readily  met,  however,  by  the 
use  of  several  "two  color"  pyrometer  systems,  each  consisting  of  a pair  of  radiometers 
for  sufficiently  separated  high  t wavelength  Intervals,  say  X^  + 1/2  AX^  and 
Xj  + 1/2  AXj . As  will  be  seen  in  Appendix  A,  the  "temperature  signals"  for  a "two 
color"  pyrometer  is  the  spectral  radiance  ratio  [I(X^,  AX^;  T)/I(Xj,  AX^ ; T) ] . A "two 
color"  pyrometer,  depending  as  it  does  only  upon  the  ratio  rather  than  upon  the  absolute 
value  of  two  input  signals,  does  not  require  absolute  calibration.  It  suffices  to 
either  have  equivalent  geometries  and  optical  and  electronic  gain  parameters  for  the 
two  systems,  or  to  be  able  to  scale  the  differences  in  these  factors  accurately.  Some- 
what wider  AX  filters  are  feasible,  since  again,  the  temperature  signal  does  not  depend 
so  critically  upon  variations  of  the  individual  received  spectral  radiance  increments, 
but  only  upon  variations  in  their  ratio.  Finally,  as  will  be  seen  from  one  of  the  curves 
shown  in  Figure  A-4,  it  is  possible  to  obtain  acceptable  accuracy  over  a temperature 
extending  from  over  8000  “K  to  less  than  3500  °K  by  selecting  one  detector  system 
centered  on  an  ultraviolet  wavelength  and  the  other  on  a near  infrared  wavelength.  The 
high  and  low  temperature  cut-offs  are  set  more  by  the  cost  and  gain  of  multiple  decade 
amplifiers  and  low  noise  optimum  detectors,  and  the  cost  of  low  noise  transmission  lines, 
than  by  the  decay  of  the  spectral  radiance  ratio  temperature  signal. 
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2. A Selected  Wavelength  Intervals  for  Two  Color  Pyrometers 
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As  was  indicated  in  sub-section  2.2,  the  choice  of  wavelength  intervals  to  be 
utilized  was  restricted  to  those  between  ~0. 25u  and  Sy  for  which  the  optical  thicknesses 
t(X,  AX;  T;  P)  could  be  predicted  to  be  ^ 3,  if  possible,  for  the  (T,  P)  regimes  of 
Interest.  Since  the  choices  of  pyrometer  wavelength  intervals  are  dependent  upon  the 
radiating  molecular,  atomic,  ionic  (or  particulate)  species  being  observed,  they  are 
different  for  the  shock  heated  air  volume  than  for  the  later  arriving  combustion  products 
which  consist  primarily  of  H2O,  CO2  and  N2  gases,  and  of  particulate  carbon  soot. 

2.4.1  Wavelength  Selections  for  Shock  Heated  Air 

After  investigation  into  the  high  temperature  optical  thickness  contributions 
Tjj(X,  AX;  T,  P)  for  X between  ~0.25y  and  Sy  to  be  expected  from  air  molecules,  atoms 
and  ions,  it  was  determined  that  only  the  ultraviolet  and  short  wavelength  visible 
spectrum  region  for  X < 0.38y  gave  any  promise  of  sufficient  opacity,  i.e.,  t > 3, 
for  shock  heated  normal  air.  It  was  determined  that  the  best  compromise  wavelength 
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intervals  for  natural  air  were  X^^  as  0.26y  and  X2  « 0.36y,  with  AXj^  w AX2  0.01  (lOOA)  . 

A somewhat  shorter  wavelength  would  have  been  preferred  for  X^^,  but  was  ruled  out  by 
optical  system  considerations.  A longer  wavelength  for  X2  was  ruled  out  by  the  rapid 
fall-off  with  increasing  wavelength  of  hot  air  opacity.  As  it  is,  the  computed  value  of 
T(0.36y;  O.Oly;  T,  P)  falls  below  the  critical  values  of  3 for  pressure  below  ~45 
atmospheres  over  much  of  the  temperature  interval  of  interest  (see  Table  1.1)  for 
DABS  II-D-1;  for  DABS  II-D-2  the  t value  will  be  five  times  greater.  The  (Xj^  “ 0.26y, 

X2  = 0.36y)  pair  are  suitable  primarily  for  temperatures  in  the  ~8000  °K  to~3000  °K 
region  and  for  pressures  above  ~40  or  50  atmospheres  (owing  to  the  diminishing  opacity 
of  the  0.36y  spectral  region). [Pressures  as  low  as  25  atmos.  suitable  for  DABS  II-D-2] 

Outside  of  some  relatively  vary  sharp  atomic  lines  from  N,  0,  Ar,  Na  and  metallic 
impurities,  it  would  appear  that  the  only  strong  air  radiation  for  temperatures  of 
less  than  8000  “K  and  wavelengths  longer  than  O.Ay  comes  from  what  are  known  as  the 
molecular  nitrogen  first  positive  band  system  [N2  (1+)]  and  from  the  atomic  oxygen 
radical  0 (free-bound)  continuum  (see  Appendix  D for  explanation  of  these  terms).  How- 
ever, computations  using  equation  (2.6)  showed  that  even  at  T = 8000  °K  and  P/P^  **  130 
atmospheres,  an  optical  thickness  t value  of  ~1  was  the  greatest  to  be  expected  in  the 
~0.55y  to~0.85y  wavelength  region  for  DABS  II-D-1]  where  this  radiation  is  the  strongest. 
[For  DABS  II-D-2, the  x values  will  again  be  approximately  five  times  greater.) 

Previous  investigations  have  used  the  0.55y  to  0.85y  wavelength  spectral  region 
for  optical  temperature  determinations  of  flames,  explosions,  shock  tube  phenomena,  etc. 
However,  they  may  have  been  looking  at  other  radiating  species  than  those  that 
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exist  In  normal  shock  heated  air,  or  at  higher  pressure  or  longer  optical  paths.  Also 
there  are  expensive  and  difficult  to  implement  optical  techniques  that  are  more  suitable 
for  laboratory  type  environments  than  for  one-shot  field  tests  using  expendable  instru- 
mentation. Some  of  these  latter  techniques  use  relatively  high  spectral  resolution 
interferometers  or  scanning  monochrometers  to  produce  spectra  of  decidedly  non-optically 
dense  gases.  The  "rotational"  or  "vibrational"  temperatures  then  are  determined  from 
the  "back  slopes"  of  the  rotational  structure  of  one  or  more  "bands",  or  from  the 
relative  intensities  of  resolved  vibrational  bands  (see  Appendix  B) . 

2.4.2  Wavelength  Selections  for  4%  CO^  Added  to  Air  [DABS  II-D-1] 

Preliminary  computations  also  showed  quite  conclusively  that  an  81  cm  path  of 
the  hot  air  volume  would  not  achieve  the  desired  opacity  in  the  near  infrared  spectral 
region  between ~0. 85  and  Sp.  The  strongest  absorption  in  this  region,  from  the  CO^  so- 
called  4.3p  wavelength  band,  still  is  too  weak  for  the  available  optical  path  owing 
to  the  very  low  concentration  of  CO2  in  normal  air.  However,  since  it  is  very  desirable 
to  make  pyrometer  measurements  at  near  infrared  wavelengths  in  order  to  measure  tempera- 
tures below  ~4000  °K,  it  was  suggested  that  ~4%  CO2  be  added  to  the  normal  pre-shock 
air,  either  throughout  the  "tube",  or  in  a thin  plastic  "balloon"  containing  the  view 
volumes.  It  is  shown  in  Appendix  E that  this  small  addition  of  CO2  will  produce 
suitable  opacity  not  only  in  the  very  strong  and  broad  CO2  vibrational  "fundamental 
band”  centered  on  4.3vi,  but  also  in  the  much  weaker  vibrational  "combination  band" 
centered  on  2.7y.  Allowing  for  the  necessary  tolerance  in  both  the  central  positions  and 
"widths"  of  relatively  inexpensive  and  quickly  available  off-the-shelf  filters,  the 
best  compromise  positions  for  turned  out  to  be  2.8p  and  4.5p  , respectively, 

with  AX2  arid  AX^  w O.lp  (see  Figure  2-1). 

The  value  of  X^  = 4.5p  rather  than  4.3p  was  selected  because  as  seen  in 
Appendix  D,  CO2  tends  to  dissociate  into  CO  and  0 for  temperatures  much  above  3000  or 
3500  °K.  The  4.5p  wavelength  optical  system  will  view:  CO2  at  the  lower  temperatures 
up  to  ~3500  °K;  CO2  and  CO  at  mid-temperatures  between  ~3500  °K  and  4000  °K;  CO  from 
~4000  °K  to  ~6500  °K;  and  both  CO  and  electron  free-free  [e  (ff)]  radiation  above  ~6500  "K 
to  ~8000  °K.  Also,  the  natural  CO2  content  of  cold  ambient  air  in  the  optical  path 
will  show  very  little  absorption  at  4.5p  and  at  2.8p, 

With  the  4%  CO2  added  to  air,  the  4.5u  wavelength  interval  will  be  "black", 
i.e.,  t^3,  over  virtually  the  full  pressure  interval  of  interest  and  for  1000  °K<,  T 
X 8000  ®K  (and  owing  to  the  e radiation,  up  to  any  higher  temperature  of  interest).  As 
a consequence,  the  4.5p  X^  radiometer  would  pair  not  only  with  the  2.8p  radiometer,  but 
also  with  the  0.26p  Xj^  radiometer  to  cover  the  temperature  interval  from  ~3500  °K  to 
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above  8000  “K.  The  lower  and  upper  temperature  cut-offs  are  set  by  considerations  of 
excess  dynamic  range,  of  transmission  line  noise,  and  of  the  expense  of  extra  high 
quality  low  noise  detectors  and  amplifiers  rather  than  by  any  lack  of  opacity  or  radiation 
signal  fall-off.  For  DABS  II-D-2  the  addition  of  A%  CO2  would  have  allowed  the  2.8u 
and  4.5p  Infrared  wavelengths  to  read  the  alrblast  temperature. 

2.4.3  Wavelength  Selections  for  the  Combustion  Products 

It  Is  relatively  easy  to  show  that  the  carbon  particle  soot  arriving  with  the 
combustion  products  should  be  "black"  at  any  near  Infrared  (and  probably  visible) 
wavelength  even  for  very  low  pressures  of  the  order  of  3 or  4 atmospheres  or  perhaps 
less.  However,  since  this  soot  will  lag  the  contact  surface  of  the  combustion  produce 
and  the  shock  heated  air  by  some  fraction  of  a millisecond  the  infrared  wavelengths  will 
be  appropriate,  thus  wavelengths  = 1.9y  and  X^  = 4.3u  were  selected,  again  with 
ax's  of  ~0.1  u.  The  1.9y  ratlometer  will  be  viewing  a relatively  rather  weak  water 
vapor  vibrational  combination  band,  which  owing  to  the  mole-fraction  of  ~0.64  H2O  in 
the  combustion  products  will  nevertheless  give  values  of  t > 3 over  virtually  the  full 
pressure  range  of  interest,  particularly  since  the  carbon  soot  should  arrive  before  any 
marked  drop-off  In  pressure  occurs. 

The  wavelength  of  X^  « 4.3y  rather  than  the  X^  value  of  ~4.5y  was  Initially 
selected  for  viewing  DABS  II-D-1  combustion  product  radiation  in  order  to  provide  a 
capability  of  measuring  temperatures  of  combustion  gases  below  ~1000  “K  before  the 
arrival  of  the  carbon  soot  particles  within  the  station  fields  of  view.  At  the  time,  it 
was  thought  that  the  lag  time  of  the  soot  particle  arrivals  behind  the  contact  surface 
might  possibly  be  a fraction  of  a millisecond.  The  CO2  absorption  coefficients 
k„«  (X,  AX;  T)  are  relatively  independent  of  temperature  for  X « 4.3y,  but  decay  by 
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several  orders  of  magnitude  as  T decreases  to  ~500  °K  for  X«  4.5y  (particularly  for 
the~0.15y  AX  Increments  which  turned  out  to  be  available  in  off-the-shelf  filters). 

Thus,  the  computed  absorption  optical  thickness  remained  well  above  the  critical  value 

of  3 throughout  the  full  pressure  range  of  interest  for  X «4.3y,  but  dropped  below  this 
value  for  T less  than  ~1000  °K,  or  possibly  ~800  °K,  for  X = 4.5y.  For  DABS  Il-D-2 

the  optical  thickness  x would  have  been  five  times  greater. 

L.U2 

An  undesirable  corallary,  however,  is  that  4.3y  wavelength  radiation,  unlike 
4.5y  wavelength  radiation,  is  strongly  absorbed  by  ambient  air  temperature  CO2.  It  is 
estimated  that  the  ~0.03%  CO2  content  of  air  in  the  ~12  foot  optical  path  would  absorb 
~20  to  30%  of  viewed  ~4.3y  wavelength  radiation.  This  loss  could  have  been  corrected 
for  by  field  calibrations,  or  measures  could  have  been  taken  to  replace  the  ambient  air 
in  the  view  path  by  non-absorbing  dry  nitrogen  gas.  In  view  of  the  short  time  available 
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for  implementation,  however,  the  latter  procedure  was  not  practical  for  the  DABS  Il-D-2 
shot.  Also,  the  purchased  off-the-shelf  filter  turned  out  to  be  centered  on~4.25y, 
where  the  ambient  air  absorption  is  even  stronger  than  for~A.3u. 

The  preceding  ambient  air  absorption  problem  was  resolved  by  a closer  evaluation 
of  the  probable  carbon  soot  lag  time.  It  is  believed  that  unlike  larger  particles, 
the  expected  sub-micron  size  soot  particles  arrive  very  close  behind  the  contact 
surface.  Since  the  arriving  combustion  products  would  be  completely  opaque  at 
all  near  Infrared  wavelengths,  there  turns  out  to  be  no  advantage  in  using  a~4.3P 
rather  than  a 4.5vi  radiometer.  It  thus  was  decided  to  substitute  the  latter  radio- 
meter for  the  DABS  II-D-2  combustion  products  temperature  measurements. 

2.4.4  Summary  of  Wavelength  and  Pyrometer  (Radiometer  Pair)  Selections 

Figure  2-1  summarizes  the  radiometer  wavelength  selections,  through  X^,  and 
their  positioning  at  either  Station  1 or  Station  2,  or  at  both  stations  for  DABS  II-D-1. 

As  seen  in  Table  1.1,  the  close-in  Station  1 for  DABS  II-D-1,  will  view  the  shock 
heated  air  for  only  ~0.15  milliseconds,  and  then  spend  the  remainder  of  its  brief 
lifetime  viewing  combustion  products.  Station  2,  on  the  other  hand,  will  see  the  shock 
heated  air,  now  at  reduced  temperature  and  pressure,  for ~ 0.85  milliseconds. 

The  hot  air  and  air  + 4%  CO2  viewing  radiometers,  filtered  to  respond  to  the 
four  wavelengths  centered  respectively  on  X^^  = 0.26u,  X2  = 0.36u,  X^  = 2.8u  and  X^  = 4.5u 
compromise  three  two-color  pyrometers  when  taken  in  the  pairs  (X^,  X2),  (X^,  X^),  and 
(Xj^,  X^)(see  Figure  A-4).  It  also  is  possible  to  interpret  pyrometer  data  from  the 
pairs  (Xj^,  X^)  and  (X2,  X^),  but  over  more  limited  temperature  intervals;  X2  and  X^ 
do  not  both  view  "black"  radiation  at  overlapping  temperature  Intervals. 

The  primary  combustion  product  viewing  radiometers  originally  were  Intended  to  be 
X,  « 1.9m  and  X,  « 4.3m,  but  for  reasons  outlined  at  the  end  of  the  preceding  sub-section, 

j b 

it  was  decided  to  substitute  X,  = 4.5m  for  X,.  In  addition  to  the  (X.,  X.)  pair  however, 
combustion  product  temperature  Information  also  can  be  extracted  from  the  (X^,  X^) 
pair,  and  (less  accurately)  from  the  (X^,  X^)  pair.  The  (X^,  X^)  pair  should  be 
particularly  useful  for  temperatures  approaching  ~500  where  the  weaker  voltage 
signal  from  the  X^  ss  1.9m  radiometer  may  tend  to  drop  too  close  to  the  noise  level. 

Figure  2-2  summarizes  the  radiometer  selection  data  for  DABS  II-D-2.  As  was 
shown  in  Table  1.1, ~0. 8 milliseconds  will  be  spent  viewing  the  shock  heated  air,  and 
the  remaining  few  milliseconds  of  the  equipment  survival  time  will  be  devoted  to  the 
presumably  cooling  combustion  products. 


2.5  Information  Presented  in  Appendices 
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A detailed  understanding  of  Black  Body  radiation  theory  and  of  molecular  spectro- 
scopy Is  not  required  for  an  understanding  of  this  report,  but  a rudimentary  under- 
standing of  the  concepts  and  nomenclatures  involved  would  be  helpful.  Appendices  A, 

B and  C present  brief  over-views  of  these  topics  for  the  reader  who  may  not  have  been 
introduced  to  them  previously,  or  who  may  wish  to  refresh  his  memory.  Appendix  A 
presents,  without  any  attempt  at  derivation,  the  rudiments  of  Black  Body  radiation 
theory,  and  of  criteria  for  the  applicability  of  optical  radiometer  measurements  for 
determining  the  radiation  temperature  of  a volume  of  hot  gases.  It  also  discusses 
the  selection  of  optimum  spectral  regions  for  observations,  and  presents  generalized 
error  analyses  criteria  for  "one  color"  and  "two  color"  pyrometer  techniques.  Appendix 
B presents  a skeletal  definition  of  some  of  the  nomenclature  of  molecular  spectroscopy. 
Appendix  C presents  a very  brief  discussion  of  the  basic  spectroscopic  parameters  of 
absorption  "cross  sections",  and  of  the  temperature  and  pressure  dependence  of  spectral 
line  contours,  line  "half  widths",  and  the  "spectral  density  of  lines".  It  also  shows 
how  all  of  these  parameters  combine  to  define  the  spectrally  smoothed  absorption 
coefficients  k(X,  AX;  T)  that  generally  are  presented  in  the  literature  and  which  were 
used  in  this  report  to  make  required  gas  volume  spectral  opacity  evaluations. 


2.5.1  Presentation  of  Radiating  Species  Information 


Appendix  D presents  a brief  discussion  of  radiating  species  parameters  for  air. 


air  plus  A%  CO2  and  the  H2O,  CO2  and  carbon  soot  combustion  products  for  the  DABS  II-D 


shots.  The  selected  radiometer  wavelength  intervals  X^  + 1/2  AX^,  were  dictated  by 
the  spectral  positions  and  strengths  of  the  available  radiating  species  absorption 
bands,  together  with  the  practical  system  considerations  mentioned  in  sub-sections 
2.2  and  2.3.  The  principal  radiating  species  parameters  of  interest  are  the  magnitudes 


and  spectral  positions  of  the  spectrally  smoothed  absorption  coefficients  k^j(X,  AX;  T) 


and  the  number  densities  N^(P,  T) , or  mole-fractions  Xj(P,  T) , for  these  specie^  As  a 


corailary,  the  equilibrium  persistence  or  decay  of  these  species  with' increasing  tempera- 
ture also  is  important,  as  is  the  effect  of  increasing  pressure  upon  relative  species 
concentrations.  All  these  topics  are  discussed  and  graphically  presented  in  Appendix  D. 

2.5.2  Presentation  of  Opacity  Data 


As  was  previously  stated,  the  somewhat  arbitrary  opacity  criterion  adopted  in 
this  report  is  that  in  the  selected  spectral  regions,  the  mean  optical  thickness 
t(X,  AX;  P,  T)  should  be  3 so  that  the  corresponding  emlssivlty  e(X,  AX;  P.  T;  L) 
would  be  ^ 0.95%.  Appendix  E graphically  presents  computed  data  to  show  that  this 
criterion  substantially  is  met  for  the  selected  X^  + 1/2  AX^^  spectral  regions  for  air. 
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air  + 4%  CO2  and  combustion  products  for  the  DABS  II-D-1  configuration.  As  previously 
stated,  the  much  longer  path  length  for  DABS  II-D-2  Insures  that  the  corresponding 
T values  for  the  same  (P,  T)  pair  will  be  a factor  of  ~5  higher  for  this  latter 
configuration.  Thus,  for  temperatures  below  ~3000  °K,  the  radiation  data  for  the 
DABS  II-D-2  configuration  could  be  used  to  predict  temperatures  for  pressures  that  are 
lower  by  a factor  of  1/5  than  the  minimum  that  would  have  been  reliable  for  the 
DABS  II-D-1  configuration.  With  the  lower  pressures  expected  at  the  DABS  II-D-2  station, 
species  dissociation  may  tend  to  occur  more  rapidly  and  the  decreased  pressure  factor 
might  only  be  ~l/3  for  temperatures  above  ~5000°K. 


3.  OPTICAL  SYSTEM  PARAMETERS  FOR  THE  DABS  II-D-2  TEST 


The  optical  system  configuration  for  the  DABS  II-D-2  test  is  shown  in  Figure  3-1. 

Note  that  molecules  located  in  the  occlusion  region  (darkly  shaded)  will  be  obstructed 
by  the  aperture  so  that  full  Illumination  of  the  lens  and  detector  is  impossible, while 

molecules  in  the  lightly  shaded  region  will  provide  full  illumination  of  lens  and  detector. 
The  circular  viewing  port  hole  into  the  radiating  gas  volume  is  2.8  cm  (1  1/8  inches) 
in  diameter.  The  objective  lenses  of  all  five  radiometers  are  at  a virtually  equal 

distance,  r,  between  ~308. 8 and  308.2  cm  (~121  1/3  inches)  from  this  aperture,  and  all 

1.6ns 

have  5.08  cm  (2  inch)  diameters,  D,  and  optical  transmissivities,  8 , equal  to 

2 2 -4 

0.90.  The  lenses  thus  subtend  solid  angles,  0,  of  (uD  /4r  ) = 2.134  • 10  steradians 

from  any  point  on  the  aperture.  Since  these  parameters  are  all  virtually  equal,  an 

Isns 

indexing  subscript  would  be  superfluous  for  r,  D,  Q,  and  g . The  spectral  filter 
transmissivities,  g^,  are  all  different,  however,  and  thus  need  the  i-indexlng  subscript. 

The  effective  radiant  power  input  to  the  l*”^  radiometer,  filtered  to  respond  to  a 
narrow  spectral  Interval  centered  on  X^,  is: 

Xi  + m^  (AX^/2) 

pinput  ^ gg^^^plens  . j ^o^^.  ^ 

Xi  - m^  (AX^/2) 

g 

The  symbol  S represents  the  unoccluded  Increment  of  gas  layer  surface  area  which  is 
actually  focused  within  the  square  responsive  area,  A*^  =(0.4  cm)^,  of  each  1^*^  radio- 
meter detector  cell.  The  measured  filter  transmissivities  are  shown  in  Figure  3-2a, 
b,  c,  d and  e.  It  was  found  that  the  relatively  inexpensive  off-the-shelf  spectral 
filters  purchased  all  had  roughly  triangular  spectral  transmissivities,  g^,  which  did 
not  all,  however,  drop  off  to  zero  as  rapidly  as  was  desirable.  The  Integral  in 
equation  (3.1)  must  therefore  be  taken  over  some  whole  or  fractional  number  m^  of  full 
widths  at  half  peak,  AX^,  with  the  number  varying  from  filter  to  filter.  Except  for 
the  two  broad  and  shallow  ultraviolet  filters  centered  on  X^  = 0.258y  and  X2  = 0.359u 
equation  (3.1)  can  be  approximated  with  sufficient  accuracy  in  the  form: 

pinput  (T)  « s®  • n • g^®"®  • I°(X^;  T)  • g^(X^)  • (A'X^/ly)  (3.1a) 


with  A'X^  an  approximation  of  the  bandwidth  that  is  sufficiently  larger  than  AX^  to 
account  for  contributions  from  the  non-zero  "wings"  of  the  filter  function  outside  of 
the  X^  + ^AX^  cut-offs  of  a hypothetical  triangular  filter;  i.e.,  so  that 


TRANSMISSION  PROPERTIES 
OF  LENS  AND  FILTER 


WAVELENGTH  (m) 
(a)  1.9^  FILTER 


WAVELENGTH  (m) 
(b)  2.8p  FILTER 


Figure  3-2.  Transmittance  Curves  for  the  Five  DABS  II-D-2  Filters 
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The  small  square  detector  cell  response  surfaces  of  0.16  cm  area,  A , define  incre- 

c c 2 

mental  "viewing"  solid  angles,  to  = (A  /v  ) . Here  v represents  the  image  distance, 
defined  by  the  geometric  optics  relation: 

1/v  = 1/F  - 1/u  (3.2) 

As  usual,  the  symbols  u and  F in  the  preceding  equation  represent  the  object  distance 
and  the  lens  focal  length.  The  focal  length  F is  equal  to  50.8  cm  (20  inches)  for  all 
five  lenses.  All  these  relationships  are  shown  in  Figure  3-1.  If  the  object  distance 
u is  taken  to  be  equal  to  r,  the  distance  from  the  viewing  port  into  the  gas  volume  to 
the  lenses,  then  v = 60.83  cm  (23.95  inches).  Although  not  necessarily  accurate,  as 
will  be  seen  later,  this  u = r assumption  was  Implicitly  assumed  in  the  earlier  defini- 
tion of  the  solid  angle  0 subtended  from  the  port  aperture  by  the  lenses. 

With  the  v-distances  equal  to  60.83  cm,  the  detector  cell  viewing  solid  angles 

-5  -3 

will  each  subtend  4.324  • 10  steradlans,  with  square  cross  sections  of  6.576  • 10 

c 

radians  on  a side.  As  is  seen  in  Figure  3-1,  the  projections  of  these  to  solid  angles 
onto  the  plane  of  the  viewing  port  are  (slightly  distorted)  squares  that  are  6.576  • 

-3 

10  r,  or  2.027  cm  on  a side.  Neglecting  optical  blurring,  the  diagonals  of  these 
squares  are  equal  to  ^/2  • 2.027  cm  = 2.866  cm,  and  thus  project  only  very  slightly 
beyond  the  2.858  cm  diameter  of  the  port  aperture.  Consequently,  only  very  small  corner 
segments  of  the  to  fields  of  view  are  blocked  off  from  entry  into  the  port  hole. 

The  radiation  to  be  observed,  however,  actually  emanates  from  molecules  within 
the  gas  volume,  rather  than  from  the  plane  of  the  port.  Radiation  originating  from 
within  the  boundaries  of  conical  surfaces  determined  by  the  lens  circumferences  and  the 
port  aperture  have  "clear  shots"  at  the  total  lens  surface,  while  radiation  from  mole- 
cules outside  of  this  conical  surface  have  at  least  some  part  of  the  lenses  obstructed. 
This  is  again  Illustrated  in  Figure  3-1. 

The  ideal  condition,  then,  is  that  the  u)  fields  of  view  not  only  pass  without 
obstruction  through  the  port  opening,  but  also  lie  completely  within  the  unobstructed 
portion  of  the  conical  surface  until  after  some  characteristic  depth  of  penetration  is 
achieved.  It  is  easy  to  show  that,  if  this  ideal  condition  had  been  achieved,  the 
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observed  radiation  signals  then  would  have  been  independent  of  the  depth  distribution 

of  radiating  molecules.  For,  the  solid  angle  subtended  by  the  lens  (f2)  and  the  incre- 

mental  detector  cell  projection  (S  ) terms  occurring  in  equation  (3.1a)  then  would 

2 2 c 2 

have  been  equal  to  nD  /4(r  + 6r)  and  w (r  + 6r)  respectively  for  molecules  at  depth 
6r.  Thus,  for  a sufficiently  large  view  port  aperture  to  allow  unobstructed  lens 
illumination  by  the  gas  molecules  within  the  viewing  solid  angles  u , the  product  of 
these  two  terms  reduce  to 

S®  . 52  - (ttD^/A)  • (3.3) 

Consequently,  equations  (3.1)  and  (3.1a)  would  become  independent  of  r + 6r  and  hence  of 
the  6r  depth  distribution  of  the  observed  radiating  molecules,  and  all  detectors  would 
have  an  optically  unobstructed  view  of  radiating  gases  across  the  width  of  shock 
tunnel  to  the  opposite  wall. 

With  the  5.1  cm  (2  inch)  diameter  lenses  of  51  cm  (20  inch)  focal  length  used,  it 
would  have  required  an  approximately  3.6  cm  (1.4  inch)  diameter  port  opening  to  collect 
unobstructed  radiation  originating  from  molecules  45  cm  deep  into  the  gas  volume.  A 
4.7  cm  (1.6  inch)  diameter  port  opening  would  have  provided  unobstructed  observation 
for  6r  < 200  cm.  The  actual  2.86  cm  (1  1/8  inch)  diameter  port  opening  utilized  was 
originally  designed  for  DABS  II-D-1  and  became  a compromise  design  for  DABS  II-D-2. 

For  DABS  II-D-2  the  smaller  than  optimum  diameter  view  port  will  not  effect  the 
validity  of  substituting  equation  (3.3)  into  equation  (3.1)  or  (3.1a)  for  the  com- 
bustion product  observations,  since  the  high  opacity  of  the  carbon  soot  will  insure 
that  virtually  all  of  the  radiation  to  be  viewed  will  originate  quite  close  to  the  view 
port,  and  moreover,  for  about  the  same  small  mean  6r  for  all  wavelengths.  Since  the 

optical  thickness  x should  be  quite  large  for  the  » 0.26p  shocked  air  radiation, 

equation  (3.3)  also  should  be  a valid  approximation  for  the  radiometer.  The  only 

real  problem  will  be  with  the  interpretation  of  the  X2  0.36)a  radiometer  data.  The 

computed  optical  thickness  x was  relatively  low  for  X=»0.36y,  particularly  at  the  lower 
air  shock  pressures  to  be  observed.  A scaling  of  Figure  E-4  shows  that  the  observed 
radiation  at  this  wavelength  may  originate  from  an  average  of  some  tens  of  centimeters 
within  the  gas  volume  when  the  shock  front  first  fills  the  viewing  port,  and  from  an 
average  of  up  to  ~ 80  to  100  cm  deep  before  arrival  of  the  combustion  products.  It 
thus  will  be  partially  obscured  from  fully  Illuminating  the  lens.  The  (S®  • 52)q 
product  may  therefore  have  to  be  adjusted  downward  from  the  value  given  by  equation  (3.3), 
particularly  as  the  pressure  decreases  from  that  of  the  initial  shock  front  passage. 
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When  equation  (3.3)  is  valid,  a substitution  of  the  numerical  values  of  the  lens 
c ”5  Ifins 

diameter,  D = 5.08  cm,  u = 4.324  • 10  steradians  and  8 =0.9  reduce  the  term 

outside  the  Integral  in  equation  (3.1)  to: 


„g  „ olens  / r.2//N  c -lens  . , -4  2 

S • n • B (ttD  /4)  • u)  • B = 3.105  • 10  cm 


th 

Thus,  the  effective  radiant  power  input  to  the  1 radiometer  becomes 


+ m^(AX^/2) 


(3.4) 


pinput  ^ 2.105  . lo"^  cm^ 


/ 

Xi  - m^(AX^/2) 


(X;  T)  B^(X)  dX 


(3.5) 


for  equation  (3.3)  valid,  and  reduces  to 


pinput  ^2-)  R.  3.105  . 10”^  cm^  . I°(X^;  T)  B^(X^)  • (A'X^/ly) 


(3.5a) 


for  equation  (3.3)  valid  and  the  b5(X)  filter  functions  having  sufficiently  sharp  cut- 

^ -4  2 

off  past  the  filter  bandwidth.  The  constant  3.105  • 10  cm  term  in  equation  (3.5)  or 
equation  (3.5a)  will  cancel  out  when  radiant  power  input,  or  voltage  signal  output 
ratios  are  taken.  Figure  3-2  shows  the  filter  parameters  ^^(X^)  and  A'X^. 

3.1  Radiometer  Output  Signals 

Let  represent  the  i*''^  radiometer  output  signal  in  units  of  either  volts  or 
amperes.  Then, 


Vi  = (T)  • (3.6) 

Here  represents  the  response  of  the  1*"^  radiometer  detector  cell  (and  its  associated 
pre-amplifier)  in  units  of  volts  per  watt  for  the  Infrared  radiometers,  and  of  ampere 
per  watt  for  the  ~0.26M  wavelength  and~0.36vi  wavelength  radiometers:  while  represents 
the.  gain  setting  on  the  corresponding  amplifiers.  It  is  seen  by  comparison  with 
equations  (3.5)  and  (3.5a)  that  when  the  respective  criteria  for  the  validity  of  these 
equations  are  met,  that: 
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(V^/Vj)  « 


I°(A^;  T)  • • (A’X^/lw)  • Rj^  ’ 

I°(Xj;  T)  • Bj(Xj)  • (A’X^/lw)  * Rj  * Gj 


(3.7) 


(3.7a) 


As  in  equations  (3.1)  and  (3.5),  the  integrals  in  equation  (3.7)  are  taken  over  the 

respective  wavelength  Intervals  of  effectively  non-zero  B^(X)  and  B^(X).  Figure  2-2 
f 1 J 

gives  the  B^(X^)  and  A^X,  for  the  five  radiometer  systems  used  in  DABS  II-D-2. 
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A.  ELECTRICAL  SYSTEM  FOR  DABS  II-D-2  TEST 

Once  the  detector  center  frequencies  and  predicted  temperature  range  were  estab- 
lished, the  detection  and  signal  conditioning  circuit  designs  were  determined  as 
follows:  two  distinct  circuit  designs  were  required,  a photovoltaic  silicone  PIN 
diode  design  for  the  two  ultraviolet  detector  channels  and  a photoconductive  lead- 
selenide  design  for  the  three  infrared  detector  channels.  Figure  A-1  shows  the  sche- 
matic of  the  total  circuit  hookup,  from  detector  to  tape  recorder. 

A.l  Photovoltaic  Circuit 

The  optical  signal  from  the  silicone  PIN  diode  is  a single  direction  pulse  whose 
output  is  amplified  by  a single  stage  current  preamplifier.  This  amplified  signal  is 
then  fed  into  a dual  stage  line  driver  amplifier  system  to  transport  the  electrical 
signal  over  the  1500  feet  of  land  line  through  RG58  coaxial  cable  to  the  instrumentation 
trailer,  where  the  tape  recording  equipment  was  located. 

The  tape  recorder  input  preamplifiers  were  set  to  provide  D.C.  offset  so  that  the 
single-direction  pulse  would  use  the  full  scale  minus  to  plus  modulation  of  the  re- 
cording band.  In  addition,  to  insure  high  resolution  results  the  input  signal  was 
split  for  dual  range  recording  of  each  channel.  This  resulted  in  a 20  to  1 difference 

I i 

i;  in  signal  levels  between  the  two  channels.  The  tape  recorder  was  operated  in  the  FM 

[j  mode.  Group  I wideband  at  120  IPS  tape  speed. 

il 

A. 2 Photoconductive  Circuit 

h 

['  The  lead-selenide  circuit  was  designed  with  a single  stage  voltage  preamplifier 

circuit.  The  preamplifier  Integrated  circuit  was  located  as  close  to  the  detector 

as  possible  to  minimize  high  frequency  roll  off  due  to  shunt  capacitance  between  the 

■i  lead-selenide  preamplifier  element  and  the  input  to  the  preamplifier.  Similarly, 

j controls  were  remotely  located  in  the  bottom  of  the  optical  detector  box.  The  atr.pllfied 

•!  signal  was  fed  into  a dual  stage  line  driver  stage.  The  electrical  signal  was  trans- 

; ported  over  1500  feet  of  land  line  to  a recording  set  up  similar  to  that  described  in 

i i 

!i  section  A.l. 

••  l 

[j  The  expected  minimum  and  maximum  signal  level  corresponding  to  estimated  tempera- 

I;  ture  ranges  for  each  channel  is  presented  in  Table  A-1.  The  dynamic  range  of  the 

t electronics  is  presented  in  terms  of  a dynamic  range  ratio  in  the  same  figure. 
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Figure  A-1.  Schematic  of  Optical  Pyrometer  Electronics 
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5.  MECHANICAL  INSTALLATION  FOR  DABS  II-D-2  TEST 

The  primary  purpose  of  the  equipment  pit  Installation  design  is  to  provide  an 
optical  window  into  the  shock  tube  which  will  protect  the  optical  radiometer/electronics 
from  damage,  and  provide  a stable  seismic  base  during  the  5 msec  measurement  envelope. 

A schematic  depicting  the  side  elevation  section  of  the  installation  is  shown  in 
Figure  5-1  (not  to  scale). 

The  main  functional  components  of  the  design  are:  1)  the  wedged  aerodynamic  wing 
section  which  provides  a window  through  the  cooler  boundary  layer  gases  into  the  main 
shock  region,  the  2.86  cm  (1  1/8  inch)  diameter  port  provides  the  main  optical  aperture 
for  viewing  the  shock  wave,  2)  the  4.57  cm  (18  inch)  long  optical  port  provides  a 
uniform  cross  sectional  area  which  allows  pressure  and  temperature  values  to  remain 
constant  for  the  first  0.15  msec  of  the  test.  This  prevents  the  optical  detectors  from 
having  to  look  through  a cooler  gas  volume  during  the  critical  first  stage  compressed 
air  shock,  3)  the  cylindrical  cavity  providing  approximately  a 100  to  1 area  expansion, 
and  along  with  the  three  baffle  plates  provide  a means  to  reduce  the  final  pressure  at 
the  instrumentation  pit  to  a level  compatible  with  equipment  survival,  4)  the  three 
foam  inserts  were  designed  to  accommodate  ground  displacement  and  attenuate  the  induced 
wave  propagation  through  the  soil,  5)  the  seismic  mass  was  designed  to  provide  an 
adjustable  stable  platform  for  placement  of  optics. 

The  pit  was  lined  with  a thin  plastic  sheet  to  keep  moisture  and  dirt  from  the 
optical  box  prior  to  test.  The  addition  of  the  lid  wrapped  with  excess  plastic  sheet 
provided  the  system  with  protection  from  the  elements  during  the  check-out  and  cali- 
bration periods  prior  to  the  test.  Finally  a sand  overburden  was  placed  over  the  lid 
and  baffle  area  to  protect  the  pit  from  debris  during  the  test. 

Figure  5-2  is  a series  of  photographs  showing  the  actual  instrumentation  used  to 
obtain  the  temperature  history  test  data.  Photograph  (a)  shows  the  complete  instrumenta- 
tion package,  the  black  optical  box  on  the  right  contains  the  lenses,  filters,  detectors 
and  detector  amplifier  stages.  The  package  on  the  right  contains  the  power  supply  and 
line  driver  amplifiers. 

Photograph  (b)  shows  the  interior  of  the  optical  box.  The  front  panel  contains  five 
lens  holders  which  can  be  adjusted  for  various  focal  lengths  to  maximize  the  signal  at  each 
detector.  The  rear  panel  has  the  five  detector/f liter  assemblies  mounted  so  that  they 
remain  parallel  to  the  front  len  plane.  The  bottom  of  the  optical  box  contains  the 
detector  amplifiers  and  associated  electronics.  Photograph  (c)  provides  a clear  view  of 
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(b)  INTERIOR  OF  OPTICAL  BOX  WITH  DETECTORS  AND  FILTERS  ON  LEFT 
LENS  ON  RIGHT  AND  ELECTRONICS  ON  THE  BOTTOM 


Photographs  of  Instrumentation  used  to  Measure 
Temperature  History  Data 


Figure  5-2 


(c)  DETECTOR  ASSEMBLIES  WITH  TWO  OF  THE  FILTERS  REMOVED 


Figure  5-2  Photographs  of  Instrumentation  used  to  Measure 
Temperature  History  Data  (Concluded) 


(d)  FIELD  INSTALLIZATION  SHOWING  EQUIPMENT  IN 
INSTRUMENTATION  PIT  AFTER  TEST 


I 


the  detector/filter  assembly.  Two  of  the  assemblies  have  their  filters  removed  to 
expose  the  detectors  inside.  The  upper  detector  is  a lead-selenlde  photoconductive 
detector,  while  the  bottom  unit  shows  a silicone  PIN  diode  photovoltaic'  detector. 

The  last  photograph  (d)  shows  the  instrumentation  pit  at  the  test  site.  This 
photograph  was  taken  after  the  test.  The  complete  equipment  package  survived  the  test 
shock  with  the  only  damage  occuring  to  the  land  lines.  These  coaxial  lines  were  pulled 
out  from  the  connectors  on  the  rear  instrumentation  package,  as  clearly  shown  in  the 
lower  left  of  the  photograph. 
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6.  ANALYSES  OF  THE  DABS  II-D-2  TEST  DATA 


The  optical  temperature  data  generated  by  the  DABS  II-D-2  test  consisted  of  two 
ultraviolet  radiometer  channels  and  three  infrared  radiometer  channels.  As  mentioned 
earlier,  the  data  obtained  from  the  three  Infrared  channels  could  not  be  interpreted 
due  to  spurious  noise  interference  during  the  test.  These  signals  saturated  the 
electronics  and  rendered  the  infrared  data  meaningless.  The  results  presented  in 
sub-section  6.1  will  consist  of  an  analyses  and  Interpretation  of  the  ultraviolet 
channel  data,  while  sub-section  6.2  will  present  a discussion  of  some  secondary,  non- 
temperature related  inferences  obtainable  from  the  temperature  history  information. 

6.1  Primary  Results:  Temperature  Histories 

The  temperature  history  results  for  the  DABS  II-D-2  test  are  clearly  shown  in 
Figure  6-1.  The  data  points  shown  on  this  figure  were  obtained  by  reading  voltage 
values  at  selected  time  intervals  from  oscilloscope  traces  which  had  been  reproduced 
from  the  master  data  tape.  Figure  6-2  displays  typical  traces  used  to  obtain  these 
voltage-time  data  points.  Next,  these  voltage  readings  were  converted  to  temperature 
measurements  by  reading  the  calibration  values  off  Figure  6-3  for  the  two  color  pyro- 
meter temperatures  and  Figure  6-4  for  one  color  pyrometer  temperatures.  All  the  0.26p 
and  0.36p  ultraviolet  wavelength  temperature  data  shown  in  Figure  6-1  were  obtained  in 
this  manner. 

The  Information  presented  in  Figure  6-1  can  be  divided  into  five  areas  of  interest. 

• Precursor  Interval  (t  < 0) 

• Air  shock  arrival  at  test  station  (t  = 0) 

• Compressed  air  shock  interval  (0^t<0.8  msec) 

• Contact  surface  arrival  at  test  station  (t  » 0.8  msec) 

• Combustion  products  Interval  (t  > 0.8  msec) 

In  addition,  theoretical  temperature  envelopes  representing  the  Rankine-Hugonlot  relation- 
ship, with  Isentropic  expansion  of  the  shock  compressed  air  are  presented  for  comparison. 
The  upper  limit  is  defined  by  pre-test  predictions  from  peak  pressure  of  83.3  atmospheres 
(1225  psia)  at  the  8.2  m (27  ft)  station.  The  lower  bound  was  obtained  from  AFWL  post- 
test data  [l]  which  presented  a average  peak  pressure  of  67.2  atmospheres  (987  psia)  at 
the  same  station.  Both  theoretical  curves  represent  a constant  gas  constant  of  y “ 1.3. 
For  convenience  the  arrival  of  the  peak  pressure  at  the  station  is  used  as  the  t = 0 
time  reference  point. 

The  precursor  region  corresponds  to  forward  scattered  radiation  arriving  before 
the  shock  front  enters  the  radiometer  fields  of  view.  A short  and  fairly  reliable 
extrapolation  of  the  forward  scattered  precursor  temperature  signal  indicates  that  the 
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(a)  0.36  M CHANNEL 


(b)  0.26  M CHANNEL 


TRACES  (a)  and  (b):  SHOCK  FRONT  ARRIVAL  THROUGH  SIGNAL  DECAY  INTO  NOISE. 

VERTICAL  SCALES  ARE  0.5  VOLTS/CM  AND  HORIZONTAL  SCALES 
ARE  312.5  fi  SEC/CM,  REAL  TIME. 


Figure  6-2.  DABS  II-D-2  Ultraviolet  Radiometer  Signal  Tracings 


(c)  0.36  H CHANNEL 


(d)  0.26  M CHANNEL 


TRACES  (c)  AND  (d):  SHOCK  FRONT  ARRIVAL  THROUGH  OPTICAL  SIGNAL  TAIL-OFF 
OWING  TO  ARRIVAL  OF  COMBUSTION  PRODUCTS.  VERTICAL 
SCALES  ARE  0.5  VOLT/CM  AND  THE  HORIZONTAL  SCALES 
ARE  125  u SEC/CM,  REAL  TIME. 


Figure  6-2.  DABS  II-D-2  Ultraviolet  Radiometer  Signal  Tracings  (Continued) 
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(e)  0.36  ^J  CHANNEL 


(f)  0.26  H CHANNEL 


TRACES  (e)  AND  (f):  TIME  SCALE  EXPANSION  OF  EARLY  PORTIONS  OF  TRACES  (c) 

AND  (d)  RESPECTIVELY.  VERTICAL  SCALES  ARE  0.5  VOLT/CM 
AND  HORIZONTAL  SCALES  ARE  62.5  h SEC/CM,  REAL  TIME. 


Figure  6-2.  DABS  ll-D-2  Ultraviolet  Radiometer  Signal  Tracings  (Continued) 
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(h)  0.26  H CHANNEL 


TRACES  (g)  AND  (h): 


A STILL  FURTHER  TIME  SCALE  EXPANSION  OF  VERY  EARLY 
PORTIONS  OF  TRACES  (c)  AND  (d)  RESPECTIVELY.  VERTICAL 
SCALES  ARE  0.5  VOLTS/CM  AND  THE  HORIZONTAL  SCALES 
ARE  31.25  H SEC/CM,  REAL  TIME. 


Figure  6-2.  DABS  II-D-2  Ultraviolet  Radiometer  Signal  Tracings  (Continued) 
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(i)  0.36  H CHANNEL 


(i)  0.26  H CHANNEL 


TRACES  (i)  AND  (j):  EXPANDED  SCALE  OPTICAL  RADIATION  PRECURSOR  SIGNALS  OWING 

TO  FORWARD  SCATTERED  RADIATION  BEFORE  ARRIVAL  OF  THE  SHOCK 
FRONT  AT  THE  OBSERVING  STATION  ( t NEGATIVE,  SINCE  ARBITRARILY 
SET  EQUAL  TO  ZERO  AT  TIME  OF  PEAK  PULSE  ARRIVAL).  VERTICAL 
SCALE  IS  0.1  VOLT/CM  AND  THE  HORIZONTAL  SCALE  IS  312.5  p SEC/ 
CM,  REAL  TIME. 


Figure  6-2. 


DABS  Il-D-2  Ultraviolet  Radiometer  Signal  Tracings  (Continued) 


(I)  0.26  M CHANNEL 


TRACES  (k)  AND  (I):  RADIATION  SIGNAL  FALL-OFFS  SHOWING  ARRIVAL  OF 

COMBUSTION  PRODUCTS.  VERTICAL  SCALE  IS  0.1  VOLT/ 
CM  AND  HORIZONTAL  SCALE  IS  125  u SEC/CM,  REAL  TIME. 


Figure  6-2.  DABS  II-D-2  Ultraviolet  Radiometer  Signal  Tracings  (Concluded) 
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Figure  6-3.  Two  Color  Pyrometer  Calibration  Curve  for  the  Ultraviolet  Radiometers 
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mean  temperature  averaged  over  the  entire  shock  front  was  of  the  order  of  ~4600  °K  + 

250  °K.  At  t = 0,  corresponding  to  the  arrival  of  the  peak  air  blast  pressure,  the  two 
color  pyrometer  temperature  (near  port)  Is  4080  °K  with  random  error  of  the  order  + 50  °K. 
As  the  shock  front  propagates  downstream,  the  mean  temperature  of  the  viewed  air  volume  in 
front  of  the  port  aperture  decreases.  Shortly  before  the  arrival  of  the  contact  surface, 
at  approximately  t = 0.78  msec  the  mean  temperature  In  the  air  volume  as  determined  from 
the  one  color  0.26y  wavelength  pyrometer  is  3410  °K  within  a random  error  of  + 100  °K. 
Similarly,  the  one  color  0.36y  wavelength  pyrometer  recorded  a mean  temperature  of  3150  °K 
(corrected),  with  a random  error  of  + 80  °K,  The  above  random  error  data  was  measured 
from  recorded  signal  noise  data  that  appears  on  the  photographs  in  Figure  6-2. 

The  ultraviolet  channels  are  incapable  of  measuring  the  low  combustion  product 
temperatures.  Since  the  infrared  data  was  lost,  a mean  absolute  temperature  for  combustion 
product  gases  cannot  be  determined.  The  one  color  ultraviolet  pyrometer  readings  shows  a 
rapid  drop  off  in  signal  level  in  the  0.85  to  1.0  msec  range,  thus  indicating  the  arrival 
of  the  cooler  combustion  product  gases.  This  value  is  within  0.05  msec  of  the  predicted 
theoretical  value. 

The  temperature  curves  exhibit  the  expected  exponential  decay  between  peak  pressure 
and  contact  surface  arrival  times.  The  data  is  in  good  agreement  with  the  theoretical 
values.  Moreover,  Figure  6-5  present  data  that  graphically  displays  the  sensitivity  of 
these  theoretical  curves  to  changes  in  gas  constant  (y)  and  pressure.  These  results 
show  that  the  use  of  y = 1.3  as  a mean  gas  constant  is  supported  by  these  test  results. 

A more  detailed  discussion  of  all  the  test  results  will  be  presented  in  the  fol- 
lowing sub-section. 

6.1.1  Detailed  Discussion  of  DABS  II-D-2  Results 

As  mentioned  in  the  preceding  section,  all  the  two  color  pyrometer  temperature 
measurements  used  the  Figure  6-3  theoretical  voltage  signal  ratio  versus  temperature 
calibration  curve  for  the  0.26y  and  0.36u  wavelength  radiometer  pair.  It  is  similar 
to  the  1^/1^  versus  temperature  curve  in  Figure  A-4  except  that  it  applies  directly  to 
the  measured  voltage  rather  than  spectral  radiance  ratios.  The  actual  voltage  ratio 
was  computed  by  numerical  convolutions  from  equation  (3.7),  since  the  transmission  of 
the  two  spectral  filters  and  particularly  the  0. 26 y wavelength  bandpass  filter,  decayed 
too  slowly  to  use  the  approximation  of  equation  (3.7a). 

The  usable  two  color  pyrometer  data  was  limited  to  the  precursor  region  (t < 0), 
the  peak  arrival  time  (t  = 0)  and  approximately  the  first  20  ysec  of  the  compressed 
air  shock  front.  The  accuracy  of  the  two  color  pyrometer  starts  to  decrease  for 
t > 20  ysec  as  will  be  discussed  later.  The  two  color  pyrometer  data  is  represenced 
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in  Figure  6-1  by  circles  around  each  data  point  corresponding  to  (V^  36  ^ 

voltage  ratios;  error  estimates  are  represented  by  the  vertical  bars. 


The  plotted  points  In  the  precursor  region  (negative  t)  indicates  a large  variance 
in  measured  real  temperature.  The  gradual  rise  of  the  two  voltage  signals  above  their 
noise  level  (see  Figures  6-11  and  6-lj)  is  almost  certainly  due  to  forward  scattered 
radiation  reflected  off  t)ie  chamber  walls  and  into  the  radiometer  fields  of  view.  The 
emitting  layer  is  probably  a little  deeper  for  the  0.36)j  wavelength  than  for  the  0.26p 
wavelength  radiation,  but  this  difference  in  optical  depth  should  be  appreciably  less 
than  for  the  directly  viewed  radiation.  Contributions  to  this  forward  scattered  radi- 
ation should  arrive  from  virtually  the  entire  shock  front  rather  than  from  the 
restricted  local  areas  within  the  radiometer  fields  of  view  as  is  true  for  the  directly 
viewed  radiations.  The  large  error  ranges  indicated  arise  from  the  low  signal  to  noise 
ratios  and,  as  a systematic  error,  from  the  difficulty  in  assessing  "zero  signal  levels" 
from  which  the  voltage  is  measured.  The  failure  of  the  indicated  temperature  to  decrease 
slightly  with  time  probably  is  due  to  measurement  errors.  It  is  unlikely  that  the  several 
hundred  degree  differential  between  the  indicated  forward  scattered  radiation  tempera- 
ture "points"  and  the  peak  response  temperature,  at  t = 0,  is  due  to  systematic  errors. 

It  is  more  probable  that  the  weight  mean  temperature  over  the  entire  shock  front  given 
by  the  forward  scattered  radiation  is  higher  than  the  temperature  of  the  gas  near  the 
port  aperture.  It  also  is  possible  that  the  diffuse  reflectivity  of  the  inner  walls  of  the 
shock  chamber  could  be  slightly  greater  near  0.26y  than  near  0.36p,  thus  increasing  the 
ratio  of  the  0.26p  signal  to  the  0.36p  signal  in  the  precursor  region  so  as  to  give  the 
appearance  of  higher  temperature  than  actually  exists.  It  is  difficult  to  obtain  ultravio- 
let reflectivity  coefficients  in  the  literature  to  check  this  latter  possibility. 

The  two  theoretical  envelopes  indicate  a dropping  temperature  after  passage  of 
the  shock  front  past  the  observation  station.  The  voltage  signal  tracing  reproduced 
in  Figures  6-la  through  6-11  exhibit  decays  consistent  with  dropping  radiant  power 
inputs,  and  hence  dropping  temperature.  The  plotted  two  color  pyrometer  data  points 
(open  circles  in  Figure  6-1),  however,  appear  to  show  somewhat  random  oscillations  about 
a rising  rather  than  falling  temperature  pattern.  This  apparent  inconsistency  readily 
was  explained  by  further  analysis. 

As  was  previously  mentioned,  shock  heated  air  exhibits  appreciably  weaker  absorption 
coefficients,  and  correspondingly  deeper  optical  penetration,  at  0.36y  than  at  0.26y. 

This  would  not  have  mattered  if  the  necessary  assumption  of  approximately  uniform 
temperature  in  gas  layers  parallel  to  a presumably  planar  shock  front  were  valid.  (This 
assumption  probably  would  have  been  more  valid  for  the  original  DABS  lI-D-1  con- 
figuration than  for  the  substitute  DABS  II-D-2  configuration  with  its  appreciably 
smaller  downstream  range  to  chamber  width  ratio).  It  is  clear,  however,  that  the 
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0.26u  wavelength  radiometer  is  more  sensitive  to  temperature  variations  occurring  very 
close  to  the  viewing  port,  while  the  0,36p  wavelength  radiometer  will  respond  more 
to  temperature  variations  deeper  into  the  shock  chamber. 

Very  soon  after  the  passage  of  the  main  shock  front  across  the  port  aperture,  a 
secondary  perpendicular  shock  front  is  propagated  down  the  45.7  cm  (18  inch)  pipe 
projecting  normally  from  the  viewing  port, (see  Figure  5-1).  It  was  Initially  expected 
that  the  air  in  this  pipe  would  not  differ  appreciably  in  temperature,  and  if  anything 
be  slightly  cooler,  than  the  air  in  the  main  chamber  adjacent  to  the  viewing  port. 
Further  analyses  revealed,  however,  that  while  this  would  be  true  for  most  of  the  pipe 
length,  the  diverted  air  will  be  heated  by  reflected  shocks  off  the  pipe  inner  surface 
during  the  first  very  few  inches  of  propagation.  This  few  Inches  is  sufficient  to 
strongly  influence  the  radiation  signal  to  the  0.26y  wavelength  radiometer,  but  would 
make  relatively  negligible  contribution  to  the  0.36u  wavelength  radiometer  signal. 

It  thus  appears  that  only  those  two  color  pyrometer  data  points  plotted  in 
Figure  6-1  for  t < 20  usee  correspond  to  actual  temperatures.  The  other  two  color 
pyrometer  data  points  correspond  merely’  to  ratios  of  received  spectral  radiance  eman- 
ating from  overlapping  gas  volumes  of  different  mean  temperatures. 

Insufficient  empirical  data  exists  for  the  accurate  absolute  calibration  of  the 
0.26u  and  0. 36 u wavelength  radiometers  as  "one  color  pyrometers".  Moreover,  it  is 
very  difficult  and  expensive  to  obtain  reliable  black  body  (or  pseudo  black  body) 
sources  for  ultraviolet  wavelengths  at  temperatures  above  ~3000  °K  to  perform  this 
absolute  calibration.  However,  fairly  reliable  relative  one  color  pyrometer  calibration 
can  be  made,  because  the  two  color  pyrometer  data  points  for  0 < t < 20u  sec  shown  in 
Figure  6-1  represent  meaningful  mean  temperature  for  the  gas  volume  inside  the  main 
chambers  and  in  the  near  vicinity  of  the  viewing  port.  The  numerator  and  denominator 
of  equation  (3-7)  were  previously  computed  for  a selected  set  of  temperatures  in  order 
to  compute  the  two  color  pyrometer  ratio  curve  shown  in  Figure  6-3.  Theoretically, 
each  of  these  quantities  is  proportional  to  the  corresponding  0.26p  and  0.36y  wavelength 
radiometer  output  voltages  for  the  corresponding  temperature.  The  unknown  constants 
of  proportionality  are  eliminated  by  dividing  each  set  of  computed  values  through  by 
its  corresponding  value  for  T = 3500  °K.  The  resulting  one  color  pyrometer  data 
evaluation  curves  are  shown  plotted  in  Figure  6-4.  They  are  separately  normalized  to 
an  arbitrary  pseudo  voltage  signal  of  unity  at  T = 3500  “K,  but  are  not  normalized 
with  respect  to  each  other. 

The  best  two  color  pyrometer  estimate  of  temperature,  evaluated  from  Figure  6-1, 
is  T = 4080  °K  at  t = 0 (peak  pulse).  It  is  seen  that  for  4080  °K,  tbe  normalized  one 
color  pyrometer  curves  of  Figure  6-4  give  (pseudo)  voltage  readings  of  5.1  "volts" 


for  the  0.361J  radiometer  and  9.3  "volts"  for  the  0.26u  radiometer.  When  expanded  to  j 

the  same  arbitrary  vertical  scale  the  signal  voltage  output  tracings  gave  t = 0 readings  j 

of  51.0  for  the  0.36p  wavelength  radiometer  and  an  average  of  28.8  for  the  0.26p  wave-  i 

I 

length  radiometer.  The  tracing  readings  for  any  later  time  t were  first  normalized  to  i 

these  respective  t = 0 values  before  entering  the  corresponding  one  color  pyrometer  i 

curve  in  Figure  6-4.  This  was  done  by  multiplying  the  tracings  vertical  coordinate  | 

i 

readings  by  (5.1/51.0)  = 0.10  for  the  0.36ij  wavelength  radiometer,  and  by  (9.3/28.8)  = i 

0. 323  for  the  0.26p  wavelength  radiometer.  The  two  resulting  computed  one  color  pyro-  i 

' ) 

meter  temperature  histories  also  are  plotted  in  Figure  6-4.  It  is  seen  that  they  give 
reasonable  appearing  temperature  decays,  with  the  apparent  mean  temperatures  differential  j 

between  the  two  radiometers  gradually  Increasing  to  ~360  °K  near  0.8  msec  before  the  i 

arrival  of  the  combustion  gases.  | 

I 

t 

There  possibly  may  be  some  fall  off  of  emissivity  below  ~1  for  the  received  0.36y  i 

wavelength  radiometer  radiation  signal  near  0.8  msec.  The  principal  non-temperature 

related  signal  decay  mechanism,  however,  is  the  partial  occlussion  of  the  radiation  ! 

collecting  lens  by  the  edges  of  the  port  aperture.  From  the  discussion  in  Section  3, 
it  follows  that  as  temperature  and  pressure  decrease,  and  the  relative  radiation  con- 
tributions from  molecules  progressively  deeper  behind  the  port  aperture  increase,  this 
partial  occlussion  will  cause  increasing  fractional  signal  loss  in  the  0.36y  wavelength  | 

channel.  It  is  quite  possible  that  this  occlussion  related  signal  loss,  which  does  i 

not  effect  the  0.26p  wavelength  channel,  is  responsible  for  some  fraction  of 

the  apparent  temper ture  differential  Increases  between  the  two  channels.  i 

6.1.2  Comparison  of  Measured  with  Theoretical  Temperatures 

The  theoretical  temperature  pressure  curves  presented  in  Figures  1-1  through 
1-4,  and  the  data  on  Table  1.1  were  based  the  somewhat  Idealized  and  simplified  Ranklne- 
Hugonlot  relationships  for  shock  front  parameters  (P^j  upon  similarly  simplified 

Isentroplc  decay  theory  for  the  relationship  of  temperature  to  pressure  decay  ratios 
behind  the  shock  front  at  a given  observation  station.  These  relationships  are  given 
by  equation  (1.1)  and  (1.2)  in  section  1.0.  The  theories  behind  both  these  simplified 
equations  utilize  the  convenient  and  empirically  justified  fiction  that  y,  the  ratio 
of  Cp  the  specific  heat  at  constant  pressure  to  C^,  that  at  constant  volume  (i.e.,  | 

Y = Cp/C^),  remains  constant  during  propagation  of  the  shock  front  and  also  during  the  i 

assumed  isentroplc  expansion  and  cooling  of  the  air  behind  the  shock  front  [2].  Moreover,  | 

the  figures  and  theoretically  computed  data  presented  in  Section  1.3  used  the  I 

additional  first  order  assumption  that  y was  equal  to  its  cold  air  value  of  1.4, 

1. e,,  the  Y“Value  in  front  of  the  shock.  Also,  from  previous  empirical  observa- 
tions it  was  assumed  that  the  contact  surf act,  between  the  shocked  air  layer  and  the 
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combustion  gases,  would  arrive  at  approximately  half  shock  front  pressure,  l.e.,  when 
(P^/Pg)  ~ 0.5.  From  Figure  1-3  and  Table  1.1  It  Is  seen  that  the  contact  surface  (or 
P^  = 0.5  P^)  was  predicted  to  pass  the  DABS  lI-D-2  observation  station  ~0.8  msec  after 
passage  of  the  shock  front.  It  can  be  seen  In  Figure  6-1  that  this  6t  prediction  was 
quite  good. 

It  Is  quite  obvious  from  theoretical  considerations  that  y really  cannot  and  does 
not  remain  constant  either  on  the  propagating  shock  front  or  behind  It.  The  partition 
functions  of  the  constituent  air  molecules  Increase  markedly  with  Increasing  tempera- 
tures, and  with  them,  the  Internal  energy  contributions  of  these  molecules  to  C^. 

Thus,  must  Increase  with  temperature,  and  since  C “ + R,  the  value  of  y must 

correspondingly  tend  to  decrease  toward  unity  for  very  high  temperatures.  Moreover, 
the  shocked  air  chemistry,  as  Is  shown  In  Figures  D-1  and  E-1,  alters  the  mole  fractions 
of  constituent  air  species  molecules,  and  thxjs  alters  the  y values  as  a function  of 
pressure  and  temperature. 

A higher  order  approximation  to  shocked  air  temperature  pressure  relations  Is  to 
still  use  equations  (1.1)  and  (1.2),  but  to  substitute  suitably  selected  effective  mean 
values  of  y for  the  first  order  y = l.A  approximation.  A 1970  DASA  report  by  H.  L.  Erode 
[2]  presents  data  from  which  such  effective  mean  shocked  air  y values  can  be  extracted 
for  varying  shock  front  pressures,  P^.  Assuming  these  y values  to  be  approximately 
correct,  and  also  that  the  shock  front  pressures  P^  are  known,  corresponding  shock  front 
temperatures,  T^,  can  be  computed  from  equation  (1-1).  Then  Isentroplc  decay  tempera- 
tures, T^,  can  be  computed  for  selected  ratios  from  equation  (1.2).  This  was 

done  to  generate  the  theoretical  temperature  curves  shown  In  Figures  6-1  and  6-5, 
using  the  previously  mentioned  additional  assumption  that  ~ 0-5  at  t ~ 0.8  msec 

after  the  shock  front  passage  by  the  observation  station. 

It  Is  seen  In  Section  1.3  that  the  predicted  shock  front  pressure  for  the  DABS 
II-D-2  observation  station  was  P^  ~ 83.3  atmospheres  (or  1225  psla) . Allowing  for 
pressure  gage  "overshoot",  this  value  Is  reasonably  consistent  with  pressure 
versus  time  charts  sent  to  TRW  by  AFWL  [ l] . These  pressure  plots  were  measured 
at  two  stations,  within  7.6  cm  (3  Inches)  of  one  another,  and  9.1  meters 
(29.86  ft)  downstream  from  the  driver  wall.  That  Is  ~.87  m (2.86  ft)  downstream  of 
the  radiometer  station  (but  on  the  opposite  wall).  The  peak  excursions  (t  = 0 values) 
of  these  two  traces  corresponded  to  pressures  of  ~88  atmospheres  (—  1300  psla) 
and  ~102  atmospheres  (~1500  psla),  but  were  subject  to  reduction  for  "overshoot",  so 
that  an  average  value  of  the  order  of  ~83  atmospheres  (1225  psla)  particularly  at  a 
station  ~. 87  m (~2.86  feet)  upstream.  Is  not  unreasonable.  However,  a more  recent 
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telephone  conversation  with  Mr,  Joe  Renick  of  AFWL  indicated  that  shock  front  velocity 

data  could  be  interpreted  to  imply  a pressure  of  ~67  atmospheres  (*^87  psla)  near  the 

TRW  radiometer  station.  The  Brode  report  [2]  (Figure  8,  page  17)  implies  effective  mean  Y 

values  of ~ 1.30  for  P ~ 83  atmospheres  (1225  psla)  and~1.32  for  P 67  atmospheres 

s s 

(~990  psla).  The  corresponding  T values  computed  from  equation  (1.1)  are  A224  ®K 
and  3641  °K  respectively.  For  P^  = 67.2  atmospheres  (987  psla)  and  y = 1.30,  the  T^ 
value  predicted  by  equation  (6-1)  is  3460  ®K. 

The  constant  y approximation  theoretical  temperature  curves  shown  in  Figures  6-1 
and  6-5  are  merely  straight  lines  connecting  T^  points  computed  from  equation  (1.1) 
with  corresponding  T^  points  computed  from  equation  (1.2)  for  equal  to  0.05 

and  positioned  at  the  predicted  station  time  Increment  of  0.8  msec  after  shock  front 
arrival  at  t = 0.  The  two  theoretical  curves  shown  in  Figure  6-1  are  both  for  y = 1.3, 
with  the  upper  curve  corresponding  to  P^  = 83.3  atmospheres  (1225  psia)  and  the  .ower 
to  P^  = (67.2  atmospheres  (987  psia).  It  is  seen  that  these  two  curves  bracket  the 
experimental  data.  They  also  serve  to  show  the  variation  in  constant  y approximation 
theoretical  temperature  predictions  corresponding  to  a change  of  assumed  shock  front 
pressure.  The  corresponding  y = 1.32,  P^  = 67.2  atmospheres  (987  psia)  curve  is  not 
shown,  since  its  variation  from  the  y = 1.30,  P^  = 67.2  atmospheres  (987  psia)  curve 
is  relatively  small  compared  to  the  limited  degree  of  validity  of  the  constant  y theory 
and  to  the  uncertainties  in  pressure  determination. 

Figure  6-5  shows  six  constant  y theoretical  curves,  for  y = 1.4,  1.3,  and  1.2, 
and  for  P^  = 83.3  atmospheres  (1225  psla)  and  67.2  atmospheres  (987  psla)  for  each  of 
these  y values.  The  y = 1.4  value  was  chosen  because  the  y value  for  the  air  in  front  of 
the  shock  often  is  used  for  first  order  temperature  predictions  (i.e.,  in  Section  1.3). 
The  1.2  value  corresponds  to  the  equilibrium  y value  computed  by  an  in-house  TRW 
program  for  air  at  4000  °K  and  80  atmospheres  (1176  psia)  pressure.  This  set  of  curves 
serve  to  show  the  effects  of  variations  of  both  shock  front  pressure  and  y value 
assumptions  on  predicted  temperatures,  and  to  emphasize  the  necessity  for  making  inde- 
pendent empirical  temperature  measurements. 

6.2  Secondary  Results 

Inspection  and  analyses  of  the  two  ultraviolet  radiometer  signal  tracings  reveal 
some  information  in  addition  to  the  computed  temperature  data  presented  in  the  preceding 
section. 

6.2.1  Mean  Shock  Front  Velocity 

Figure  6-6a  and  b,  for  the  0.36iJ  and  0.26vi  wavelength  radiometers  respectively, 
show  the  delay  time  of  these  optical  signals  with  respect  to  the  fiduciary  pulse  which 
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(b)  0.26  H RADIOMETER  SIGNAL 

TRACES  (a)  AND  (b):  THE  VERTICAL  SCALE  IS  0.5  VOLT/CM  AND  THE 

HORIZONTAL  SCALE  IS  312.5  |a  SEC/CM,  REAL  TIME. 


Figure  6-6.  Delay  of  DABS  Il-D-2  Optical  Signals  With  Respect  to  Fiduciary  Pulse 


started  the  oscilloscope  sweep.  This  fiduciary  pulse  also  corresponds  to  the  actual 
field  electrical  pulse  which  served  to  ignite  the  explosive  that  acted  as  the  "shock 
piston".  The  vertical  scales  for  these  tracings  are  0.5  volts/cm,  while  the  horizontal 
time  scale  is  312. 5u  sec/cm,  real  time.  The  time  from  the  power  pulse  to  initiation 
of  the  main  optical  signal  pulses  is  2.875  msec.  Information  received  from  AFWL 
personnel  indicated  a 1.2  msec  delay  between  the  ignition  pulse  and  the  beginning  of 
shock  front  propagation.  If  this  1.2  msec  is  subtracted  from  the  2.875  msec  optical 
signal  delay  time,  then  the  computed  mean  velocity  to  the  8.2  m (27  ft)  downstream 
station  is  ~ (8.23m/1.68  • 10  ^ sec)  = ~4.9  • lO^m/sec  ( 1.6  • 10^  ft/sec). 

The  preceding  value  of  mean  velocity  during  the  first  8.2  meters  (27  ft)  of  shock 

front  propagation  appears  to  be  theoretically  correct.  The  speed  of  ignition  in  the 

4 4 

high  explosive  material  itself  should  be  of  the  order  of  6.6  • 10  m/sec  (2  • 10  ft/sec), 
and  the  indicated  initial  pressure  of  the  combustion  gases  acting  as  the  shock  driver 

is  of  the  order  of  24,000  atmospheres.  The  initial  shock  front  velocity  of  this  order 

3 3 4 4 

of  5.8  • 10  to  8.5*  10  m/sec  (1.9  • 10  to  2.8  • 10  ft/sec)  are  not  unreasonable. 

Also,  available  theoretical  curves  appear  to  indicate  an  instantaneous  shock  front 

4 4 

propagation  speed  of  the  order  of  3.3  • 10  m/sec  (1  • 10  ft/sec)  for  a pressure  of 
~81.6  stmospheres  (~1200  psla)  and  temperature  of  ~4000  “K. 

6.2.2  Lag  Time  Between  the  Two  Ultraviolet  Radiometers 

As  has  been  mentioned  previously,  the  fields  of  view  of  the  0.36y  and  0.26y 
wavelength  radiometers  are  separated  laterally  by  15.2  cm  (6  inches)  and  produce  a 
parallax  of~0.0206  radians  upstream  and  downstream  of  the  approaching  shock  front 
(See  Figure  2-2).  It  was  expected  that  the  0.36m  wavelength  radiometer  would  detect 
the  approaching  shock  front  slightly  before  the  0. 26u  wavelength  radiometer.  The  actual 
relative  delay  was  found  to  be~7.5  ysec.  This  is  shown  in  the  expanded  tracing  segments 
presented  as  Figure  6-7  a,  b.  The  time  scale  for  these  two  segments  is  12.5m  sec/ cm, 
real  time  [0.2  msec/cm  oscilloscope  sweep],  the  signal  amplitudes  are  not  to  scale. 

4 

If  the  theoretically  indicated  shock  front  propagation  speed  of  ~3.3  • 10  m/sec 

4 

(~1  • 10  ft/sec)  is  assumed,  the  7.5Msec  relative  delay  time  indicates  that  the  shock 
front  entered  the  leading  edge  of  the  0.36m  wavelength  radiometer  field  of  view  ~ 2. 3 cm 
(~0.9  inches)  upstream  of  where  it  entered  that  of  the  0.26m  wavelength  radiometer.  If 
the  assumption  were  made  that  the  shape  of  the  shock  front  were  such  that  the  depths 
of  penetration  of  these  two  fields  of  view  into  the  shock  chamber  were  the  same  for 
the  initial  detections,  then  a simple  computation  implies  that  this  hypothetical  first 
contact  depth  was  *-55.9  cm  ( ~ 22  inches).  The  virtually  certain  irregularities  in  the 
shock  front,  however.  Increase  the  uncertainty  associated  with  the  above  55.9  cm  (22 
inch)  Initial  contact  depth  value. 
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THE  HORIZONTAL  SCALE  FOR  BOTH  TRACES  IS  12.5  SEC/CM, 
REAL  TIME,  (0.2  MSEC/CM  OSCILLOSCOPE  SWEEP) 


LOWER  TRACE: 
UPPER  TRACE: 


THE  0.36  H RADIOMETER  SIGNAL  AMPLITUDE  HAS  A 
VERTICAL  SCALE  OF  1 VOLT/CM 

THE  0.26  H RADIOMETER  SIGNAL  AMPLITUDE  HAS  A 
VERTICAL  SCALE  OF  0.5  VOLT/CM 


Figure  6-7.  Relative  Lag  Time  Between  Arrivals  of  the  Two  DABS  lI-D-2  Ultraviolet 

Radiometer  Signals 
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While  some  of  the  semi-oscillatory  structure  seen  in  Figure  6-7  probably  can  be 
attributed  to  noise  or  to  instrument  imperfections,  much  of  it  undoubtedly  is  real. 

It  apparently  indicates  spatial,  temporal  and  thermal  irregularities  in  the  propagating 
shock  front.  In  particular,  the  pronounced  hump  occurring  during  the  rise  of  the  O.Sbp 
wavelength  pulse,  and  some  of  the  oscillations  past  this  peak  value,  do  not  look  at 
all  like  the  passage  of  a thin  uniform  planar  shock  front  oriented  normal  to  the  tunnel 
centerline.  It  will  be  remembered  that  the  0.36p  wavelength  radiometer  receives 
radiation  from  regions  appreciably  deeper  into  the  shock  chamber  than  does  the  0.26u 
wavelength  radiometer.  It  seems  evident,  then,  that  the  above  data  provides  another 
fairly  strong  indication  of  shock  front  irregularities. 

6.2.3  Radiometer  Signal  Rise  Times 

It  also  is  seen  in  Figure  6-7  a,  b that  the  rise  times  for  the  main  pulse  appear 
to  be  approximately  28  ysec.  for  both  ultraviolet  radiometers.  As  is  shoxm  in  Figure 
6-8,  the  rise  time  of  the  equivalent  electronic  systems  used  on  both  radiometers  is 
~ 12  psec.  A root  mean  square  correction  using  this  system  response  time  modifies 
the  observed  rise  times  of  the  main  pulses  to  ~25  psec.,  corresponding  to  complete 
filling  of  the  radiometer  field  of  view.  Some  attempts  were  made  to  correlate  this 
rise  time  with  shock  front  velocity,  but  too  many  uncertainties  are  Involved  to  allow 
for  meaningful  correlation.  These  uncertainties  Include  shock  front  non-planar  geo- 
metrical and  thermal  structure  and  the  unknown  varying  depths  of  optical  penetration 
as  the  shock  front  first  enters  and  eventually  fills  the  radiometer  fields  of  view. 


THE  HORIZONTAL  SCALE  IS  12,5  SEC/CM,  REAL  TIME 
(0.2  MSEC/CM  OSCILLOSCOPE  SWEEP).  THE  RISE  TIME 
IS  SEEN  TO  BE~12hSEC. 


Figure  6-8. 


Eler-ttonlc  System  Response  to  a Fast  (1  psec)  Calibration  Pulse 


7.  CONCLUSIONS 


This  study  has  developed  the  necessary  theoretical  background  for  showing  the 
feasibility  of  using  properly  designed  one  and  two  color  pyrometer  systems  to  measure 
the  temperature  history  of  explosively  generated  shock  wave  gases.  The  experimental 
data  from  the  DABS  II-D-2  field  test  were  analyzed  and  showed  excellent  correlation 
with  theoretical  results.  The  information  obtained  from  the  design,  calibration, 
installation  and  operation  of  this  system  provides  a sound  technical  basis  for  design- 
ing a second  generation  optical  pyrometer  system  that  is  compact  in  size,  easy  to 
install  and  simple  to  operate. 

The  main  objective  of  this  program  was  to  demonstrate  the  practical  application 
of  using  optical  radiometers  to  measure  airblast  temperature  histories.  This  goal 
was  clearly  demonstrated.  The  principle  results  can  be  summarized  as  follows: 

• Successful  operation  of  the  ultraviolet  and  infrared  channels  was  demonstrated 
in  the  field. 

• Air  shock  temperature  data  was  accurately  measured  by  ultraviolet  pyrometers 
and  correlated  with  theoretical  predictions 


• The  mechanical  design  of  the  open  port  installation  provided  protection 
against  damage  from  the  shock  wave. 


8.  RECOMMENDATIONS 


8.1  Identification  of  Problem  Areas 

The  Implementation  of  a new  system  usually  uncovers  certain  problem  areas  that 
cannot  be  fully  anticipated  prior  to  the  calibration,  fielding  and  testing  of  the  system. 
The  DABS  II-D-2  Temperature  Measurement  experiment  was  no  exception.  The  only  major 
problem  Involved  the  loss  of  all  three  infrared  channels.  The  other  problems  were 
minor  in  nature  and  dealt  with  Improved  design  changes  or  convenience  changes  for  im- 
proved field  operation. 

8.1.1  Electrical  Interference  on  Infrared  Channels 

The  Infrared  channels  had  been  checked  out  and  calibrated  under  full  field 
conditons.  Data  from  a pulsed  light  source  was  recorded  in  the  instrumentation  van 
and  observed  to  be  in  proper  working  order  (see  Figure  8-1).  The  dry  run,  the  day 
prior  to  the  test,  and  a quick  check  three  hours  before  detonation  showed  that  the 
system  was  in  perfect  operating  conditions.  Trouble  occurred  during  the  electrical 
surge  produced  by  the  fiduciary  pulse  as  sho^-n  in  Figure  8-2.  Close  proximity  of  these 
cables  to  our  experimental  station  and  land  line  hook-up  could  have  been  a factor. 

A preliminary  evaluation  indicates  that  the  cause  of  the  problem  could  have  been 
the  unshielded  electrical  leads  on  the  infrared  channels  within  the  optical  detector 
box.  A one  foot  section  of  cabling  from  the  voltage  preamps  to  the  preamp  gain  con- 
trols was  unshielded.  This  cabling  may  have  acted  as  an  antenna  for  the  3000  volt 
detonation  pulse.  It  is  felt  that  proper  shielding  and  grounding  of  all  electrical 
circuits  should  be  checked  and  corrected  where  necessary. 

8.1.2  Optical  Parallax 

The  two  ultraviolet  channels  were  located  at  either  extreme  in  the  five  lens 
cluster.  The  lateral  distance  between  channel  3 and  channel  5 was  12.7  cm  (5  Inches) 
at  the  lens  plane  and  approximately  15.2  cm  (6  inches)  at  the  detector  plane.  Thus 
allowing  each  detector  to  be  focused  on  a different  point  in  the  shock  front.  This 
phenomenon  was  apparent  from  the  different  arrival  times  as  discussed  in  Section  6.0. 
While  the  parallax  allowed  for  calculation  of  the  optical  thickness  of  the  0.36p  wave- 
length channel,  data  within  the  shock  wave  was  not  observed  at  a point  perpendicular 
to  the  optical  axis.  Repositioning  the  detectors  with  respect  to  the  optical  axis 
will  solve  any  parallax  problems. 


HORIZONTAL  TIME  TRACE  FOR  ALL  CHANNELS  IS  1.25  MSEC/CM, 
REAL  TIME  (20  MSEC/CM  OSCILLOSCOPE  SWEEP) 

TRACE  (a):  THE  1.9  ^ INFRARED  RADIOMETER,  CHANNEL  1 

TRACE  (b):  THE  4.5  ^ INFRARED  RADIOMETER,  CHANNEL  4 

TRACE  (c):  THE  2.9  ^ INFRARED  RADIOMETER,  CHANNEL  2 

TRACE  (d):  THE  0.26  p ULTRAVIOLET  RADIOMETER,  CHANNEL  5 

TRACE  (e):  THE  0.36  ^ ULTRAVIOLET  RADIOMETER,  CHANNEL  3 


Figure  8-1.  Pre-Test  Optical  Calibration  Curves  Recorded  at  the  DABS  II-D-2  Site 


HORIZONTAL  TIME  TRACE  IS  300  u SEC/CM, 
REAL  TIME  (0.5  SEC/CM  OSCILLOSCOPE  SWEEP) 
VERTICAL  SCALE  IS  0.1  VOLT/CM 


Figure  8-2.  Electrical  Noise  Recorded  on  Typical  Infrared  Channel  During  the 

DABS  H-D-2  Test 
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8.1.3  Optical  Alignment 


The  installation  and  set-up  of  the  baffle/pit  system  was  a time  consuming  job. 

The  design  and  construction  of  an  alignment  jig  would  greatly  alleviate  this  problem. 

Also  the  redesign  of  the  system  as  proposed  in  subsection  8.2  would  eliminate  this 
problem. 

8.1.4  Electronic  Packaging 

The  main  inconvenience  associated  with  the  electronics  was  the  need  to  open  the 
pit  and  remove  the  lid  from  the  optical/detector  box  in  order  to  activate  the  electronics 
power  supply.  Each  opening  increased  the  chance  of  damaging  the  electronics  and  op- 
tics from  dust,  moisture  or  debris.  This  problem  can  be  easily  eliminated  by  repack- 
aging the  electronics  so  that  all  switches  necessary  to  activate  the  system  are  lo- 
cated in  a remote  control  box  outside  the  main  pit.  Then  the  optical  box  and  pit  can 
be  sealed  once  the  equipment  is  Installed  and  aligned. 

8.2  Recommended  Hardware  Improvements 

The  problems  discussed  in  the  previous  subsection  can  be  corrected.  The  infrared 
channels  can  be  protected  by  shielding  all  electrical  cabling  between  preamplifier  cir- 
cuit and  controls.  Optical  parallax  can  be  corrected  by  rearranging  the  position  of 
the  optics,  the  optical  alignment  made  simpler  by  design  of  an  alignment  jig.  The 
inconvenient  arrangement  of  electrical  switches  can  be  alleviated  by  providing  for 
external  switching  of  electronics. 

But  the  main  recommendation  for  improvement  of  overall  performance  would  be  to 
reduce  the  package  size.  The  key  to  this  reduction  would  be  the  addition  of  a window 
(quartz  or  sapphire)  optically  transparent  in  the  required  wavelength  region  through 
which  the  shock  wave  would  be  observed.  The  addition  of  the  window  would  provide 
protection  from  pressure  without  obscuring  the  radiation  at  the  wavelength  of  interest. 
Without  pressure  to  worry  about,  the  detectors  can  be  brought  in  close  to  the  viewing 
port.  A lens  will  still  have  to  be  used  to  collimate  the  light  source.  Care  must  be 
taken  to  limit  the  transmission  loss  through  window  and  optics.  The  detector/optic 
package  would  be  a small  sealed  unit  already  prealigned  for  easy  installation.  All 
electronics  not  required  to  be  in  proximity  with  a detector  can  be  remotely  located  to 
improve  reliability  and  eliminate  interference  from  shock  propagated  through  structure. 

This  new  design  would  reduce  installation  and  calibration  time,  allow  placement 
of  more  units  in  a greater  variety  of  locations,  and  many  of  the  problem  areas  dls- 
^'uased  In  subsection  8.1  would  become  moot. 
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APPENDIX  A:  BLACK  BODY  RADIATION  THEORY  j 

Black  Body  radiation  exhibits  a unique  relationship  of  absolute  temperature  T to  1 

Black  Body  hemispheric  spectral  radiant  emittance  Wt>(X;  T),  and  to  Black  Body  spectral 

radiance  I°(X;  T).  That  is,  within  any  sufficiently  narrow  spectral  interval, 

+ 1/2  AAi,  there  is  a very  useful  one  to  one  correspondence  between  WO(A,AA;  T)  and  T, 

and  between  Io(A,AX;  T)  and  T.  Also,  optical  systems  using  photo-detector  cells  give  a 

voltage  signal  output,  V^,  directly  proportional  to  the  radiant  power  input  impinging  upon 

c 

the  objective  element  of  the  optical  system  and  directed  within  a field  of  view,  “ , that 
is  focused  upon  the  detector  cell.  Thus,  a spectrally  filtered  radiometer  optical  system 
responsive  only  to  a preselected  Xi  +1/2  AX^  spectral  Interval  and  viewing  Black  Body 
radiation  provides,  in  theory  a direct  measure  of  the  spectral  radiance  Io(X,AX;  T),  and 
hence  of  T. 

Planck's  original  form  of  the  Black  Body  spectral  radiation  law  was  specified  in 
terms  of  total  hemispheric  spectral  radiant  emittance  Wo(X;  T)  rather  than  the  directed 
spectral  radiance  I°(X;  T) . Planck's  Law  states  that: 

WO  (X;  T)  = 211  c2h  X-5[exp  (|^)  - 1]'^  (A.l) 

For  present  purposes,  however,  the  formulation  in  terms  of  directed  spectral  radiance, 

Io(X;  T)  is  more  applicable.  For  a "Lambertian  surface",  or  for  (as  in  the  present  appli- 
cation) "optically  thick"  layers  of  radiating  gases: 

I°(X;  T)  = 11-1  wO(X;  T)  = 2 c2h  X-5  [exp  - I]"!  (A. 2) 

In  both  equations  (A.l)  and  (A. 2)  c,h,  and  k as  usual  symbolize,  respectively,  the  speed 
of  light,  Planck's  constant  and  Boltzmann's  gas  constant.  The  spectral  wavelength,  X,  Is 
to  be  expressed  in  units  of  centimeters.  Equation  (A. 2)  often  is  exhibited  in  the  con- 
venient form: 

I°(X;  T)  = Cj^  X”5  [exp  (C2/XT)  - 1]'^  (A. 2a) 

with  the  first  and  second  radiation  constants,  Cj^  and  C2  defined  by  the  relations: 

c-^  = 2 c2h  = 1.1910  • 10”12watts  cm2  ster”!  (A.  3) 

C2  H = 1.4388  (cm°K)  (A. 4) 

The  definition  given  in  equation  (A. 4)  is  universally  used  for  C2,  but  the  literature 
often  contains  alternative  definitions  for  Cj^  depending  upon  geometrical  considerations 
and  also  depending  upon  whether  the  radiation  quantities  are  expressed  in  wavelength  units, 
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as  is  done  here,  or  in  wavenumber  units  as  is  often  more  convenient.  In  any  case,  the  c^ 
constant  shortly  will  be  eliminated  from  further  consideration. 

Typically,  Io(A;  T)  is  presented  in  units  of  watts 
[cm2  . ster  • (micron-wavelength) ]~1 , although  occasionally  units  of  watts 
[cm2  . ster  • (cm-wavelength) or  watts  [m2  . ster  • (micron-wavelength) ]~1 , 
etc.  are  used.  For  application  to  finite  spectral  bandwidth  measurements: 


(A. 2b) 


I°(A,AX;  T)«lO(X;  T)  • (^)  = c^  X [exp  (^  - 1]"^  • (|^) 

In  equation  (A. 2b)  it  was  assumed  that  I0(X;  T)  is  given  in  units  of  watts 
[cm2  . ster  . 


Figure  A-1  presents  log-log  plots  of  Io(X;  T)  against  wavelength  X and  also 
against  wavenumber  v = [X  (in  cm)]~^  for  selected  values  of  absolute  temperature  T.  It 
is  seen  that  in  log-log  plots,  the  maximum  values  of  I°(X;  T)  for  each  T lie  on  a straight 
line.  The  equation  of  this  parametric  maxima  curve,  known  as  Wien's  displacement  law, 
is : 

^max  = 0.28978  (cmOR)  = 2897.8  (p  °K)  (A. 5) 

where  r;^^  is  defined  below  in  equations  (A. 7),  (A. 8)  and  (A. 9). 

Also,  the  maximum  values  of  I°(X;  T)  for  any  T are  given  by  the  expression: 

I°(Xmax;  T)  = 4.098  • 10”^^  (T/l^K)^  watts  [cm2  . ster  • p]~^  (A. 6) 

Equations  (A. 5)  and  (A. 6)  are  obtained  by  setting  dl°  (X;  T) 

dX  “ ^ 

and  solving  for  first  and  then  I°(Xm3x>  T) . Before  differentiating,  it  is  convenient 

to  re-write  equation  (A. 2a)  in  the  form: 


1°  (x)  = cj^  (T/c2)^  [eX  - 1]-1, 


The  parameter  r^^  (which  equals  4.9651)  that  occurs  in  equation  (A. 5)  is  the  positive 
root  of  the  equation; 

X = 5 [1  - exp  (-x)]  (A. 9) 

It  is  evident  from  equation  (A. 7)  that  isothermal  spectral  radiance  ratios 

^2 

I®(X2^;  T)/I°(X2;  T)  are  functions  only  of  the  dimensionless  parameters  xi  = — - and 
Co  '^1^ 

X2  = . There  exist  several  convenient  condensed  tables  of  radiametrlc 

quantities  that  tabulate  only  I°(5f)/I°(^niax’  corresponding  WO(x)/W'^(X^^j^;  T)  etc., 

together  with  a short  table  of  I°(Xn,_x»  T)  etc. 
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A.l  Optimum  Spectral  Regions  for  Temperature  Determination 

In  radiation  temperature  measurements  it  would  be  preferable  to  choose  wavelengths 
X in  spectral  regions  where  the  monochromatic  spectral  radiance  changes  relatively 
rapidly  with  a change  in  temperature.  An  expression  for  the  wavelength,  X^^^  diff’ 
at  which  the  maximum  change  in  I°(X;  T)  occurs  for  a given  change  in  T is  found  by 
solving  the  triple  derivative  expression: 


3^  T) 

3^X3T 


(A. 10) 


The  resultant  condition  is: 


T 


^max  diff 


= 2404  (u  "K) 


(A. 11) 


Thus,  for  example  a near  infrared  wavelength  of  X = 2.4u  would  be  an  excellent  choice 
(other  factors  permitting)  for  measuring  radiation  temperatures  near  1000  “K,  while 

O 

ultraviolet  wavelengths  near  X « 0.24y  (=  2400  A)  would  be  preferable  for  temperatures 

near  10,000  °K.  A comparison  of  equations  (A. 5)  and  (A. 11)  shows  that,  for  a fixed 

falls  on  the  short  wavelength  side  of  X 
max  axil  max 

For  radiation  temperature  measurements  of  non-Black  Bodies  such  as  the  volume  of 

hot  radiating  gases  which  represent  the  current  problem,  however,  it  will  be  seen  that 

both  physical  and  practical  system  constraints  apply  a severe  limitation  upon  the 

wavelength  regions  available  for  use. 

A. 2 Wien's  Law  and  Raylelgh-Jeans  Law  Approximations  for  Black  Body  Radiation  Expressions 

There  are  two  widely  used  approximations  for  Black  Body  radiation  expressions  that 
find  application  to  different  (X,  T)  regimes.  Where  they  apply,  both  are  useful  for 
quick  calculations  and  for  analyses.  The  first  of  these  approximations  is  valid  when 
X ^ 2.5.  It  will  be  recalled  that  the  useful  dimensionless  parameter  x is  equal  to 
(C2/XT)  [or  (hc/kXT],  with  X measured  in  cm  - units  (1  cm  = lO^y)  for  C2  = 1.4388 
(cm  "K) . Under  these  conditions,  equation  (A. 2a)  reduces  to  the  form: 

I°(X;  T)  =«  c^x“^  exp  (-C2/XT)  ; for  x ^ 2.5  (A. 12) 

Equation  (A. 12)  or  its  equivalent  in  terms  of  hemispheric  radiant  emittance  W°(XT)  is 
known  as  "Wien's  Law".  The  other  approximation,  known  as  the  "Raylelgh-Jeans  Law"  is 
valid  when  x:i.0.15.  It  then  follows  that  equation  (A.  2a)  reduces  to  the  approximate 
expression: 
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I°(X;  T)  = (XT/C2)  = (c^/c2)  X”'^  T 

(A. 13) 

« Cj^(Tx/c^)^  • X ^ ; for  x <,  0.15 


A. 3 Error  Analysis  for  "One  Color"  Radiometric  Temperature  Determination 

It  can  readily  be  shown  that  the  fractional  error  in  temperature  determination  by 
"one  color"  radiometer  techniques  is  never  greater  than  the  corresponding  fractional 
error  in  the  measured  spectral  radiance.  For,  from  equation  (A. 7)  rewritten  in  the 
form: 


In  I°(x)  = In  [Cj^X  - In  (e^  -1) 


(A. 7a) 


and  recognizing  that. 


dx 

dT 


= - xT 


-1 


(A. 14) 


XT 


it  follows  that 


dT 

T 


dl!  r.exj.  (X)  - ll  ^ ^ . 

jO  X exp  (x)J  jO 


(A. 15) 


f 


il 

T 


[exp  (-x) ] 


dI°(X;  T) 
I°(A;  T). 


=»  [1  - x] 


dI°(X;  T) 
. I°(A;  T) 


; 0 < X ^ 0.15 

AT  ^ 9.6  • 10^  (y  . -K) 


(A. 15b) 


In  the  Intermediate  (A;  T)  regime  of  0.15  C x < 2.5,  the  factor  [E(x)  = [e^  - 1)/ 

(x  e ) ] still  remains  less  than  unity  and  monotonically  (^ecreases  with  increasing  x 
to  finally  approach  x ^ asymptotacally,  as  is  seen  from  Table  A-1. 


TABLE  A.1  PYROMETRIC  ERROR  FACTOR  E(x) 

for  one  color  pyrometer 


C2 
- AT 

E(x) 

We in's  Approx. 

E(x)  ~ X ^ 

Rayleigh- Jeans* 
Approx, 

E(x)  (l-x) 

0.01 

0.995 

0.99 

0.05 

0.975 

0.95 

0.10 

0.952 

0.90 

0.20 

0.906 

0.80  "1 

O.AO 

0.82A 

0.60 

0.60 

0.787 

0.50 

0.80 

0.688 

1.25 

0.20 

0.00 

0.632 

1.0 

1.50 

0.518 

0.667 

2.00 

0.432 

] 0.500 

A. 00 

0.245 

0.250 

8.00 

0.125 

0.125 

Table  A-1  and  the  preceding  discussion  illustrate  again,  as  was  pointed  out  in  the  earlier 
discussion  of  equation  (A. 11)  [concerning  the  wavelength  A^^^  diff  greatest 

change  in  I°(A;  T)  occurs  for  a given  change  in  T]  that  as  the  temperature  regime  of 
interest  Increases  it  would  be  preferable  (if  feasible)  to  utilize  progressively 
shorter  wavelengths  A so  as  to  increase  x and  thus  decrease  the  error  factor  E. 
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Figure  A-2  Is  a plot  of  the  "one  color"  pyrometric  error  factor  (E(x)  against  T 
for  several  selected  values  of  wavelength,  X,  ranging  from  O.lp  to  lOp • It  clearly 
shows,  as  does  Table  A-1,  that  the  measurement  error  (6I/I)  is  always  attenuated  rather  j 

than  amplified  in  determing  fractional  errors.  j 

< 

I 

A. 4 Two-Color  Pyrometer  Systems 

Equation  (A. 2b)  could  be  used  to  plot  a I°(X,  AX;  T)  vs.  T curve  for  a pre- 
selected fixed  X and  AX,  which  in  turn  could  be  used  in  theory  as  a measure  of  T as 
a function  of  received  radiation  signal  (i.e.,  output  voltage  from  the  X ± 1/2  AX 
sensitive  detector)  if  an  absolute  calibration  is  first  made.  However,  in  order  to 
avoid  the  difficulties  and  Inherent  errors  in  an  absolute  calibration,  and  also  the 
necessity  and  excessive  cost  of  perhaps  keeping  AX  extremely  narrow,  it  is  often  pref- 
erable to  go  to  a more  readily  implemented  "two  color"  (or  "multi-color")  system. 

These  systems  utilize  two  or  more  detector  systems,  each  with  similar  geometry  and  each 
filtered  to  be  responsive  to  different  X ± 1/2AX  spectral  intervals.  The  temperature 
indicators  then  are  the  dimensionless  spectral  radiance  ratios 

I°(X^,  AX^;  T)/I°(Xj,  AXj ; T) . For  reasons  and  compromises  that  will  be  explained  in 
Appendix  D,  six  different  wavelength  Intervals  X^  + 1/2  AX^  eventually  were  selected 
for  use  in  DABS  Il-D-1.  Five  of  these  same  wavelength  intervals  will  be  utilized  in 
DABS  II-D-2.  Figure  A-3  exhibits  I°(X^,  AX^;  T)  vs.  T curves  for  five  of  these  six 
selections;  the  X^  = 4.3y  and  X^  = 4.5vi  curves  (selected  for  use  in  observing  combustion 
products  and  4%  CO2  in  air,  respectively)  lie  too  close  together  to  exhibit  separately 
on  the  chosen  semi-log  scale.  The  vertical  lines  superimposed  on  Figure  A-1  also 
Indicate  the  selected  X^  values.  Figure  A-4  exhibits  the  four  "temperature  signal" 
ratios  of  [I°(X^,  AX^;  T)/I°(Xj,  4Xj ; T) ] selected  as  most  useful  for  determinations  of 
radiometric  temperature  history  in  air,  in  air  plus  4%  CO2  (presuming  this  addition 
would  be  allowed),  and  in  the  combustion  gases. 

A. 5 Applicability  of  "Two-Color  Pyrometer"  Techniques 

As  has  been  previously  stated,  a "two-color  pyrometer"  consists  of  a pair  of 
radiometers  for  which  filters  and  detector  cells  have  been  selected  to  confine  response 
to  pre-selected  spectral  increments  X^  + 1/2  AX^^  and  X^  + 1/2  AX^,  respectively. 

Assuming  that  thermal  equilibrium  is  a valid  assumption  for  the  radiating  surface  or 
volume  being  observed,  it  is  clear  that  the  "two-color  pyrometer"  can  provide  a useful 
estimate  of  the  radiation  temperature  under  any  of  the  following  sets  of  conditions: 
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I. 


Radiating  object  is  a Black  Body 


T)  - e(Xj,  AXj ; T)  - 1 for  all  A and  all  T 


II.  Radiating  object  is  a Gray  Body 

e(A,,  AA  : T)  = e(A.,  AA,;  T)  = e < 1 
1 i j J gray 


e = some  value  < 1,  but  essentially  constant  for  all  A,  AA,  and  T 
gray 

within  some  specified  bounds. 

III.  The  radiating  volume  is  not  "Black  and  Not  Gray,  and  has  a spectrally 
variable  emissivlty  e(A,AA;  T) , but  it  is  known  that  for  the  selected 
A . and  A , Intervals 

1 j 

c(A.,  AA;  T)  « e(A  , AA;  T)  < 1 
1 i J J 

Over  some  T-interval  of  Interest. 

IV.  The  radiating  volume  has  a spectrally  variable  emissivlty 

e(A,  AA;  T)  , and 
e(A^,AA^;  T)  ^ 

e(A.,AA.;  T)  < 1 and 
1 1 

e(Aj,AAj;  T)  < 1 


But  both  of  the  e-values  are  accurately  known  or  can  be  accurately 
computed  over  the  T-interval  of  interest. 


V.  The  Radiating  volume  again  has  spectrally  variable  emmlslvity, 
but  for  the  pre-selected  A^^-  and  Aj-  intervals,  the  respective 
T-values  are  both  ^ 3,  so  that. 


c(A^,  AA;  T)  « (A^,  AA^ ; T)  « 1 


As  was  pointed  out  in  the  Introduction,  Cases  1 and  II  are  not  applicable  to  gas  volumes 
while  Cases  III  and  IV  generally  are  very  difficult  to  implement.  Therefore,  the  only 
case  of  current  interest  is  Case  V.  It  will  be  shown  in  Appendix  E that  the  set  of  six 
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wavelength  Intervals  selected  in  Section  1.4.4  and  its  sub-sections  correspond  to 
suitably  large  t-values  >3.  A brief  error  analysis  for  "two  color"  pyrometers  is 
presented  in  tue  subsequent  sub-section. 


A. 6 Error  Analysis  for  "Two-Color"  Pyrometer  Systems 

It  was  stated  in  sub-section  A. 4 that  the  temperature  signal  for  a "two-color" 
pyrometer  is  the  ratio  of  measured  spectral  radiances: 


R^j(T)  = I(X^,  AX^;  T)/I(Xj,  AX^;  T) 
or  = Ijl. 


(A. 16) 
(A. 16a) 


where  the  second  form  is  a convenient  shorthand.  Then  since, 


In  R.  . = In  I - In  I, 

i j 


(A. 17) 


it  follows  that  the  fractional  error  in  determining  R^^  is  given  formally  by  the 
expression: 


(A. 18) 


Here  the  (61^/1^)  and  (51^/1^)  terms  represent  the  fractional  errors  in  measuring 
and  Ij , respectively.  It  is  clear  then,  that  if  the  errors  were  entirely  random, 
that  the  predicted  fractional  error  in  R^^  would  be  given  by  the  expression: 


6R 


1/2 


(A. 18a) 


(predicted  for  random  measurement  errors) 


If  the  major  errors  in  measuring  and  are  systematic  rather  than  random, 
then  different  types  of  error  expressions  would  have  to  be  used.  If  it  were  known 


same 

sign. 

then 

6R,  . 

61,  61, 

Ij 

_A 1 

R. 

I.  1. 

Ij 

i j 

(A. 18b) 


(for  systematic  errors  with  (51^) 


(«I^)  > 0) 
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If,  on  the  other  hand,  it  were  known  that  the  systematic  errors  had  opposite  signs 
then 


“l.1 

+ 

61^ 

"ij 

'il 

^3 

(A. 18c) 


(for  systematic  errors  with  (51^) (61^)  < 0) 


Now  each  of  the  (61/ I)  - quantities  are  related  to  (6T/T)  by  the  relation  given  in 
Equation  (A. 15),  e.g.. 


(61. 


-1  r±  ® 

^ (6T/T) 


((5T/T) 


(A. 15a) 


e ^ -1- 


where  = 1,4388  x 10^  (p  . °K)/(X^  - T) 


(A. 8a) 


when  is  specified  in  microns  rather  than  in  centimeters. 


It  follows,  by  combining  equations  of  the  forms  (A. 18)  and  (A. 15),  that  if 
the  61  measurements  errors  each  could  be  resolved  into  three  increments  consisting 
respectively  of  random  errors  and  systematic  errors  of  the  two  types  characterized 
in  equations  A. 18b)  and  (A. 18c),  then: 


T 


-2  -2 
E + E 

-1/2 

1 

•H 

1 

^li 

(A. 19a) 


for  random  measurement  errors  in  I^  and  Ij  and  with 


(6R. 


given  by  Equation  (A. 18a) 
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for  systematic  measurement  errors  in  and  and  with  (61^) (61^)  > 0 
and  with  [ '^^i.1  \ given  by  equation  (A.  18b); 

V 


(A. 19b) 


6T 

-1 

T 

Ei  +Ej 

(A. 19c) 


for  systematic  measurement  errors  in  and  1^  and  with  (61^) (61^)  < 0 


and  with 


given  by  Equation  (A. 18c) 


Figure  A-5,  A-6,  and  A-7  present  plots  of  the  three  "two-color"  pyrometer 
error  factors  given  respectively,  by  equations  A. 19a,  A. 19b,  and  A. 19c  for  selected 
Aj)  pairs.  It  is  seen  that  the  two  error  factors  given  by  equations  A. 19a  and 
A. 19c  are  always  error  attenuators;  i.e.,  the  predicted  (6T/T)  increments  are  always 
less  than  the  corresponding  random  and  systematic  fractional  error  increments  of  (6R/R) . 
However,  the  error  factor  given  by  Figure  A-6  for  the  equation  A. 19b  relation  for 


systematic  error  Increments,  such  that  (61.) 

’ 1 systematic 


(<SI^)  ^ > 0»  sometimes 

j systematic 


becomes  an  "error  amplifier"  rather  than  an  "error  attenuator"  (i.e.,  the  error  ;Eactor 
becomes  greater  than  unity).  This  is  seen  to  happen  for  IROO'K,!^  = 2.8p  and 
AjWSp.  For  the  (X^,  X^)  pair  of  (2.8p,  4.5u),  this  "error  amplification"  would  occur 
for  T ^ 1500°K.  This  "error  amplification"  is  not  particularly  serious,  however,  since 
the  error  quantity  being  "amplified"  is  really  of  the  type. 


6^ 

61, 

e,  61, 

i 

i 

J j 

1 ^1 

9A 


6.  Two  Color  Pyrometer  Error  Factor  (Equation 
Systematic  Measurement  Errors  of  the  Same  Sign 


Figure  A-7.  Two  Color  Pyrometer  Error  Attenuation  Factor  (Equation  A. 19c) 
For  Systematic  Measurement  Errors  of  Opposite  Sign 


I 


where  the  0^^  and  6j  represent  those  fractions  of  61^  and  61^  which  are  both  systematic 
rather  than  random  and  also  have  the  same  sign.  It  is  clear  that  this  preceding  type 
of  measurement  error  Increment  could  be  vanishingly  small,  and  in  general  is  probably 
smaller  than  the  measurement  error  increments  multiplied  by  the  curves  shown  in  either 
Figure  A-5  or  A-7.  Also,  pairs  of  near  Infrared  radiometer  channels  are  of  use  primarily 
for  relatively  low  temperature  measurement  rather  than  at  the  higher  temperatures  where 
the  amplification  factor  becomes  quite  large. 


F 


f 

I 
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APPENDIX  B:  A BRIEF  INTRODUCTION  TO  NOMENCLATURE  OF  MOLECULAR  BAND  SYSTEMS  [3]  < 

The  designations  of  molecular  electronic-vibration-rotation  band  systems  are  his- 
torical, being  either  named  after  one  or  more  of  the  original  investigators,  or  using 
names  applied  by  an  original  investigator.  Electronic  energy  levels  are  designated,  more 
systematically,  with  the  electronic  ground  state  (or  for  a few  molecules,  the  state  that 
Initially  was  mistakenly  taken  to  be  the  electronic  ground  state)  being  designated  by  the 
symbol  X.  Then,  the  symbols  A,B,C,  etc.  designate  successively  higher  lying  electronic 
states  all  of  which  have  the  same  "multiplicity"  as  the  X-state  and  thus  are  capable  of 
"permitted"  electric  dipole  transitions  with  the  X-state  or  with  each  other.  "Forbidden" 
transitions  can  occur  but  are  very  much  weaker  than  "permitted"  transitions.  Electronic 
states  with  "multiplicity"  differing  from  that  of  the  ground  state  generally  are  desig- 
nated by  a lower  case  letter,  a,b,c,  etc.,  again  in  order  of  increasing  electronic  exci- 
tation. 

Each  electronic  state  of  a molecule  has  associated  with  it  a sequence  of  quantized 
vibrational  energy  sub-levels,  and  each  of  these,  in  turn  has  associated  with  it  a sequence 
of  quantized  rotational  energy  levels.  In  general,  the  energy  level  separation  is  suc- 
cessively much  greater  between  electronic  states  than  between  vibrational  sub-levels  of 
these  states,  and  between  these  latter  sub-levels  than  between  the  rotational  sub-sub-levels. 

For  a given  electronic  transition  the  corresponding  upper  and  lower  energy  levels 
are  tagged  by  their  corresponding  electronic,  vibrational  and  rotational  quantum  numbers 
designated  by  the  symbols,  n,v,  and  J respectively.  Sometimes  additional  rotational  quan- 
tum numbers,  K,L,  etc.  also  are  required.  For  polyatomic  molecules,  v must  be  replaced  by 
a set  of  vibrational  quantum  numbers,  Vj^,  V2,  v^,  etc.  each  designating  the  degree  of  vi- 
brational excitation  of  the  Indexed  normal  vibrational  modes  1,2,3,  etc.  For  convenience, 
single  and  double  "primes"  are  used  to  designate  the  respective  upper  and  lower  states  in 
a transition.  Thus  (neglecting  K,L.  etc.  and  "fine"  structure)  the  most  general  transit- 
ion corresponding  to  a single  spectral  line  of  a diatomic  molecule  would  be  of  type 
(n",  v",  J")  ->•  (n* , v' , J'),  where  the  lower  arrow  represents  absorption  and  the  upper  e- 
mlsslon.  For  transitions  between  two  specified  electronic  states,  say  A and  X,  these 
state-symbols  are  substituted  for  the  n'  and  n"  symbols,  i.e.,  (X,  v",  J")  ^ (A,v',J'). 

The  totality  of  all  such  "lines"  could  be  known  as  the  A-X  band  system,  etc. for  the  speci- 
fied molecule.  Instead,  however,  the  band  systems  generally  are  designated  by  some  his- 
torical name.  Individual  vibrational  8ub-"bands"  within  an  electronic  system  correspond 
to  all  the  allowed  J"  ^ J'  transitions  for  some  specified  pair  of  (v",  v')  levels.  These 
J-lines  are  restricted  by  "selection  rules"  to  either  only  AJ  - J'  - J"  = + 1 or  only 
+ 1,  0.  The  nominal  "center"  of  each  band  system  is  designated  as  (in  cm“l  units) 
and  corresponds  to  the  sometimes  hypothetical  transition 
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(n",  v"  “ 0,  J"  “ 0)  ->■  (n*  v'  • 0,J'  “ 0).  The  relative  Intensities  of  the  various  vi- 
brational sub-bands,  and  of  the  Individual  rotational  "lines"  within  a band  depend  upon 
quantum  mechanical  considerations  and  also  upon  the  relative  population  levels  of  each 
energy  sub-state.  The  population  levels,  in  turn,  are  very  strongly  temperature  de- 
pendent through  the  Boltzmann  distribution  law. 


Most  Infrared  spectra  consist  of  transitions  be  veen  different  vibrational-rota- 
tional sub-levels  of  the  same  electronic  state,  i.e.,  n"  h n' . Except  at  very  high  tem- 
peratures , , the  overwhelmingly  predominant  contributions  are  usually  from  the  ground 
electronic  state.  Again,  the  totality  of  "allowed"  transitions  between  all  the  rotation- 
al sub-levels  of  a pair  of  vibrational  states  are  said  to  constitute  a single  vibration- 
al band.  For  diatomic  molecules  this  would  be  called  a (v"  ■+  v')-band.  At  low  and  mod- 
erate temperatures,  only  transitions  with  low  v"  (I.e.  0,  1,  2 ..)  are  Important.  But 
as  the  temperature  Increases,  higher  vibrational  levels  become  populated,  and  "hot  band" 
transitions  originating  from  higher  v"  lower  states  begin  to  make  Increasingly  Important 
contributions.  The  spectral  regions  for  the  same  Av  = v'  - v"  tend  to  partially  over- 
lap, but  as  the  temperature  Increases,  the  overall  resultant  band  structure  generally 
broadens  and  the  region  of  maximum  absorption  or  emission  can  shift  very  appreciably. 


B.l  Infrared  Vibrational  Bands  for  Polyatomic  Molecules 


A polyatomic  molecule  consisting  of  n atoms  generally  has  three  non-quantized 
transitional  modes  (I.e.  In  arbitrary  x,  y and  z directions),  three  normal  rotational 
modes,  and  3n-6  normal  vibrational  modes.  "Linear  molecules"  such  as  CO2  however,  have 
one  less  normal  rotational  mode  and  one  additional  normal  vibrational  mode,  since  their 
moment  of  Inertia  about  their  figure  axis  Is  essentially  zero. 


The  water  vapor  molecule,  H2O,  for  example,  has  the  symmetry  of  an  Isosceles 
triangle  with  the  0-atom  at  the  vertex.  Its  three  normal  modes  consist  of  symmetric  and 
anti-symmetrlc  "bond-stretching"  vibrations  labeled  and  V3,  and  a bond-bending  vibra- 
tion \>2-  These  three  fundamental  vibrations  are  symbolically  represented  below; 


( 1,0,0) 


(0,1,0) 


(0,0,1  ) 


I 


r 


T 


Any  arbitrary  vibrational  mode  is  designated  by  the  vibrational  index  triplet  (vj^,V2,V3) 
with  each  v Integer  quantum  number  denoting  the  degree  of  excltational  contribution  by 
the  corresponding  fundamental  mode.  As  was  indicated  below  the  preceding  figures, 

(1,  0,  0),  (0,  1,  0)  and  (0,  0,  1)  denote  excitation  of  the  first,  second  and  third  fun- 
damental vibrational  modes.  Similarly,  (v,  0,  0),  (0,  v,  0)  or  (0,  0,  v),  where  v is 
any  non-zero  Integer,  would  indicate  excitation  of  the  v^h  harmonic  [or  (v-l)bh  overtone] 
of  the  corresponding  fundamental  mode.  If  two  or  more  of  the  v-lntegers  are  non-zero, 
i.e.  (vj^,  V2,  0),  (vi,0,  V3) , (0,  V2,  V3) , or  (vj^,  V2,  V3),  then  the  indicated  excita- 
tion is  said  to  be  a "combination  mode".  Since  diatomic  molecules  have  only  a single 
fundamental  mode,  they  can  have  fundamental  and  harmonic  mode  excitation,  but  no  combi- 
nation modes. 

The  linear  CO2  molecule  (more  accurately  represented  as  0 « C “ 0 really  has  four 
fundamental  vibrational  modes,  consisting  also  of  symmetric  and  anti-symmetric  bond- 
stretching modes  and  V3  and  two  equivalent  (doubly  degenerate)  bond-bending  modes 

both  represented  by  the  symbol  V2  and  Indexed  by  V2,  where  the  i-superscript  represents 
a pseudo  rotational  quantum  number  resulting  from  a combination  of  degenerate  vlbra- 
tlonal  modes.  The  most  general  mode  designation  is  (v^^,  V2,  v^),  with  fundamental,  har- 
monic, and  combination  modes  defined  as  before.  The  fundamental  vibrational  modes  are 
symbolized  below.  The  designations  for  harmonic  and  combination  modes  are  similar  to 
those  for  water  vapor. 


o- 

o 


-o- 
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(1,0°  0) 


.♦-o ^1^3 

t ( 0,  0°  1 ) 

? 4 

/ y 

^ 

(0,1°  0) 


The  very  strong  CO2  absorption  centered  near  4.3)j  wavelength  arises  from  a "funda- 

mental  band".  At  low  temperatures  the  absorption  (or  emission)  is  mostly  from  rotational 

lines  of  the  fundamental  U3,  or  (0,  0°,  0)  -*•  (0,  0°,  1)  transition.  As  the  temperature 

i"  *■  t' 

Increases,  however,  transition  of  the  type  (v^,  V2,  v^)  (v^,  V2  , V3  + 1)  with 
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(vj^  ^2  ^ ^^3^  ~ ^ make  Increasingly  Important  contributions.  The  much  weaker  CO2  ah- 
sorbtion  centered  near  2.7p  wavelength  is  a combination  band.  At  low  temperatures  it  a- 
rlses  primarily  from  the  "accidentally  degenerate"  (and  "Fermi-Resonant")  pair  of  tran- 
sitions (0,  0°,  0)  X (0,  2°,  1)  and  (0,  0°,  0)  ^ (1,  0°,  1).  As  the  temperature  increa- 
ses, there  are  increasing  contributions  from  "hot  band"  transitions  of  the  type, 

(vj^,  V2,  V3)  ^ (vj^,  V2,  + 2,  V3  + 1)  and  (vj^,  V2,  V3)  ^ (v^  + 1,  V2,  V3  + 1),  [ the  SL" 

and  I'  superscripts  are  omitted  for  convenience]  with  (v^^  + V2  + V3)  > 1. 

The  relatively  very  weak  H2O  absorption  band  centered  near  1.87  is  also  a combi- 
nation band.  At  low  temperatures  it  arises  primarily  from  the  (0,  0,  0)  ■*  (0,  1,1) 
vibrational  transition.  With  increasing  temperature,  there  are  "hot  band"  contributions 

of  the  type  (Vj^,  V2,  V3)  ^ (vj^,  V2  + 1,  v^  + 1),  again  with  (vj^  ^^2  ^ 
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APPENDIX  C:  A BRIEF  INTRODUCTION  TO  SOME  RELEVANT  SPECTRAL  PARAMETERS 

It  was  mentioned  in  Section  1.4.1  that  the  absorption  coefficients  k^j(X;  T) 

are  arbitrarily  defined  parameters.  It  will  be  recalled  that  the  indexing  subscripts 

designate  the  band  system  [see  Appendix  B]  of  the  i^^  molecular  (or  ionic)  species 

which  may  exhibit  spectral  activity  in  the  X + | AA  spectral  region  of  interest.  [For 

an  atomic  species,  the  j subscript  would  index  an  individual  atomic  "line".]  The  truly 

fundamental  spectroscopic  parameters  are  the  "absorption  cross  sections", 

®ljm^^*  Xjij) , for  the  individual  m^h  Une  contours,  each  centered  on  some  X = Xj^  within 

the  band  system.  In  the  wavenumber  domain,  which  is  more  convenient  for  analyses,  the 

cross  sections  would  be  represented  as  with  the  line  center  wavenumber  po- 

-1  J m 

sltlon  if  X is  given  in  centimeters.  In  either  representation,  the  a's  are 

given  in  units  of  cm^  per  molecule  of  the  i^^  species  that  is  in  a proper  electronic 
lower  state  to  give  rise  to  absorption  in  the  band  system.  [see  Appendix  B.]  The 
k to  a relationship  will  be  discussed  in  Sub-Section  C.l. 

Other  very  important  parameters,  which  together  with  the  absorption  cross  sections 
determine  the  spectrally  smoothed  mean  absorption  coefficients,  R^j(X,AX;  T),  are  the 
spectral  widths  and  shapes  of  the  individual  "line"  contours  and  the  separations  of  ad- 
jacent spectral  lines.  These  considerations  will  be  discussed  in  Sub-Section  C.2. 

C.l  The  Relation  Between  Absorption  Cross  Sections  and  Absorption  Coefficients 


The  absorption  coefficients,  , are  related  to  the  fundamental  cross  sections. 


®ljm’  defining  equation: 


klj(X;  - ”o  ■ (T)  J 

Here  N"yju  represents  the  population  number  density,  in  units  of  cm“^,  for  the  fractions 
of  the  i^^  species  molecules  whose  state  of  internal  excitation  corresponds  to  that  lower 
energy  state  from  which  absorption  transitions  can  produce  the  spectral  line  (cen- 
tered on  X^)  of  the  band  system.  As  in  Section  1.4.1,  N^  represents  the  2. 67 "10^^ 
molecules  cm"^  corresponding  to  an  ideal  gas  at  STP  conditions,  while  Nj^  represents  the 
total  number  density  of  i^^  species  molecules  in  cm~^.  An  explicit  representation  of 


N"  (T) 
ijm 


,(X,  X ) 


the  pressure  dependence  of  N^  and  has  been  omitted  in  equation  (C.l)  because  the 

(N"ijm/N^)  ratio  is  not  pressure  dependent.  The  summation  in  equation  (C.l)  extends 


The  spectrally  "smoothed"  parameters  (X,  AX;  T)  result  from  averaging  the 

quantities  expressed  by  equation  (C.l)  over  some  suitable  X + | AX  spectral  interval.  It 

is  seen  that  a substitution  of  these  E^j-values  into  equation  (1.6)  would  cancel  out  the 

(N,/N  )-term  occuring  in  the  latter  equation,  leaving  the  more  fundamental  lower  state 
1 o 

population  density  terms 


For  the  assumption  of  thermodynamic  equilibrium,  the  (N"£jnj/N^)  ratios  in  equation 
(C.l)  exhibit  the  familiar  Boltzmann  exponentially  dependent  temperature  distribution. 

This  Boltzmann  dependence  causes  the  k(X;  T)  [and  consequently  the  corresponding 
k(X,  AX;  T)]  values  given  by  equation  (C.l)  to  exhibit  very  strong  temperature  variation 
when  the  lower  state  energy  levels  for  the  spectral  "Band  system"  belong  to  a highly 
excited  electronic  state  [see  Appendix  B].  Let  represent  the  generally  dominant 

electronic  portion  of  the  total  energy  of  the  lower  state.  Then  the  (N"j^jnj/N£)  are  very 
small  when  kT  <<  E^j  (the  k in  the  kT  terms  represents  the  Boltzmann  gas  constant),  and 
Increase  exponentially  with  increasing  T as  kT  first  approaches  and  then  exceeds  ; 
it  is  seen  from  equation  (C.l)  that  the  k^j(X)  absorption  coefficients  also  increase 
proportionately. 


The  increase  of  the  k^^  (X;  T)  parameters  with  T is  much  less  marked  when  the 
lower  states  for  the  absorption  "band  system"  belong  to  the  ground  electronic  state  (at- 
is  virtually  always  true  for  Infrared  vibration-rotation  bands  where  both  upper  and  lower 
electronic  states  are  the  same  electronic  ground  state).  Even  here  however,  increases  in 
temperature  increase  the  populations  of  more  highly  excited  vibration  (and  rotation)  levels 
and,  as  was  discussed  in  Appendix  B,  tends  to  Introduce  absorption  in  the  so-called  "hot 
bands".  That  is,  the  summation  in  equation  (C.l)  will  include  significant  contributions 
from  a larger  number  of  m-indexed  "lines". 


C.2  Spectral  Line  Contour  and  Line  Density  Considerations 

The  absorption  coefficients  K (X,  AX;  T)  used  in  the  present  computations  are  all 
spectrally  smoothed  mean  values  rather  than  spectrally  fine  structured  monochromatic 
values.  The  smoothing,  however,  is  generally  over  25  to  'v  50  cm~^  which  are  much 
narrower  sub-intervals  than  the  0.01  and  O.ly  AX-values  assumed  in  this  report. 

Also,  wherever  these  smoothed  values  exhibit  sharp  fluctuations  within  a AX-interval  of 
interest  either  a minimum  or  at  least  a low  average  was  taken.  Even  so  one  need  be 
concerned  with  the  possibility  of  some  of  the  actual  k(X;  T)  values  within  a AX-radio- 
meter  filter  transmission  region  falling  appreciably  below  the  lc(X,  AX;  T)-values  used. 
This  could  occur,  generally  only  at  low  temperature  and  low  pressures,  when  the  "line 
packing  factor"  (y/d)  drops  below  unity.  Here  y represents  the  individual  spectral  "line" 
contour  half  widths  (at  half-amplitude),  while  d is  the  distance  between  the  centers  of 
adjacent  spectral  "lines". 
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In  the  (P,T)-reglme  of  current  Interest,  the  spectral  lines  profiles  are  com- 
pletely dominated  by  collision  rather  than  Dopper  effects.  The  relative  line  profiles, 
as  functions  of  wavenumber  frequency  v in  the  vicinity  of  v = are  given  by  the  collis- 
ion profile  relation: 

A(v,  v^)  » y[v-v^j)^  + Y^]“^  (C.2) 

The  A(v,  Vjjj)  parameter  introduced  here  is  directly  proportional  to  the 

[N"^j^(T)/Nj^(T)  ] '^ni^  that  would  occur  in  equation  (C.l)  if  it  were  rewritten  in 

terms  of  wavenumber  rather  than  wavelength.  It  is  seen  that  the  maximum  A(v,  v ) value 

-1 

occurs  when  (v-v^)  = 0,  and  is  equal  to  Y . A(v,  Vu|)-values  of  (1/2y)  , (1/5y)  , (l/10*f) , 

(1/17y)...  occur  respectively  when  v-Vjjj  = + Y,  + 2y,  + 3y,  + 4y The  pressure 

broadened  "half  widths"  y have  relatively  complicated  dependence  upon  temperature  and 
species  partial  pressures  of  each  of  the  constituent  gases  [4].  Quite  often,  however,  a 
sufficiently  good  approximation  is  given  by  a relation  of  the  form: 


YaY  (P/P„)  • (T„/T)°-^  (C.3) 

o o o 

More  accurate  expressions  (see  Reference  [4],  pages  218,  219  and  223)  use  the  partial 
pressures  of  the  various  constituent  gases,  with  various  y^j  rather  than  Yq  and  additional 
powers  of  (T^/T)-terms  which  differ  from  0.5.  It  was  found  that  the  more  complicated  re- 
lations were  not  needed  in  the  current  problem  except  for  application  to  H2O.  The  very 
strong  dipole  moment  of  this  molecule  required  the  use  of  a (T^/T)^  term  as  well  as  sever- 
al terms  multiplied  by  mole-fractions  of  CO2,  N2  etc. 

Figure  C-1  shows  mean  values  of  the  "line  nacklne"  factor,  d ^ for  the  4.3  and 
2.7p  centered  CO2  bands  at  600  and  3000  °K,  the  1.87u  centered  H2O  band  at  the  same 
temperatures,  and  a portion  of  the  CO  - band  centered  at  4.9y  (but  extending  past  4.5p). 
Figure  C-2  gives  corresponding  y-values  [computed  from  the  more  correct  Reference  3 
equation  rather  than  the  equation  (c.3)  approximation]  for  CO2,  CO,  and  H2O  both  as 
they  would  occur  in  4%  CO2  + air  and  in  the  DABS  II-D  combustion  products. 

Figure  C-3  gives  schematic  representations  of  the  relative  line  profiles 
A(v,Vjj|),  and  of  "line  packing"  for  hypothetical  sets  of  constant  intensity  equally  spaced 
"lines"  for  which  the  (y/d)  factor  is  equal  to  1,  0.5,  and  0.25  respectively.  The  latter 
three  sub-figures  also  show  k(v)-curves,  with  k(v)  (proportional  to  the  absorption  co- 
efficient) given  by  the  relation: 
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Figure  C-3.  Spectral  Line  Profiles  and  Resultant  Emission  (or  AbsnTption)  Patterns 


00 


'^(v)  = A(v,  Vj||) 

in  B ^00 

A mean  value  horizontal  k (v)  line  also  Is  drawn  in  for  the  (y/d)  =0.5  and  0.25  repre- 
sentations, but  is  too  close  to  k(v)  to  show  for  the  (y/d)  » 1 representation.  The 
fluctuation  of  k(v)  is  equal  to  ± 8.5%  for  (y/d)  = 0.5,  and  to  "v  ± 60%  when  (y/d)  » 0.25. 

It  should  be  realized  that  the  y-values  are  directly  proportional  to  (P/P  ) 

o 

if  species  mole-fractions  remain  relatively  constant.  Thus  an  isothermal  two-fold 
increase  in  pressure  would  convert  a (y/d)  = 0.25  spectral  region  into  a (y/d)  « 0.5 
spectral  region,  and  the  latter  into  a (y/d)  = 1 spectral  region.  Similarly,  an 
appropriate  Increase  in  temperature,  while  it  would  tend  to  decrease  y,  could  excite 
enough  additional  "hot  lines"  to  produce  the  same  effect.  Thus,  as  we  stated  previously, 
the  concern  about  low  (y/d)-values  and  a consequent  drop  of  portion  of  k(X;  P,  T) 
below  the  corresponding  k(X,  ^X;  P,  T)  is  generally  relevant  only  for  low  T and  low  P 
values.  For  this  reason  it  was  not  considered  worthwhile  to  expend  any  of  the  limited 
time  and  funds  available  on  a literature  search  for  high  temperature  0^  and  NO,  etc., 
y-  and  d-values. 

When  (y/d)  drops  below  unity,  an  upward  adjustment  of  minimum  t-values  is 
required,  since  the  actual  t(X,  AX;  T,  P)-values  fluctuate  from  their  mean  value  t in 
direct  proportion  to  the  fluctuation  of  the  k-values  about  their  k.  Also,  in  real 
spectra,  any  computed  (y/d)  value  Implies  the  possibility  of  still  lower  spectrally 
fine  grained  (y/d)  values  within  the  6X(or  6v)  smoothing  Interval.  Figure  C-3b,  c and 
d also  give  the  listings  of  the  increased  x-values  (over  t ^ 3)  which  are  estimated 
to  be  required  for  the  corresponding  (y/d),  and  (y/d)  values  for  the  figures.  It 
is  estimated  that  (y/d)  “ 0.5  and  (y/d)  = 0.5  are  acceptable  for  x-values  3.3  and 
Z 5 respectively;  while  (y/d)  • 0.25  and  (y/d)  = 0.25  are  acceptable  for  respective 
T values  ^ 6 and  Z 10. 

After  investigation,  it  turned  out  that  the  Increased  minimum  t criteria  for 
low  (y/d)  was  pertinent  only  to  the  1.9m  water  vapor  spectral  region  of  the  X^-radlo- 
meter.  The  reason  for  this  was  the  proverbial  "ounce  of  prevention";  i.e.,  care  was 
taken,  whenever  possible,  to  select  the  radiometer  AX-increments  where  d ^ line  density 
did  not  fall  too  low.  The  Increased  t confidence  Interval  for  the  1.9m  band  is  hatched 
in  Figure  E-lOb  of  Appendix  E. 
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APPENDIX  D:  CONCENTRATIONS  AND  SPECTRAL  DATA  FOR  RADIATING  SPECIES  MOLECULES 

Table  D-1,  extracted  from  Reference  [5],  shows  the  accepted  normal  composition 
of  dry  atmospheric  air  near  sea  level.  Water  vapor  concentration  is  not  shown,  since 
it  is  very  highly  variable  both  spatially  and  temporarily.  An  appropriate  mean  value 
of  H2O  concentration  would  be  of  the  order  of  0.002  or  0.003%,  but  the  relative  vari- 
ations are  so  large  as  to  make  this  mean  value  of  very  little  use.  It  is  seen  that  the 
major  normal  constituents  of  atmospheric  air  are  N2  and  O2,  followed  by  Argon  and  CO2. 

As  air  temperature  is  increased  however,  first  O2  (and  CO2),  and  then  N2  tend  to 
dissociate,  forming  NO,  0,  N and  other  species.  Figure  D-1,  extracted  from  Reference  [6], 
gives  mole-fractions  of  major  dry  air  constituents  for  temperatures  of  from  1000°K  to 

8000°K,  and  for  a pressure  of  (P/P  ) equal  to  10  (standard  sea  level)  atmospheres. 

o 

As  will  be  seen  from  data  to  be  presented  in  Appendix  E,  the  mole-fraction  distributions 


Table  D-1  Normal  Composition  of  Clean,  Dry  Atmospheric  Air  Near  Sea  Level 


Constituent 

Gas 

Content 

Molecular 

Gas 

Symbol 

(%  by  volume) 

Weight* 

Nitrogen 

N 

78.084 

28.0134 

Oxygen 

0? 

Af 

20.9476 

31.9988 

Argon 

0.934 

39.948 

fCarbon  dioxide 

CO 

0.0314 

44.00995 

Neon 

Ne 

0.001818 

20.183 

Helium 

He 

0.000524 

4.0026 

Krypton 

Kr 

0.000114 

83.80 

Xenon 

Xe 

0.0000087 

131.30 

Hydrogen 

ci 

0.00005 

2.01594 

fMethane 

0.0002 

16.04303 

Nitrous  oxide 

N,0 

0.00005 

44.0128 

f Ozone 

Summer:  0 to  0.000007 

47.9982 

Winter:  0 to  0.000002 

47.9982 

fSulfur  dioxide 

SO 

0 to  0.0001 

64.0628 

fllitrogen  dioxide 

NO2 

0 to  0.000002 

46.0055 

f Ammonia 

nh:: 

0 to  trace 

17.03061 

t Carbon  monoxide 

co'^ 

0 to  trace 

28.01055 

flodine 

^2 

0 to  0.000001 

253.8088 

12  12 

*0n  basis  of  C Isotope  scale  for  which  C equals  12.0000. 
fThe  content  of  these  gases  may  undergo  significant  variations  from 
tine  to  time  or  from  place  to  place  relative  to  the  normal. 
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are  strongly  dependent  upon  pressure  as  well  as  upon  temperature.  The  equilibrium 
concentrations  of  O2,  CO2  and  N2  (and  later  CO  and  NO)  decay  appreciably  less  rapidly 
as  the  pressure  Increases.  It  Is  evident  from  Figure  D-1,  however,  that  NO  becomes  a 
major  constituent  of  high  temperature  air  and  that  any  CO2  present  tends  to  dissociate 
to  CO  (and  0)  for  temperatures  much  above  ~ 3000°K. 

At  normal  temperatures  and  pressures,  atmospheric  air  is  relatively  very 

transparent  in  the  visible  and  near  ultraviolet  spectral  regions.  This  fortunate 

lack  of  opacity  occurs  because  at  normal  atmospheric  temperatures  all  the  really  strong 

N-  and  0,  absorption  bands  occur  in  the  "vacuum  ultraviolet"  wavelengths  shorter  than 
^ o 

~ 2100  A.  Below  about  2000  A,  a few  centimeters  of  air,  at  one  atmospheric  pressure, 

are  completely  opaque  owing  to  the  N2  and  O2  absorption.  The  earth's  high  altitude 

ozone  layer  serves  as  protection  against  near  ultraviolet  radiation  in  the  ~ 2100  to 
o 

~ 3000  A region,  transmitting  only  the  less  hazardous  "sunburn"  ultraviolet  wavelengths 
above  ~ 3OO0X.  However,  it  can  be  s&en  from  Table  D-1  that  near  sea  level  (and  up  to  ~ 20 
kilometers  altitude)  the  normal  (but  highly  variable)  "parts  per  million"  ozone  concen- 
tration is  relatively  low.  In  the  relatively  short  path  lengths  of  current  interest, 
then,  the  ozone  absorption  spectral  region  is  available  for  use.  Except  for  some 
relatively  very  weak  "forbidden  transition"  O2  absorption  near  0.78ii,  N2  and  O2  also 
are  transparent  in  the  infrared.  Virtually  all  of  the  atmospheric  infrared  absorption 
and  re-emission  of  solar  radiation,  so  important  to  the  earth's  "heat  budget",  is  due 
to  the  minor  atmospheric  constituents,  and  primarily  to  H2O  and  CO2.  Figures  D-2  and 
D-3,  also  extracted  from  Reference  [5],  present  solar  spectra  (spectra  of  the  sun  as 
transmitted  through  earth's  atmosphere)  for  spectral  regions  extending  respectively 
from  ~0.2  to  3.2ia  and  from  1 to  16|j.  Figure  D-3  also  shows  laboratory  absorption 
spectra  for  H2O  and  CO2  and  some  of  the  other  minor  atmospheric  constituents  listed  in 
Table  D-1. 

D.l  Radiating  Species  for  Shock  Heated  Air 

For  high  temperature  radiation  or  absorption  by  shock  heated  air,  an  investi- 
gation has  demonstrated  that  the  three  most  important  molecular  band  systems  are  those 
historically  designated  as  02-Schumann-Runge  [©2  (S  - R) ] , NO  (B)  and  NO  (y).  In  the 
notation  presented  in  Appendix  B,  these  are  respectively  the  (B  Z X)  system  for  the  O2 
molecule  and  the  (B  Z X)  and  (A  Z X)  systems  for  the  NO  molecule.  Their  -positions 
are  49,163  cm  45,440  cm  ^ and  44,138  cm  ^ respectively. 

It  was  found,  in  spite  of  the  relatively  large  amounts  of  N2  present,  that 
N2-band  systems  made  only  relatively  very  minor  contributions  to  air  radiation  or  air 
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Figure  D-3.  Comparison  of  the  Near- Infrared  Solar  Spectrum  with  Laboratory  Spectra  of 

Various  Atmospheric  Gases 
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opacity  within  the  spectral  region  and  temperature  Interval  of  Interest.  The  N-  Second 
Positive  [N2  (2+)]  with  = 29,670  cm  Is  a type  system  and  the  highly  excited 

lower  electronic  state  Is  only  relatively  sparsely  populated  even  at  8000°K.  As  was 
described  In  Appendix  C,  the  absorption  coefficients  kj^(X,  AX;  T)  actually  are  propor- 
tional to  the  lower  state  populations  rather  than  the  total  number  density 

population  Nj  for  all  states  of  a molecule.  Thus  when  the  electronic  energy  of  the 

lower  state  Is  relatively  high,  as  for  the  B-state  of  N_.  (NV.  /N.)  Is  strongly  temp- 

^ Ijm  1 

perature  dependent  and  there  Is  relatively  little  radiation  or  absorption  until  temp- 
eratures with  kT  equivalent  to  the  excitation  energy  are  approached. 


The  N2  First  Positive  [N2  (1+)  ] system  with  equal  to  9557  cm  ^ Is  the  most 
Important  hot  air  species  radiator  In  the  visible  and  near-infrared  spectral  region 
from  ~ 0.55y  to  ^ ly.  It  Is  a (B-<->-A)  type  system;  Its  upper  state  B Is  the  same  as 
the  lower  state  of  the  N2  (2+)  system.  The  9557  cm  ^ value,  equivalent  to  ~ 1.2 
electron  volts,  thus  also  Is  the  difference  between  the  lower  state  excitations  of  the 
two  systems.  This  appreciable  drop  In  lower  state  excitation  energy  assures  a markedly 
higher  lower  state  population,  N"  = (of  N2)  than  for  the  N"  = Ng  (of  N2)  of  the  N2  (2+) 
system.  However,  the  (N"/N)  ratio  for  N2  (1+)  still  remains  relatively  small  and 
strongly  temperature  dependent  even  at  8000°K. 


High  temperature  air  contains  relatively  small  amounts  of  Ionized  oxygen  atoms 
0 and  NO2  molecules.  It  was  found  however  that  the  electron  attachment  continue  for 
0-atoms,  represented  as  0 (free-bound),  or  0 (fb)  for  short,  and  the  NO2  (fb)  contlnua 
both  made  appreciably  greater  contributions  to  shock  heated  air  opacity  In  the  (P,T)- 
reglme  of  Interest  than  did  either  of  the  N2  (2+)  or  N2  (1+)  systems. 

D.1.1  O2  (S-R),  NO  (B)  and  NO  (y)  Radiation 

For  obvious  reasons,  no  O2  (S-R)  or  NO  (B)  and  NO  (y)  spectra  appear  In  the  solar 
spectra  shown  In  Figures  D-2  and  D-3.  Normal  atmospheric  O2  molecules  are  not  nearly  hot 
enough  to  exhibit  very  much  absorption  at  wavelengths  longer  than  ~ 2400  X,  and  more- 
over the  ©2  (S-R)  absorption  region  Is  masked  by  the  much  stronger  (for  cold  molecules) 
absorption  by  the  high  altitude  0^  layer  In  the  ~ 2100  to  ~ 3000  A region.  Also,  normal 
atmospheric  air  contains  at  most  only  very  minute  proprotlons  of  NO  (see  Table  D-1). 

Figure  D-4,  extracted  from  Reference  [7]  presents  both  computed  curves  for 
radiation  from  Individual  hot  air  species  and  empirical  data  resulting  from  measurements 
of  this  radiation  In  shock  tube  experiments.  The  authors  of  the  referenced  report  [7] 
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plot  spectral  intensity.  In  units  of  [watts  (cm  • ster  • p)  ] for  shocked  air  heated 
to  7000°K  as  a function  of  four  wavelength  Intervals  and  of  air  density 

ratio  (p/p^).  The  predicted  spectral  Intensity  contributions  from  various  air  species 
are  shown  as  solid  lines  In  sub-figures  a,  b,  c and  d,  as  Is  also  the  predicted  total 
air  spectral  intensities  for  corresponding  AX-lntervals . The  empirical  measurements, 
made  by  looking  at  reflected  shock  waves  along  the  shock  tube  axis,  are  shown  by  small 
open  circles.  As  can  be  seen,  the  match  between  theory  and  experiment  is  relatively 
very  good.  The  computations  made  here  at  TRW  Systems,  the  results  of  which  will  be 
presented  In  Appendix  E,  are  based  upon  the  same  equations  and  basic  molecular  parameter 
values  as  were  used  by  the  authors  of  this  reference. 

The  solid  lines  in  the  sub-figures  D-4a,  b,  c and  d present  computed  spectral 
Intensity  contributions  for  NO  (B),  NO  (y),  Oj  (S-R) , N2  (2+),  N2  (1+) , N2'*'  (1-) , and  the 
0 (fb)  [and  the  weaker  0 (free-free)]  continue.  The  dashed  horizontal  lines  represent 
a factor  of  0.1  times  the  corresponding  Black  Body  spectral  intensity  at  7000°K  for  each 
AX-lnterval.  Since  the  effective  radiating  air  sample  thickness  is  given  as  ~ 1.0  cm, 
the  values  corresponding  to  any  (p/p^)  ratio  abslssa  value  and  any  geometrical 

path  L can  be  readily  scaled.  For  any  (p/p^)  abslssa  value  shown,  the  corresponding 
ordinate  values  for  total  air  radiation  and  (1/10)  Black  Body  radiation  are  first  read 
and  then  their  ratio,  r(X,  AX;  7000'^K,  p/p^)  calculated.  Then: 

t(X,  AX;  7000'’K,  p/p^)  w 0.1r(X,  AX;  7000°K,  p/p^)  L (D.l) 

It  follows  directly  from  the  basic  definitions  of  the  mass  density  p,  molecule 
number  density  N,  and  mean  molecular  weight  m,  that; 

(N/N  ) = (p/p  ) • (m  /m)  (D.2) 

o 00 

where  the  o subscript  represents,  as  before,  STP  conditions.  Equation  (D.2)  is 
Independent  of  the  validity  of  the  Ideal  gas  approximation.  When  he  latter  appro?:-- 
imatlon  is  valid,  as  it  is  for  the  (P,  T)-regime  of  current  interest,  a comparison  with 
equation  (1.8)  shows  that: 

(p/p  ) = (P/P  ) • (T„/T)  • (m/m  ) (D.3) 

Unfortunately,  Reference  [7]  gives  no  (P/P^)  or  (m^/m)  data.  It  is  known,  however, 
that  the  (p/p^)  domain  of  current  Interest  lies  somewhere  between  ~ 2.5  and  ~5  for 
T ” 7000°K  (DABS  II-D-1).  Also  from  the  shock  tube  oriented  data  of  Reference  [8],  it 
is  known  that  at  ~ 7000°K  and  (P/P^)  = 113  atm.,  (p/p^)  is  ^ 3.7.  For  p/p^  « 3.7,  the 
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application  of  equation  (D.l)  to  Figures  D-4a  and  b give  x-values  of: 


t(0.25u,  0.02ij;  7000  K,  113  atnO  ~32 
t(0.33m,  0.02y,  7000°K,  113  atm)=r9.6 


The  corresponding  T-values  are  seen  to  drop  below  3 for  the  longer  wavelengths. 


(D.A) 


Figure  D-4  also  shows  that  while  the  N2  (2+)  and  N2  (1-)  radiation  is  dominant 

for  very  low  (p/p^)-ratlos  (l.e.,  at  low  pressures)  for  some  wavelengths  in  the  interval 

of  interest,  that  their  relative  contribution  decreases  with  increasing  p/p  (and  P/P  ) 

0 o 

as  the  mole-fractions  of  NO  and  O2  Increase  relative  to  those  of  N2  (and  of  N2  ) • The 
dominance  of  N2  (1+)  radiation  at  wavelengths  > 0.625y,  and  also  the  implication  of 
relatively  low  t values  is  shown  in  Figure  D-4b. 

D.l. 2 Free  Electron  [e  (ff)]  Radiation 

Free  electrons  in  gas  volumes  both  radiate  and  absorb  radiant  energy  when  their 
momentum  is  altered  by  collisions  with  the  other  gas  particles  (molecules,  atoms,  ions, 
etc.).  The  free-free  [e  (ff)]  absorption  coefficient  generally  is  given  in  terms  of 
actual  geometrical  path  L rather  than  of  atmosphere-centimeters  [which  as  has  been  shown 
previously  for  molecular  absorption,  really  would  represent  (N^  ^ rather  than 

(P/P^)L  . The  geometric  absorption  coefficient  for  e“(ff)  absorption,  to  be  multiplied 
by  L to  obtain  is  given  by  the  expression: 


k'^_«(l/c) 


''e“ 


(D.5) 


0)  + 


Here  v _ , the  collision  frequency  of  the  electrons  with  gas  particles  is  within  two 

® 11  -1 
significant  figures  equal  to  1.0  • 10  (P/Pq)  sec  . As  usual,  c represents  the  speed 

of  light,  equal  to  cgs  units  to  2.997925  • 10^®  cm  sec”^.  The  "plasma  frequency",  co  , 

-1  P 

in  units  of  radians  sec  , is  given  by  the  expression: 

2 


2 -1 
Att  N _ (e)  (m  ) 


(D.6) 


The  symbol  u represents  the  radian  frequency  of  the  optical  radiation  absorbed  or  emitted. 
That  is: 


= 2ii  c X 


-1 


X in  cm 


(D.7) 
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In  the  optical  frequency  and  pressure  regime  of  interest,  i.e.,  2.5  x 10  cm 

_4  2 2 
X ;<  5 X 10  cm  and  (P/P^)  < 200  atmospheres.  It  Is  found  that  (o  >>  and 

equation  (D.5)  simplifies  to  the  form: 


k' 

e~ 


p e (e)  1.0  X 10  sec  /„ /„  n •,  2. 


(c 


(P/P  ) x"n  _ 
o e 


IT  m _ (c)' 
e 


or  for  A in  p[ljj  = 10  cm]  rather  than  in  cm. 


(e)^  1.0  • 10^  (P/P  ) X^(cm^  sec“^)  N _ 

k'_  » ^ — ; X in  w 

e 


(D.8) 


TT  m (c)3 
e~ 


Substituting  numerical  (cgs)  values  for  the  electron  charge,  e~,  the  electron  mass 
and  speed  of  light,  c one  obtains: 

k'_  « 2.99  X lO"^^  (P/P  ) X^  • N _ (cm^  ' y"^) 
e o a 

-3 

X in  (J,  in  cm 

For  the  81  cm  optical  path  length  assumed  for  DABS  II-D-1 


-19  2 3 -2 

T _ = k^_  X 31  cm  « 2.42  x 10  (P/P„)  ^ x N^_  x (cm  x y ) 


(D.8a) 


(D.8b) 


e e 


X in  y,  in  cm 


-3 


(D.9) 


As  an  example.  Reference  [8]  gives  an  N _ concentration  of  6.87  • 10^^  free 

electrons  cm  ^ at  T = 8131°K  and  (P/P^)  = 135.6  (standard)  atmospheres.  By  substituting 

into  equation  (D.9),  one  finds  a value  of  45.7  at  X^  = 4.5y.  Owing  to  the  X^- 

dependence,  however,  the  contribution  at  the  same  temperature  and  pressure  is  only 

0.15  for  X,  = 0.26y  and  0.29  at  X_  = 0.36y.  The  N concentrations  increase  quite 
1 2 e 

strongly  with  both  temperature  and  pressure  [8].  Even  at  the  higher  pressures  and  longer 
wavelengths  of  interest,  the  contributions  become  negligible  for  T much  below  6000°K. 

D.1.3  Infrared  Radiation  from  CO2  and  H2O 

In  the  infrared  spectral  region,  virtually  all  observed  spectral  activity  arises 
from  transitions  between  vibration-rotation  levels  of  the  same  electronic  state,  and 
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mostly  of  the  ground  electronic  state.  The  NASA  Handbook  of  Infrared  Radiation  From 
Combustion  Gases  [4]  gives  very  convenient  spectrally  averaged  absorption  coefficients, 
k(X;  T)  for  the  principal  Infrared  vibration-rotation  bands  of  CO2,  H2O,  CO,  NO,  etc., 
for  selected  temperatures  between  600  and  3000°K.  It  is  readily  seen  from  Figures  D-2 
and  D-3  that  CO2  and  H2O  are  the  most  important  infrared  absorbing  (and  emitting)  species 
in  normal  atmospheric  air.  As  temperatures  Increase,  the  absorption  bands  tend  to 
broaden,  shift  contour  centers  and  intensify  as  transitions  from  successively  higher 
lower  state  levels  become  Important  contributors.  It  was  found,  however,  that  even  the 
strongest  infrared  band,  that  of  the  CO2  fundamental  centered  near  4.3y,  could  give  a 

T-value  ^ 3 for  the  DABS  II-D-1  configuration  only  for  (P/P^)  ^ 120  atm.  This  is  un- 

fortunate since  it  is  desirable  to  have  at  least  one  reliable  near-infrared,  radiometer, 
preferably  near  5u,  for  measurement  of  temperatures  below  4000°K. 

D. 1.3.1  Infrared  Radiation  From  Air  + 4%  CO2 

The  addition  of  4%  CO2  to  the  ambient  air  in  the  shock  tube  (before  the  shock 
occurs)  would  readily  solve  the  preceding  problem.  As  will  be  shown  by  graphs  to  be 
presented  in  Appendix  E,  the  4%  CO2  (which  tends  to  dissociate  to  form  CO  at  temperatures 
above  ^ 3000°K) , would  provide  for  x-values  of  ^ 3 not  only  near  4.3vi  wavelength,  but 
also  within  the  contours  of  the  much  weaker  CO2  "combination  band"  centered  near  2.7p 

wavelength.  Wavelengths  of  X x;  2.8vi  ± O.OSy  and  X^  ± 4.5y  ± G.OSy  were  selected  for 

observations  of  the  air  shock  with  4%  CO2  added.  The  selections  were  Influenced  by  the 
ready  availability  of  relatively  inexpensive  off-the-shelf  filters  (which  may  be  somewaht 
broader,  however,  at  4.5y  wavelength).  The  X^-  centered  at  4.5y  rather  than  4.3y  wave- 
length was  chosen  so  as  to  be  able  to  also  observe  the  high  temperature  CO  radiation 
which  occurs  as  progressively  more  CO2  dissociates  with  increasing  T.  It  will  be  seen 
from  Figure  E-9,  that  the  total  x-value  at  ~ 4.5y  wavelength  should  be  ^ 3 over  virtually 
the  total  (P,  T)  regime  of  interest  from  less  than  1000°K  to  over  8000°K.  The  lower 
temperature  emission  is  from  CO2  and  then  CO2  plus  CO,  and  then  only  CO.  As  T increases 
^6500°K,  the  e (ff)  emission  discussed  in  Section  D.1.3  begins  to  first  supplement, 
and  then  take  over  from  the  eventually  dissociating  CO  molecules. 

As  was  seen  in  Figure  A-3,  the  4.5y  wavelength  radiometer  can  be  combined  not 
only  with  the  X^  = 2.8y  radiometer  for  useful  relatively  low  temperature  measurements, 
but  also  with  the  Xj^  ” 0.26y  radiometer  to  provide  good  coverage  of  the  temperature 
Interval  from  3500  to  above  8000°K.  As  remarked  earlier  in  Section  1.4. 4. 3,  the  upper 
and  lower  temperature  boundaries  are  set  by  practical  rather  than  intrinsic  limitations. 
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D.2  Infrared  Radiation  From  Combustion  Products 


Preliminary  Indications  are  that  the  gaseous  combustion  products  for  the 
DAB  II-D  shots  would  be  Initially  ~ 64%  H2O,  ~ 13%  CO2  and  ^ 22%  N2.  Both  the  0^  and 
CO2  would  begin  to  dissociate  with  temperatures  above  ~ 2500°K  but  Indications  are  that 
there  Is  very  little  possibility  of  exceeding  this  temperature  even  If  considerable 

afterburning  of  solid  explosive  fragments  follows  downstream;  a maximum  temperature  of  ! 

o ^ 

less  than  1000  K Is  much  more  likely.  i 

) 

I 

! 

Radiometers  with  X,  centered  on  1.9p  and  X,  on  4.3p  were  selected  as  the  primary 

J D 

pyrometer  pair  to  view  combustion  products.  As  can  be  seen  from  Figure  D-3,  the  l.Sy 
wavelength  water  vapor  "combination  band"  to  be  viewed  by  the  X^  radiometer  Is  much 
weaker  than  (about  10%  as  strong  as)  the  less  conveniently  located  2.7y  and  6.3p  wave- 
length H2O  "fundamental"  bands.  However,  with  the  65%  of  H2O  present,  the  computed 
t(1.9vi,  O.lp;  T,  P)  values  are  ^ 3 over  virtually  all  of  the  range  Interest,  particularly 
for  T <,  2000°K.  The  2.8y  wavelength  radiometer  channel  also  will  be  available  to  team 
with  the  4.3vi  wavelength  channel  as  a (X_,  X,)  pyrometer  useful  at  the  lower  combustion 
product  temperatures.  The  X^  value  of  4.3vi  rather  than  a repeat  of  the  4.5vi  X^  value 
was  selected  because  It  Is  more  suitable  at  lower  temperatures  If  there  Is  a time  lag 
In  the  arrival  of  the  carbon  soot  particles,  and  because  the  high  temperature  CO 
radiation  at  4.5y  wavelength  Is  Irrelevant  to  the  combustion  products;  Figure  E-lOb 
shows  the  computed  1.9y  and  4.3p  wavelength  t(P,  T) -values  for  the  combustion  products 
for  the  DABS  II-D-1  configuration.  The  computed  x-values  are  very  high  for  the  4.3y 
wavelength  channel  (as  they  also  would  be  at  4.5vi  wavelength  except  for  the  lowest 
temperature  values).  At  the  same  temperature  and  pressure,  x-values  for  the  DABS  II-D-2 
configuration  are  five  times  as  large  as  those  for  DABS  II-D-1. 


APPENDIX  E RESULTS  OF  OPACITY  COMPUTATIONS 


All  the  opacity  computations  presented  In  the  following  sub-sections  are  for  the 
DABS  II-D-1  configuration.  For  the  same  temperature  and  pressure,  the  DABS  II-D-2 
T-values  would  be  a factor  of  5 larger,  owing  to  the  longer  path  length  of  radiating  gas 
being  viewed. 

E.l  Air  Species 

The  opacity  determining  optical  thickness  parameters,  t(A,  AA;  P,  T)  were  computed, 
with  approximate  assumptions  from  equations  (1.6)  or  (1.6a)  for  the  six  selected  wave- 
length intervals  ± 1/2AX^.  The  required  thermal  equilibrium  mole-fraction  values, 

(T,  P),  were  computed  for  temperatures  up  to  5000°K  from  an  available  in-house  program. 
Unfortunately,  this  program  has  not  yet  been  extended  to  Include  chemistry  reaction  rates 
for  temperatures  above  5000°K.  For  temperatures  above  5000°K  and  for  P/P^^  > 10  atms. 
species  number  density  values  N^  (T,  P) , rather  than  mole-fraction  values,  were  extracted 
from  extensive  computed  theoretical  "shock  tube  chemistry"  tables  given  In  a JPL  Report  [7]. 
These  tables  listed  N^  - values  corresponding  to  different  presumed  temperature  and 
pressure  initial  conditions  for  an  "air"  shock,  and  for  a presumed  resulting  "moving 
shock",  "standing  shock",  and  "reflected  shock".  The  temperatures,  pressure  ratios  (P/P^) , 

density  ratios  (p/p  ),  enthalpies,  etc.,  are  listed  for  each  N. -column.  However, 
o 1 

systematic  distributions  of  temperature  values,  or  of  (P/P^) -values  for  any  temperature 
are  not  available,  and  the  desired  data  must  be  gleaned  by  checking  through  each  set  of 
tables  to  find  If  any  (T,  P/Pp)  pairs  of  Interest  are  Included.  Although  the  data 
extracted  by  this  gleaning  process  was  appreciably  less  detailed  than  was  desirable,  there 
was  still  sufficient  data  to  make  useful  interpolations  and  extrapolations. 

Figure  E-1  presents  the  resulting  computed  mole-fraction  data  for  N2,  O2,  NO, 

CO2  and  CO  in  air  for  T varying  from  1000  to  ~8000°K  and  for  pressure  intervals  of  10 
atms  £ — 125  atms.  The  10  atmosphere  data  was  taken  directly  from  the  Langley 

Research  Center  NASA  Technical  Note  by  J.  S.  Lee  and  P.  J.  Bobbitt  from  which  Figure  D-1 
was  extracted  [6].  It  is  seen  that  the  dissociation  of  N2,  O2  and  CO2,  and  the  consequent 
formation  of  NO  and  CO,  tends  to  be  lowered  as  pressure  Increases.  Thus  0^,  N2  and 
CO2  radiation  tends  to  Increase  with  increasing  pressure.  At  the  lower  temperatures,  the 
reverse  is  true  for  CO  and  NO  radiation.  As  temperature  Increases,  however,  dissociation 
also  tends  to  occur  in  NO  and  CO  molecules,  and  this  in  turn  is  slowed  by  increase  in 
pressure. 


Mean  spectral  absorption  coefficients  for  Oj  (S-R),  NO  (6),  NO  (y),  N-  (2+) , 
N2  (1+),  0 (ff),  the  NO2  continua,  and  N2  (1-)  were  computed  from  formulae  and  data 


MOLE-FRACTIONS,  X,  IN  AIR 


Figure  E-1.  Temperature-Pressure  Dependence  of  Mole  Fractions  In  Air 
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given  (mostly  in  Chapter  V.  121  of  Volume  I)  in  the  excellent  text  by  Zel'dovich  and 

Ralzer  [8],  The  reader  should  be  warned,  however,  that  there  is  a confusion  in  definition 

of  their  c^  "concentration"  parameter  occurring  on  the  formulae  given  on  page  334  (and 

elsewhere)  of  this  reference.  lor  physical  consistency,  their  c^  cannot  possibly  be  equal 

to  (N^/N^)  as  appears  to  be  implied  by  the  somewhat  obscure  wording  of  their  definition, 

perhaps  owing  to  a mls-translation  from  the  original  Russian  text.  Rather,  their  entire 

(c.  p/p  ) term  must  be  replaced  by  either  (N  /N  ) or  by  its  equivalent 
10  g 1 o 

(Xj^  N/N^)  = Appendix  D].  The  data  presented  by  Zel'dovich  and  Raizer 

was  largely  taken  direc?ly  or  adapted  from  the  previously  referenced  paper  [7]  from  which 
Figure  D-4  was  extracted. 

Figure  E-2  presents  the  absorption  coefficients  for  O2  (S-R) , NO  (6)  and  NO  (y) 
contributions  computed  from  the  Zel'dovich  and  Ralzer  equations  for  temperatures  of 
10,000,  8000,  6000,  and  4000°K.  It  is  seen  that  the  O2  (S-R)  absorption  (or  emission) 
is  the  strongest,  and  that  the  NO  (y)  radiation,  is  not  only  weaker  than  the  NO  (B) 
near  0.2600iJ,  but  drops  off  much  more  rapidly  with  increasing  wavelength. 

Figure  E-3  presents  similar  data,  but  with  5000  rather  than  4000  “K  curves,  for 
the  N2(2+)  and  N^  (1-)  band  systems.  It  Is  seen  that  the  N2  (2+)  absorption  is  more 
than  two  orders  of  magnitude  weaker  than  that  for  NO  (B)  near  8000  “K,  and  drops  off 
very  rapidly  with  decreasing  temperature.  Therefore,  in  spite  of  the  relatively  large 
N2  concentrations,  as  shown  in  Figure  E-1,  the  N2  (2+)  contributions  to  the  computed 
t's  were  very  minor.  As  previously  mentioned  (see  Figure  D-4),  the  N2  (1+)  spectral 
activity  occurs  mostly  in  the  longer  wavelength  visible  region  and  turned  out  to  be 
too  low  to  give  t's  > 1 for  the  DABS  II-D-1  configuration. 

The  N2  (1-)  absorption  coefficients  are  relatively  very  large  as  is  seen  (on 
the  right  hand  ordinate  scale)  in  Figure  E-3.  In  the  (T,  P)  regime  of  interest,  however, 
the  number  density  of  N^  ions  was  insufficient  to  provide  more  than  very  minor  contri- 
butions to  t(0.36p  , O.Oly;  T,  P). 

Absorption  coefficients  for  the  0 (fb)  and  NO2  contlnua  are  not  shown  since 
their  t-  contributions  were  computed  directly  without  intervening  x^  ^ computations 
from  relatively  very  simple  formulae  and  data  given  by  Zel'dovich  and  Ralzer. 

Figure  E-4  presents  the  results  of  opacity  calculations  for  2gp 

O.Olp;  T,  P)  and  t^^^(0.36vi,  O.Olp;  T,  P) . As  previously  stated,  the  contributions 
were  mostly  from  O2  (S-R),  NO  (6)  and  NO  (y)  at  the  shorter  wavelength,  and  from  the  first 
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two  of  these  for  0.36ii  wavelength.  The  NO2  and  0 (fb)  T-contributlons  added  up  to 
as  large  as  ~8%  for  some  temperatures  and  pressures.  It  is  seen  that  the  0.26u  wave- 
length radiometer  will  be  viewing  the  equivalent  of  virtually  block  body  radiation 
throughout  the  full  range  of  the  (T-P)  regime  of  interest  for  both  Station  1 and  Station 
2 of  DABS  II-D-1.  This  is  not  completely  true  however,  of  the  0.36vi  wavelength  radio- 
meter for  Station  2,  where  the  air  shock  is  predicted  to  arrive  at  -<82  atm.  (1200 
psla)  and  decay  to  ~41  atms.  (600  psla)  in  1 msec.  Figures  1-3  and  1-4,  and  Table 
1.1  show  that  the  corresponding  predicted  air  shock  temperature  interval  is  from~5100  “K 
to  ~ 4100  °K.  Figure  E-4  shows  that  a pressure  of  ~50  atms,  is  required  to  keep 
T(0.36y,  0.02u;  T,  P)  above  ~ 3 near  4100  “K.  [t's  5 x greater  for  DABS  II-D-2] 

E.2  Air  + 4%  CO2 

Figure  E-5  shows  that  temperature  and  pressure  dependent  mole-fractions  of  CO2 
and  CO  computed  for  4%  CO2  added  to  the  air.  The  computation  for  these  cuirves  was 
done  on  the  in-house  TRW  program  which  is  useful  for  T ^ 5000  °K,  with  small  extra- 
polations to  higher  temperatures  being  guided  by  the  10  atms  data  of  Reference  [6]. 

Figure  E-6  gives  k(X,  AX;  T)  values  for  CO2,  CO  and  H2O  (for  later  combustion  gases 
use)  taken  from  graphs  and  tables  given  by  Reference  [4]. 

Figures  E-7  and  E-8  show  the  respective  values  computed  for  t(2.8u,  O.lp;  T,  P) 
and  t(4.5m,  O.ly;  T,  P) . These  plotted  curves  are  not  critically  dependent  upon  the 
radiometer  band  centers  2.8u  and  X^  « 4.5p  , nor  upon  the  assumed  AX  Intervals 

of  O.ly  as  long  as  the  spreads  X + 1/2  AX,  remain  within  desired  intervals.  It  is 
seen  that  the  X^,  4.5  y,  radiometer  should  be  looking  at  equivalent  Black  Body  radiation 
for  virtually  all  pressure  values  of  Interest  and  for  temperatures  from  —1000  to  > 6000“K. 
The  1000  "K  cutoff  occurs  because  of  the  selection  of  X^  =s  4.5y  rather  than  4.3y  wave- 
length. The  Xy  2.8y  , radiometer,  however  begins  to  lose  its  usefulness  for  T ^ 3500  ®K 
at  the  lower  pressure  values.  This  is  not  particularly  serious  because  the  (X^,  X^) 
pair  of  radiometers  is  of  primary  usefulness  for  T ^ 3500  “K. 

The  T^-value  shown  in  Figure  E-8  are  not  complete,  since  they  do  not  Include 

e“(ff)  contributions,  as  computed  from  equation  (D.9).  Figure  E-9  shows  superimposed 

plots  of  Tj^(0.26y,  O.Oly;  T,  P)  and  T^(4.5y,  O.ly;  T,  P)  extending  from  —1000  “K  to 

—8000  “K,  with  the  e“  absorption  contributions  now  Included  in  the  latter  curves.  The 

upper  and  lower  (P/P  ) pressure  boundaries  for  the  t-reglons  shown  are  125  and  30  atms. 
o 

respectively.  It  readily  is  seen  that  the  (X^^,  X^)  pair  both  are  looking  at  effectively 
"Black  Body"  radiation  over  virtually  the  full  (P,  T)  regime  of  interest.  However, 
signal  fall  off  and  transmission  noise  and  dynamic  range  considerations,  and  the  non- 
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Figure  E-7.  Computed  Optical  Thickness  Near  2.8;j  Wavelength  for  CO^  In  Air 
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E-8.  Computed  Optical  Thickness  for  CO 
to  Air  - A.5y  CO-  + CO  Absorption  Region 
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Figure  E-9.  Computed  Optical  Thickness  at 
).26u  and  at  4.5u  for  A%  CO.  Added  to  Air 


expendabillty  of  higher  priced  optimum  detectors  and  low  noise  amplifiers,  etc.,  would 
tend  to  restrict  the  useful  field  operations  range  of  the  radiometers  to  T ^ 3000  °K. 

E.3  Combustion  Products 

Figures  E-lOa  shows  computed  values  of  mole-fractions  of  CO2  and  of  H2O  in  the 
combustion  products.  The  corresponding  k(X,  AX;  T)  values  were  shown  in  Figure  E-6. 
Figure  E-lOb  shows  the  computed  values  of  T^(1.9p,  O.lp;  T,  P)  and  t^(4.3u,  O.lu;  T,  P) . 
The  cross-hatched  portion  at  the  lower  left  of  the  1.9p  wavelength  region  Indicates 
an  "area"  where  the  low  density  spectral  "line  packing"  in  low  temperature  H^O  creates 
some  doubt  of  reliability.  However,  the  arrival  of  carbon  soot  following  closely  be- 
hind the  air  combustion  gas  contact  surface  would  soon  insure  virtually  complete  opacity. 
It  now  appears  that  the  carbon  soot  lag  time  is  negligible. 


The  temperature  limits  for  the  combustion  products,  assuming  no  afterburning,  is 
not  predicted  to  rise  much  above  800  °K  at  most.  However,  these  predictions  are  for 
Isentropic  expansion,  and  do  not  Include  afterburning,  which  possibly  could  raise  these 
temperatures  to  ~2500  “K.  For  this  reason.  Figure  E-10  calculations  were  carried  out 
to  3000  “K. 

E.3.1  Carbon  Soot 

Figure  E-11  presents  empirical  near  infrared  spectral  absorption  coefficients 
k^  for  carbon  particle  soot  at  temperatures  ranging  from  300  to  2600  “K.  The  figure 
was  extracted  from  the  NASA  Handbook  of  Infrared  Radiation  From  Combustion  Gases  (page 
228)  [4].  As  is  indicated  by  its  title,  the  ordinates  represent  mass-density  rather 
than  number  density  coefficients.  In  order  to  obtain  the  desired  values,  therefore, 
these  ordinates  must  be  multiplied  by  the  product  of  optical  path  length  L with  the 
mass-density,  p^,  of  carbon  soot,  rather  than  by  (N^/N^)  ■ L as  was  true  for  the  molecular 
gas  absorption  coefficients.  This  is  shown  in  equation  E.l  below. 
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Figure  E-10.  Combustion  Products  Curves 


Figure  E-12  presents  the  computed  values  for  the  carbon  soot  from ~ 500  to 
2500  °K  for  » 1.9p  and  X^  »4.3vi.  These  curves  can  be  extrapolated  somewhat  to 
higher  and  lower  temperatures  and  also  to  neighboring  wavelength  Intervals,  as  Is  se^n 
from  an  Inspection  of  Figure  E-11.  It  should  be  noted  that  the  bottoms  of  the  hatched-in 
areas  go  down  to  (P/P^)  = 1 atmosphere  rather  than  only  to  10  or  30  or  40  atmospheres 
as  for  gaseous  t-  contributions. 

After  the  carbon  soot  particles  arrive  at  an  observation  station,  now  presumably 
with  a negligible  lag  time  behind  the  air  shock  combustion  gases  contact  surface,  it 
would  appear  that  the  "Black  Body"  condition  of  e(X,  AX^;  P,  T)  « 1 would  certainly  be 
satisfied  for  the  near  infrared  radiometer  AX  intervals.  Very  probably,  from  the 
appearance  of  Figures  E-10  this  would  also  be  true  at  0.26u  and  X^^O.lbp.  At 

the  expected  low  combustion  product  temperatures  however,  the  signals  from  these  latter 
two  radiometers  probably  would  be  too  weak  to  utilize. 

It  would  appear  that  streaming  of  cold  carbon  soot  into  the  observation  port 
could  cause  a serious  problem,  not  only  in  data  interpretation,  but  also  in  data  ob- 
servation if  the  lenses  become  coated  with  soot.  It  would  therefore  appear  imperative 
to  combat  this  possible  carbon  soot  streaming  and  coating  by  port  design  and  also 
possibly  with  a counter  stream  of  dry  nitrogen  gas. 
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Figure  E-12.  Computed  Optical  Thickness,  x . of  Carbon  Partical  Soot  Carried  by  Combustion  Gases 
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NOMENCLATURE 


A)  By  Cy  ••• 


A 


c 


C 

c 


^2 


CO 

CO2 

D 

d,  d^j(X;  T) 
d,  d^^(X,  AX;  T) 


e 


e 

e-(ff) 


eV 


the  first,  second,  third,  etc.,  electronically  excited 
quantum  states  of  a molecule  or  molecular  ion,  in  order 
of  Increasing  excitation  energy  above  the  ground  electronic 
state  (represented  by  the  symbol  X) 

the  Angstrom  unit  of  spectral  wavelength  measure, 
ik  = 10~®  cm 

the  responsive  area  of  the  detector  cells,  a square, 

0.4  cm  on  a side,  or  0.16  cm^,  for  all  five  detectors 
used  in  DABS  II-D-2 


the  carbon  atom 

the  speed  of  light  {2.997925  • 10^*^  cm/sec] 

the  'first  radiation  constant',  here  (hut  not  universally) 
defined  as  Cj^  “ 2 c^  h ■ 1.1910  • 10”  watts  cm^  ster"^ 

the  second  radiation  constant,  universally  defined  as 
- (hc/k)  - 1.43880  (cm  “K)  - 1.43880  • 10^  • (y  “K) 

the  carbon  monoxide  molecule 


the  carbon  dioxide  molecule 

the  diameters  of  the  lenses  » 5.08  cm  (2  inches) 

the  strongly  temperature  dependent  distance,  in  cm  ^ 
(wave^^^mber)  units,  between  adjacent  spectral  lines  (for 
the  j band  system  of  the  1^  species) 


the  spectral  mean  value  of  d over  a X + 1/2  AX  interval 

pyrometer  error  factors  (mostly  "measurement  error 
attenuators").  See  Appendix  A,  equations  (A. 15)  and 
(A. 19a,  b,  c). 


the  (negative)  charge  on  an  electron  « 4.80298  • 10 


the  electron 


the  electron  "free-free"  absorption  and  emission  spectral 
continuum,  important  for  high  temperatures  and  long 
wavelengths 

-12 

an  electron  volt.  1 eV  1.60210  • 10  erg,  and  is 
equivalent  to  1.16049  • 10  “K 
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F 

(ff)  ' 


(fb) 


»2° 

1°,  I°(X;T) 


I,  I(X;  P,T;L) 


1°,  I°(X^;  T) 

I^,  i°(X^,  AX^^;  T) 


J, 


J’, 


J" 


*K 


k 


k^j(X;  T) 


k’ 


e 


L 


n 

o 


the  focal  length  of  the  lenses  = 50.8  cm  (20  inches  for 
each  lens  used ) 

Indicating  a free-free  spectral  transition  during  non- 
capture, but  inelastic,  collisions  between  electrons  and 
atoms,  molecules  or  ions,  or  between  molecules  and  ions, 
etc. 

an  ionization,  dissociation  or  capture  transition  lead- 
ing to  a spectral  continuum 

the  amplifier  gain  factor  for  the  l*”^  radiometer  system 

-27 

Planck  Constant  = 6.6256  • 10  ergs 
the  water  vapor  molecule 

Black  Body  spectral  radiance,  in  units  of  watts 
[cm^  ster  u]“  or  the  equivalent 

spectral  radlsince  from  a surface  element  of  a layer  of 
hot  gases,  in  the  same  units  as  for  I . 

I = 1°  e(X;  P,  T;  L) 

Black  body  spectral  radiance  for  a specific  set  of 

X,  , X., , ...  X values 
1 z n 

Black  Body  spectral  radiance  mean  value  averaged  over 
the  interval  X + 1/2  AX^  for  a specific  set  of  X^^, 

X_,  . . . , X values 

n J 

the  principal  rotational  quantum  number  of  a molecule, 
ion  or  atom.  J'  and  J"  correlate  to  the  respective 
upper  and  lower  quantum  states  in  a spectral  transition 

degrees  Kelvin;  0 °K  = -273.15  °C 

the  Boltzmann  Constant  « 1.38054  • 10  erg  (“K)  ^ 

the  "spectral  absorption  coefficient"  in  units  of 
inverse  atmosphere-centimeters  as  conceptually  reduced 
to  STP  conditions  [(atm  cm  - STP)”^].  See  Appendix  C, 

Section  C.l  for  clarification  of  the  meaning.  The  li 
dummy  indices  refer  to  the  band  system  of  the  1^ 
spectrally  active  species. 

the  spectral  mean  value  of  the  above  parameter,  as 
averaged  over  a X + 1/2  AX  spectral  interval 

the  (geometric)  absorption  coefficient  for  e~(ff)  • 

spectral  activity  ] 

optical  (geometrical)  path  length  through  an  absorbing  I 

and  emitting  media  (here  a layer  of  gas)  j 

! 

mean  molecular  weight  of  gas  particles  at  STP  conditions  j 
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m _ 
e 

m,  ni(P,  T) 

N 

N,  N(P,  T) 

N^,  N^(P,  T) 
N 

o 


Nljk 


N , N (P,  T) 
e 


N2(-1+) 


N2(2+) 


N^d-) 


NO 


—28 

electron  rest  mass  « 9.1091  • 10  gm 

molecular  (or  atomic)  weight  of  the  l*"^  species  of 
spectrally  active  molecule,  atom  or  Ion 

mean  molecular  weight  of  gas  particles  at  pressure  P 
and  temperature  T 

the  nitrogen  atom 


number  density  of 
cm~^  [l.e.,  total 
per  cm^] 


gas  species  particles,  In  units 
number  of  particles,  of  all  gas 


of 

species. 


number  density  of 


species  particles. 


In  units 


of  cm 


-3 


number  density  of  gas  particles  at  STP  conditions 


N = 2.67  • 10^^  cm' 


the  number  densities  of  1 species  particles 
the  respective  upper  and  lower  quantum  states 
the  spectral  line  of  the  band  system. 
Appendix  C,  Section  C.l] 


excited  to 
to  generate 
[See 


number  density  of  electrons.  In  units  of  cm 


the  nitrogen  molecule 


the  singly  Ionized  nitrogen  molecule 

the  "first  positive"  spectral  band  system  of  the  N2 
molecule,  corresponding  to  a (B  A)-  type  electronic 
transition  and  with  Its  "origin"  near  l.OSy 

the  ’second  positive'  spectral  band  system  of  the  N2 
molecule,  corresponding  to  a (C  B)-type  electronic 
transition,  and  with  Its  "origin"  near  0.3370y 

the  "first  negative"  spectral  band  system  of  the  molecular 
nitrogen  Ion  N2,  corresponding  to  a (B  ^ X)=type 
electronic  transition,  and  with  Its  "origin"  near  0.391y 


the  nitric  oxide  molecule 


NO  (3) 


NO  (y) 


the  "beta-band  system"  of  the  nitric  oxide  molecule, 
corresponding  to  a (B  ^ X)-type  electronic  transition, 
and  with  Its  "origin"  near  0.2201y 

the  "gamma-band  system"  of  the  nitric  oxide  molecule, 
corresponding  to  an  (A  t-  X)-type  electronic  transition, 
and  with  Its  "origin"  near  0.2266y 

the  nitrogen  dioxide  molecule 
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N02(fb) 


the  spectrally  very  broad  nitrogen  dioxide  free-bound 
continuum 


n'  and  n" 


0 

O" 

0"(fb) 
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(P/P  ) 
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Input 

^ j 
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ij 
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6r 


S 
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T 

T 


the  principal  electronic  quantum  numbers  correlated  to  the 
respective  upper  and  lower  quantum  states  of  a spectral 
transition 

the  oxygen  atom 

the  negative  oxygen  ion,  with  a single  additional  attached 
electron 

the  broad  free-bound  spectral  continuum  of  the  negative 
oxygen  atomic  ion 

gas  pressure,  in  units  of  standard  atmospheres  or  pounds 
per  square  inch,  absolute  (psia) 

one  standard  (sea  level)  atmosphere  = 760  mm  Hg  = 14.7 
psia 

gas  pressure  in  units  of  standard  atmospheres 

shock  front  pressure  (at  time  of  initial  arrival  of  any 
specified  observation  station) 

predicted  isentroplc  decay  pressure  at  some  time  after 
passage  of  shock  front  (Section  1.3) 

combustion  products  pressure 

the  effective  radiant  power  input  to  the  i*"^  radiometer 
system  (See  equation  3.1) 

the  response  factor  of  the  i^'^  radiometer  cell,  in  amps 
output,  per  watt  input  for  = 0.26v  and  0.36vi  , and 
in  volts  output  per  watt  input  for  Xj^  = 2.8u,  4.5y  and 
1.9y 

the  ratio,  (I  /I  ) « for  both  and  e near 

unity.  This  R. .^ratio  is  the  basic  "temperature  signal 
input  for  a two-color  pyrometer  system 

the  distance  from  a lens  to  the  viewing  port  aperture 

the  distance  of  a viewed  radiating  molecule  behind  the 
viewing  port  aperture 

the  incremental  projected  area  of  the  radiating  gas  layer 
actually  focused  within  the  responsive  area  AC  of  a 
detector  cell  [See  Section  3) 

gas  temperature  in  units  of  °K 

standard  temperature  * 273.15  °K 
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6t 


At 


(v^.  V2.  V3...,  v^) 


;rd  v" 


shock  front  temperature  (at  time  of  initial  arrival  of 
any  specified  observation  station;  corresponds  to  P^) 

predicted  Isentroplc  decay  temperature  at  some  time 
after  passage  of  shock  front  (Section  1.3,  corresponds 
to  P^) 

time 

time  Interval  between  arrival  of  shocked  air  front  and 
combustion  products  contact  surface  at  any  station 

available  measurement  time  Interval  between  initial 
arrival  of  shocked  air  front  and  destruction  of  measuring 
apparatus  at  any  station  [Figure  2-5] 

the  geometrical  optics  object  distance  with  weighted 
mean  value  of  (r  + 6r)  varying  slightly  for  each  radiometer 

the  geometrical  optics  image  distance  (in  practice,  the 
adjustable  set  distance  between  a detector  cell  and  its 
lens) 

the  i^^  radiometer  output  signal  in  units  of  volts  for 
the  infrared  radiometers  and  units  of  amperes  for  the 
0.26vi  and  0. 36p  radiometers 

the  vibrational  quantum  number  for  a diatomic  molecule 
(or  diatomic  molecular  ion) . v takes  on  the  values 
0,  1,  2,  3 ...  depending  upon  the  degree  of  vibrational 
excitation 

a set  of  vibrational  quantum  numbers  corresponding  to 
each  of  the  m normal  vibrational  modes,  of  a polyatomic 
molecule  (or  molecular  ion).  Each  v.  can  assume  the 


values  0,  1,  2,  3, 


depending  upon  the  degree  of 


excitation  of  the  corresponding  normal  mode  within  an 
arbitrary  combinational  mode  [see  Appendix  B,  Section  B.l] 

the  vibrational  quantum  numbers  corresponding  to  the 
respective  upper  and  lower  excitation  states  for  a spectral 
transition  in  a diatomic  molecule 


(v,',  v' . . .v'),  (v",  vL',...v")  the  sets  of  vibrational  quantum  numbers  corresponding  to 
X ^ in  X ^ m . ” ~ , , - • , g. 

the  respective  upper  and  lower  excitation  states  for  a 

spectral  transition  in  a polyatomic  molecule 

X,  x^  the  analytically  and  computationally  useful  dimension- 

less parameter  equal  to  c.^l\T  or  02/^^  T 
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a(A;  P,  T;  L) 


a(X,  AA;  P,  T;  L) 

b[(x) 

glens 

Y 


Y,  Yij  (P.  T) 


(Y/d) 


(y/d)  or  (y/d) 

A,  6 

e>  e(A5  P»  Tj  L) 

e,  e(X,  AX;  P,  T;  L) 

X, 

AX,  AX^ 


A'X 


i 


spectral  absorptivity  of  a layer  of  absorbing  gases  of 
thickness  L,  temperature  T and  pressure  P.  The  a-para- 
meter  also  depends  very  strongly  on  the  proportions  of 
different  gas  species  In  the  layer 

the  spectral  mean  value  of  the  preceding  a-parameter 
over  a X + 1/2AX  spectral  Interval 

the  spectral  filter  transmissivity  at  wavelength  X,  for 
the  arbitrary  1^^*  radiometer  filter  used  on  DABS  II-D-2 

the  spectrally  flat  transmissivity  of  the  lenses  used 
on  DABS  II-D-2.  glens  r-  0.90  for  each  of  the  five  lenses 

ratio  of  specific  heat  at  constant  pressure  to  specific 

heat  at  constant  volume;  y = C /C 

P V 

the  value  of  y at  STP  conditions 

the  half  width  at  half  peak.  In  units  of  cm  of  a 
(conceptually)  Isolated  spectral  "line".  See  Appendix  C, 
Section  C.2,  and  equation  (C.3)  and  adjacent  discussion 
for  P-  and  T-dependence  of  y 

the  strongly  temperature  and  pressure  dependent  dimen- 
sionless "packing  factor",  or  "packing  ratio",  of 
spectral  lines 

the  spectral  mean  value,  over  a suitable  X + 1/2  AX 
Interval,  of  the  above  "packing  factor" 

Indicating  Incremental  values 

spectral  emlsslvlty  of  a layer  of  emitting  gases  of 
thickness  L,  temperature  T,  and  pressure  P (also  very 
strongly  dependent  upon  proportions  of  different  gas 
species  In  the  layer).  From  Klrchoff's  Law,  e(X;  P,  T;  L) 
- a(X;  P,  T;  L) 

the  spectral  mean  value  of  the  preceding  e-  parameter 
over  a X + 1/2  AX  Interval 

optical  radiation  wavelength,  measured  In  units  of 
microns  (p),  angstrom  units  (A),  or  centimenters  (cm) 

The  dummy  subscript  1 Indicates  one  of  a set  of  specific 
wavelengths 

Incremental  wavelength  Intervals  centered  on  X or  X^, 
usually  referring  to  optical  filter  bandpass  or  to 
spectral  regions  of  strong  absorption  or  emission. 

the  effective  width,  somewhat  larger  than  AX^,  for  a 
spectral  filter  bandpass  that  falls  to  decay  as 
rapidly  as  Is  desirable. 
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V,  V. 


''e- 


Av,  Av. 


P,  P(P,  T) 
P, 


o 

T,  T 


t^j(A,AX;  P,  T;  L) 


T^jf  AX,  X^P » 


Xi..  Xi(P,  T) 


P 

c 


n - 8(D,  r) 


a = n(D,  r + 6r) 


primes,  ' and  " 


, 4 * -A 

one  micron  (or  micrometer.  Ip  « 10  A - 10  cm  - 

10"®  m 

optical  radiation  wavenumber  spectral  positions.  When 
X Is  measured  In  cm  units,  then  v » X“1  and 

the  "origin"  of  an  electronic  band  system 

the  collision  frequency  of  electrons  with  gas  particles. 
In  units  of  sec 

Incremental  wavenimiber  Increments  centered  on  v or 

_3 

gas  density.  In  units  of  gm  cm 

gas  density  at,  or  reduced  to,  STP  conditions 

optical  thickness  and  spectral  mean  value  of  optical 
thickness 


contribution  to  spectral  mean  optical  thickness  by  j 
band  system  of  !*■  species  gas  particle 


th 


the  spectrally  averaged  optical  thickness  contribution, 
in  the  X + 1/2  AX  spectral  in^rval,  from  the 
spectral  band  system  of  the  i'"  species  gas  particles. 

The  x^P  product  in  the  abstract  functional  representation 
of  can  be  replaced  by  the  corresponding  species 
partial  pressure  P^,  or  species  number  density 

the  mole  fraction  of  1^^  species  of  gas  particles 
(molecules,  atoms  or  Ions) 

the  "optical  frequency"  of  spectral  radiation  (used 
only  in  Appendix  D,  Section  D.1.2  for  e“(ff)  radiation), 
with  u)  “ 2ir  c X 

the  "plasma  frequency",  again  used  only  in  Appendix  D, 
Section  D.1.2  [see  equation  D.6)] 

c c ■*2 

the  detector  cell  solid  angle  field  of  view,  u • A v , 
where  v is  the  geometrical  optics  "image  distance", 
l.e.,  the  distance  between  the  lens  and  the  detector  cell 

the  solid  angle  subtended  by  a lens  f^om  ^ny  point  on 
the  viewing  port  aperture,  fj  “ (irD  /4r  ) w 2.134 
* 10  steradians  for  each  lens  used. 

the  solid  angle  subtended  by  a lens  from  a radiating 
gas  molecule  a distance  6r  behind  the  port  aperture. 
n(D,  r + 6r)  - [itD^/4(r  + 6r)^] 

single  and  double  primes,  ' and  ",  are  used  to  indicate 
parameters  belonging  to  the  respective  "upper"  and 
"lower"  quantum  states  in  any  "bound-bound"  spectral 
transition 
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PART  B 


NOVEL  MEASUREMENT  TECHNIQUES 
1.  INTRODUCTION 

1.1  Objectives 

The  objectives  of  the  theoretical  study  described  herein  were  to  establish  the 
conceptual  feasibility  of  deploying  novel  instrumentation  to  measure  temperature 
histories  and  density  histories  within  the  interior  of  a nuclear  driven  shock 
tunnel.  This  study  will  only  address  the  basic  physics  concepts  rather  than  present 
hardware  design  solutions.  Specifically,  it  was  considered  desirable  to  investigate 
procedural  techniques  to  measure  the  density  and  temperature  history  of  the  shocked 
air  before,  during  and  after  the  passage  of  the  shock  front,  and  the  density  history 
of  any  lofted  debris.  It  is  virtually  certain  that  the  tunnel  walls  will  undergo 
transitory  displacement  during  the  passage  of  the  shock.  Since  this  displacement,  if 
not  considered,  could  readily  distort  the  interpretation  of  the  data,  it  is  desirable 
to  have  measurements  of  the  local  wall  displacement  history.  To  sum  up,  the  objectives 
of  these  proposed  novel  measurements  are  to  obtain: 

• Density  and  temperature  histories  of  shocked  air 

• Density  history  of  lofted  debris 

• Translational  displacement  history  of  tunnel  walls  in  the  tunnel  radial  direction 

1.2  Background 

These  types  of  measurements  have  not  been  made  previously.  During  the  Marvel 
experiments  [1,  2 and  3],  a series  of  measurements  were  made  for  various  parameters 
associated  with  the  shock  generated  by  the  detonation  of  a 2.2  KT  nuclear  explosive  at 
one  end  of  a one  meter  diameter  horizontal  tunnel  122  meters  long.  The  only  experi- 
ments looking  directly  along  the  axis  into  the  shocked  air  viewed  the  advancing  shock 
front  with  luminosity  instrumentation  placed  at  the  far  end  of  the  tunnel.  In  addition, 
a series  of  light-pipe  photodetectors  were  placed  in  side  alcoves  to  measure  time  of 
arrival  of  the  shock  front.  The  former  instrumentation  saw  mostly  the  precursor  radiation 
emitted  by  the  optically  dense  radiation  heated  air  well  in  advance  of  the  actual  shock 
front.  Owing  to  the  extremely  high  temperature  and  pressure  air  in  and  behind  the 
shock  front,  the  time  of  arrival  photodetectors  only  responded  to  the  radiation  emitted 
within  a few  centimeters  of  their  immediate  vicinity.  Thus  no  real  measurements  were 
made  of  state  parameters  within  the  interior  of  the  shocked  air. 

To  date,  measurements  of  the  type  proposed  in  this  report  have  not  been  accurately 
made  in  such  extreme  environments  (temperatures  above  12,000  °K  and  pressures  of  ~1000 
atmospheres) . 
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1.3  Proposed  Techniques 
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This  report  will  discuss  two  techniques  which  will  provide  the  means  for  obtaining 
the  needed  measurement  histories  presented  above.  The  proposed  techniques  are: 

• Gas  bag  method 

• Radioisotope  method 

The  first  technique,  the  "gas  bag  method",  is  not  a measurement  technique  per  se, 
but  a method  for  providing  an  optical  window  into  the  interior  region  of  the  shock  front. 
This  would  permit  the  use  of  passive  two  color  (or  multi-color)  optical  pyrometers 
for  determining  temperature  histories  of  the  shocked  air. 

The  second  technique,  the  "radioisotope  method"  is  an  active  X-ray  system  using 
two  or  more  different  frequency  radioactive  X-ray  sources  with  detectors  to  measure 
areal  density  histories  and  tunnel  wall  translational  displacement.  The  density  history 
would  consist  of  measurements  of  the  initially  arriving  shocked  air  and  the  later 
arrival  of  shocked  air  plus  lofted  debris.  These  measurements  along  with  the  wall 
displacement  histories  will  be  used  to  resolve  the  separate  areal  density  contributions 
of  air  and  lofted  debris. 

1.3.1  The  Gas  Bag  and  Pyrometer  Technique 

The  concept  of  two  color  or  multi-color  optical  pyrometers  has  been  fully 
discussed  and  their  usefulness  demonstrated  in  Part  A of  this  report.  In  a nuclear 
driven  shock  tube  experiment,  however,  temperatures  and  pressures  are  both  much  higher 
than  in  the  DABS  II-D  experiments,  and  more  hostile  to  Instrumentation  system  survival 
and  performance.  Moreover,  the  higher  temperatures  and  pressures  result  in  greatly 
Increased  radiation  thus  causing  even  a few  centimeter  layer  of  shocked  air  to  be 
completely  opaque.  Also  precursor  radiation  tends  to  arrive  well  in  front  of  the  shock 
owing  to  the  intense  radiation  heating  of  the  downstream  air.  This  radiation  heated 
downstream  air  tends  to  be  "black"  and  thus  completely  obscures  the  actual  shock  front 
radiation  from  any  optical  Instrument  looking  upstream. 

The  "gas  bag  method"  proposes  to  use  a very  thin  and  transparent  walled  plastic 
bag  of  rectangular  cross-section,  filled  with  a noble  gas  to  act  as  a window  into  the 
shocked  air  region  that  is  to  be  viewed.  Theoretical  calculations  show  that  neon 
cooled  to  ~220  “K  is  the  best  candidate  material  for  this  window  to  date.  Neon  will 
be  suitably  transparent  and  still  match  the  air  parameters  sufficiently  well  to  minimize 
any  reflected  shock  heating. 
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1.3.2  Radioisotope  Source  Technique 

It  Is  well  known  that  X-ray  absorption  measurements  can  be  used  to  provide  a 
sensitive  measure  of  the  areal  mass  density  of  the  material  In  the  path  of  a collimated 
X-ray  beam.  One  approach  for  Implementing  this  technique  would  be  to  use  an  expensive 
and  relatively  fragile  X-ray  generating  machine  to  provide  a source  of  X-rays.  A 
second  approach,  which  Is  considerably  less  expensive  and  easier  to  Implement  In  very 
hostile  nuclear  generated  shock  conditions,  Is  to  utilize  commercially  available  radio- 
active X-ray  sources  Imbedded  In  a rugged  ceramic  matrix.  Owing  to  the  large  radio- 
active source  strengths  desired,  these  will  not  be  off  the  shelf  Items,  but  will  be 
prepared  to  order  In  the  size  and  shape  desired.  They  will  still  be  orders  of  magnitude 
cheaper  than  manufactured  X-ray  generating  devices.  Also  commercially  available 
radioactive  sources  generally  radiate  In  one  or  more  discrete  X-ray  or  y-vay  frequencies 
rather  than  over  a very  broad  spectral  band  as  do  the  X-ray  machines . Since  both  the 
radioactive  sources  and  detectors  are  so  compact,  there  Is  no  difficulty  In  using  them 
In  pairs,  triplets,  etc. 

1.4  Conceptual  Feasibility 

The  conceptual  feasibility  of  the  preceding  two  techniques  will  be  established  by 
the  results  presented  In  Section  3.  The  feasibility  of  designing  operating  hardware 
systems  to  Implement  the  techniques  looks  promising  but  still  requires  further  study. 

Some  specific  areas  that  require  further  study  are  outlined  In  Section  4,  Many  of 
these  problem  areas  could  be  resolved  by  Implementing  scale  tests.  Also,  the  theoretical 
study  should  be  pushed  further  by  Including  the  Important  effects  of  radiative  heating 
upon  the  shocked  and  pre-shocked  air  parameters.  First,  however,  a brief  overview  of 
the  physical  principles  Invoked  for  each  of  the  two  techniques  Is  presented  In 
Section  2 below. 


2.  THEORY:  A BRIEF  OVERVIEW 
2.1  Gas  Bag  - Pyrometer  Experiment 

In  a nuclear  driven  air  shock  the  temperatures  and  pressures  of  the  shocked  air 
can  be  expected  to  exceed  those  encountered  in  the  DABS  II-D  tests  by  at  least  an  order 
of  magnitude.  The  peak  temperatures  and  pressures  to  be  encountered  depend  of  course 
upon  the  magnitude  of  the  explosion  and  upon  the  downstream  position  of  the  observation 
station.  But  in  the  HYBLA  GOLD  test  for  example,  even  for  moderately  far  downstream 
stations  it  is  expected  that  the  temperatures  will  exceed  12,000  'K  and  the  pressures 
will  exceed  1000  atmospheres.  For  these,  and  considerably  lower  pressures  and 
temperatures,  the  air  behind  the  shock  is  expected  to  be  opaque  to  optical  radiation. 

This  opacity  is  due  to  bound-free  electron  transitions,  inverse  bremsstrahlung 
absorption,  and  the  pseudo  continuum  resulting  from  the  conjunction  of  closely  spaced 
pressure  broadened  discrete  spectral  lines  [see  Part  A,  Appendices  B and  C].  There 
also  will  be  very  appreciable  radiative  heating  of  air  well  downstresim  from  the 
propagating  shock  front.  The  air  heated  by  this  precursor  radiation  may  itself  become 
opaque,  and  radiate  as  a black  body  at  some  considerably  lower  temperature  than  that  of 
the  air  immediately  behind  the  front.  Thus  a new  technique  is  required  to  make  the 
shocked  air  temperature  measurements  penetrate  further  than  a few  centimeters  from  the 
tunnel  walls. 

The  conceptual  procedure  for  looking  deeper  into  the  shocked  air  was  to  design 
a "gas  bag"  to  act  as  an  optical  window  into  the  region  to  be  observed.  The  gas  bag 
is  to  consist  of  a thin  transparent  plastic  membrane  filled  with  an  inert  gas  that  would 
remain  transparent  to  the  spectral  wavelengths  to  be  observed.  The  noble  gases,  helium, 
neon,  argon,  etc.,  and  mixtures  thereof,  were  selected  for  inspection  since  their 
ionization  and  excitation  potentials  are  high,  and  under  ordinary  conditions  they  are 
transparent  in  broad  spectral  regions.  It  was  desired  to  match  conditions  at  the  bag- 
air  shock  interface  so  that  there  would  be  very  little  super  heating  by  shock  reflection. 
Computer  computations  were  made  with  the  TRW  ACE  code  and  with  an  acquired  JPL  code  [4]. 
It  was  found  that  the  computed  electron  density  predicted  for  the  shock  passing  through 
most  of  these  gas  bag  candidate  gases  was  much  higher  than  desirable.  These  predicted 
high  electron  densities  would  lead  to  opacity  owing  to  inverse  bremsstrahlung  absorption. 
The  surviving  candidate  gas  is  cooled  neon.  Calculations  made  with  the  JPL  code  indicate 
that  a transition  from  air  to  neon  cooled  to  ~ 220  °K  will  suffer  very  little  reflected 
I shock,  and  that  for  shocked  air  temperatures  below —9, 000  “K  and  pressures  below  ~ 500 

I atmospheres  the  neon  will  remain  transparent. 

j It  is  expected,  of  course,  that  the  plastic  membrane  containing  the  neon  will 

I vaporize  when  hit  by  the  shock.  Preliminary  approximate  diffusion  and  turbulent  mixing 
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comp.  :lons  show  however  that  the  shock  front  should  pass  the  observing  station  well 
before  any  significant  mixing  of  the  neon  and  air  occurs  to  "fog  the  window".  An 
evaluation  of  the  possible  range  of  alterations  In  state  parameters  caused  by  this 
vaporized  material  Is  one  of  the  problem  areas  mentioned  In  Section  4 for  future  Investl- 
gat Ion . 

One  of  several  possible  configurations  for  the  neon  filled  bag  "window"  Is  shown 
schematically  In  Figure  2-1.  For  this  configuration  the  radiometers  would  be  receiving 
precursor  radiation  from  the  radiation  heated  air  In  front  of  the  shock  front  until 
the  shock  front  enters  Lhelr  fields  of  view.  It  Is  assumed  that  the  radiometers  would 
be  utilizing  available  spectral  windows  where  there  Is  at  most  negligible  neon  radiation 
or  absorption.  One  such  spectral  window  useful  for  high  temperatures  exists  for  wave- 

O 

lengths  shorter  than  2595  A.  Other  useful  gaps  In  the  neon  spectrum  exist  near 

O O 

7000  A and  near  5625  A wavelengths . 

2.2  Radioisotope  Source  Experiment 

In  this  experiment  It  Is  conceived  that  two  or  more  semi -collimated  radioisotope 
sources  emitting  selected  discrete  frequency  X-rays  (or  y-rays)  will  be  slightly 
recessed  Into  one  side  wall  of  the  shock  tunnel,  while  correspondingly  filtered 
monochromatic  detectors  with  collimating  shields  will  be  placed  In  line  on  the  opposite 
wall.  This  configuration  Is  shown  schematically  In  Figure  2-2. 


CROSS-SECTION 


Figure  2-2.  Schematic  of  Proposed  Radioisotope  Source  Experiment 
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A main  consideration  of  the  design  will  concern  the  interaction  between  X-rays 
and  matter  which  takes  the  form  of:  (1)  photoelectric  absorption;  (2)  Compton 
scattering;  and  (3)  elastic  scattering.  The  second  process  results  in  a photon  of 
lower  energy,  and  thus  lower  frequency,  than  the  original  X-ray  photon  and  generally 
traveling  in  an  altered  direction.  Even  if  these  photons  are  deflected  back  into  the 
detector  by  subsequent  collisions,  they  will  not  be  passed  by  the  spectral  filter  on 
the  detector.  Multiple  elastically  scattered  X-ray  photons,  however,  are  subject  to 
being  accepted  by  the  detector  just  as  if  they  were  primary  beam  photons.  The  Compton 
scattered  photons  also  would  register  on  the  detector  unless  this  detector  system 
is  designed  for  frequency  (or  photon  energy)  discrimination.  Thus  it  is  preferable  for 
diagnostic  purposes  to  use  a monochromatic  source  system  with  a relatively  restricted 
solid  angle  output  with  a spectral  filtered  detector  system  with  a narrow  acceptance 
angle.  The  acceptance  of  any  appreciable  proportion  of  multiply  scattered  radiation, 
either  elastic  or  Compton,  would  be  treated  just  as  if  it  were  directly  transmitted 
radiation  by  the  detector  and  make  interpretation  of  the  data  difficult. 

The  attenuation  of  a beam  of  X-rays  upon  passage  through  matter  is  primarily 

dependent  upon  the  areal  density,  or  mass  per  unit  projected  areas  of  the  matter 

traversed.  For  this  reason  the  X-ray  absorption  coefficient,  B,  is  expressed  in  units 

2 -1 

of  inverse  areal  density,  usually  cm  gm  . The  areal  mass  density  parameter  generally 

symbolized  by  y (not  to  be  confused  with  wavelength  nomenclature  used  in  Part  A)  is 

-2 

given  in  units  of  gm  cm  . When  a monochromatic  beam  of  X-rays  traverses  through 
matter  of  areal  density  y,  the  Intensity  of  the  X-ray  photons  left  in  the  origin.il 
beam  obey  the  familiar  Lambert-Beer  law: 

I = I (2.1) 

o 

Here  is  the  original  beam  Intensity,  6 is  the  inverse  areal  density  absorption 
coefficient  for  X-rays  of  that  particular  energy  (or  photon  frequency),  and  y is  the 
areal  density  of  the  matter  penetrated  by  the  beam.  But,  as  previously  indicated,  an 
X-ray  detector  with  a large  acceptance  angle  will  receive  not  only  the  direct  un- 
scattered radiation  but  also  some  portion  of  multiply  scattered  radiation.  The 
corresponding  effective  attenuating  law  would  have  to  be  solved  by  iterative  tech- 
niques, but  not  in  as  simple  a relationship  as  given  by  equation  (2.1) 
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3.  FEASIBILITY  ANALYSIS 
3.1  Gas  Bag  and  Pyrometer  Feasibility 

In  order  to  determine  the  feasibility  of  the  neon  bag  "window"  a survey  of  the 
neon  spectrum  was  made  to  show  the  existence  of  wide  windows  in  the  optical  region. 
Next  the  inverse  bremsstrahlung  absorption  in  the  shock  heated  neon  was  computed. 

The  computer  code  results  are  presented  which  show  the  feasibility  of  matching  con- 
dition at  the  air-bag  Interface  so  that  there  are  no  reflected  shocks. 

From  the  spectrum  search  some  attractively  large  gaps  in  the  spectrum  of  neon 

o ■> 

were  found.  A wavelength  of  7700  A is  quite  far  from  the  nearest  line,  7839  A of 

O 

neon  which  is  quite  weak.  Another  good  wavelength  appears  to  be  5625  A;  nearest  neon 

O O 

lines  are  5589  A and  5653  A,  both  of  which  are  weak.  Monatomic  nitrogen  (which  should 

O 

be  expected  in  the  shock  heated  air  behind  the  gas  bag)  has  a weak  line  at  5617  A. 

Also  to  be  found  in  air  and  possibly  as  an  impurity  in  neon  is  argon  which  has  a weak 
line  at  7724  A« 

Line  broadening  and  Doppler  broadening  are  not  expected  to  cause  any  appreciable 

O O 

absorption  at  5625  A and  definitely  not  at  7700  A.  This  can  be  shown  as  follows: 

20  -3 

The  particle  density  of  the  shock  heated  neon  is  6 • 10  cm  . Assuming  a collison 

”16  2 5 

cross-section  of  2 • 10  cm  and  a mean  speed  of  4 • 10  cm/sec  for  neon,  one  gets  a 
collison  frequency  of  5 • 10^^  sec 
AX  « X^  — . Setting  Av  = 5 • 10^^,  X 


This  will  provide  the  line  broadening  Av,  and 

-4  5 -8 

10  cm,  one  finds  AX  = — • 10  cm  or 


about  2 A.  Thus,  one  can  expect  negligible  contributions  by  line  broadening  of  adjacent 

O O 

lines  to  absorption  at  5625  A and  7700  A. 

The  next  concern  is  the  absorption  due  to  Inverse  bremsstrahlung.  The  electron 
density  is  given  by  the  equation 


X^N 

1-X 


(gA+>  8e 
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2t\  mkT^ 


I 

kT 


v-3 


where  N is  the  number  density  (cm)  of  species  to  be  ionized  and  X its  fractional 

ionization  potential,  g^+  , g^  and  g^  are  the  statistical  weights  of  the  ionized  atom, 

the  neutral  atom  and  the  electron.  Using  the  values  appropriate  for  neon,  I • 21.47  eV, 

20  — —6 

kT  - 1 eV,  g g.+/g.  « 12,  N » 6 • 10  cm  , one  finds  X = 3 • 10  and  electron  density 
16^  3 ” 

n^ = 7 • 10  cm  . Now  in  order  to  complete  the  computation  of  inverse  bremsstrahlung 
an  estimate  of  the  electron-atom  collision  frequency  is  needed.  From  gas  discharge  data 

The  electron  speed  at  1 eV 


”16  2 

one  finds  the  electron-neon  cross  section  »2  • 10  cm 
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is  6 • 10^  cm/sec.  Then  v,  the  collision  frequency,  is  6 • 10^  cm/sec  • 2 • 10 
cm^  • 6 • 10^°  - 7 • 10^^  sec“^. 


The  absorption  coefficient  in  Inverse  bremstrahlung  is  given  by, 
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An  early  analysis  assumed  neon  in  the  gas  bag  shocked  to  1000  atmospheres  and 
12,000  °K.  Under  these  conditions,  the  gas  was  essentially  transparent  in  wide  optical 
spectral  regions  and  Inverse  bremstrahlung  was  negligible.  A more  careful  look  at  the 
Interface  conditions  between  the  bag  and  the  shock  in  air  shows  that  on  the  basis  of 
ideal  gas  theory , the  gas  in  the  bag  is  heated  to  75,000  °K.  Such  a temperature  would 
produce  excessive  ionization  and  make  the  concept  of  the  gas  bag  invalid.  A JPL  shock 
tube  computer  program,  [4]  has  been  employed  to  survey  some  gas  candidates  for  the  bag. 
The  program  treats  the  gases  as  non-ideal  and  computes  the  pressure,  temperature  and 
composition  behind  the  shock  once  the  upstre2un  conditions  and  shock  velocity  are 
specified.  The  aim  is  to  match  conditions  at  the  gas  bag-air  shock  interface  such  that 
there  is  little  shock  reflection.  This  can  be  accomplished  if  the  pressures  and 
particle  velocities  behind  the  shocks  in  the  gas  bag  and  in  air,  match.  The  program 
was  employed  to  see  how  well  this  can  be  done  with  a variety  of  gases;  helium,  neon  and 
air.  The  results  of  the  calculations  are  shown  in  Figures  3-1,  3-2  and  3-3.  Down- 
stream pressure,  temperature  and  electron  density  are  plotted  as  a function  of  shock 
velocity  for  each  upstream  gas. 

Figure  3-2  shows  that  for  neon  at  300  “K  there  is  still  a difference  in  matching 
the  Interface  conditions  (pressure  and  gas  velocity  continuous  across  interface) . A 
relationship  found  on  Page  139  of  Reference  [5],  suggests  a better  match  would  occur 
if  the  sound  velocity  in  neon  were  reduced.  This  can  be  done  by  reducing  the  tempera- 
ture. Accordingly,  the  JPL  shock  program  was  run  with  neon  at  a temperature  of  232  °K. 
The  results  as  shown  in  Figures  3-4  through  3-6  show  a much  better  interface  match 
between  air  and  neon  and  shows  that  the  temperature  is  not  too  far  off  for  a very 
good  match. 
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Figure  3-2.  Downstream  Temperature  for  Moving  Normal  Shock  Upstream 
Temperature  ■ 300  ®K,  Upstream  Pressure  " 1 Atmosphere 
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Figure  3—A  Downstream  Pressure  for  Moving  Normal  Shock  Upstream 
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Figure  3-5.  Downstream  Temperature  for  Moving  Normal  Shock  Upstream 
Temperature  ■ 300  ®K  for  Air;  232  *K  for  Neon,  Upstream  Pressure  “ 1 Atmosphere 
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The  neon  bag  system  might  be  maintained  at  this  low  temperature  In  the  field  by 

bleeding  liquid  neon  Inlzo  the  bag  and  venting  the  bag  at  the  proper  rate.  In  addition 

Figure  3-A  through  3-6  show  that  In  order  to  keep  the  electron  density  in  neon  to 

18  “3 

reasonable  levels  (<  10  cm  ) the  downstream  pressure  must  be  kept  below  500  atmospheres. 

3.2  Radioisotope  Source  System  Feasibility 

Assume  that  a high  Intensity  air  shock  is  to  be  driven  down  a tunnel  by  a nuclear  | 

explosion.  It  Is  desired  to  measure  the  density  of  the  shocked  air  and  the  mass  of  j 

lofted  material  entrained  in  the  air.  It  is  expected  that  the  air  will  be  compressed 

3 

by  a factor  of  six  or  so  in  density  and  the  amount  of  matter  lofted  will  be  ~.25  gm/cm  . 

The  time  history  of  the  resultant  areal  mass  density  is  to  be  monitored  by  X-ray  ‘ 

absorption,  with  a time  resolution  of  At,  nominally  10  ^ sec.  The  tunnel  diameter  is  = 

i 

1 meter,  but  may  change  significantly  during  the  passage  of  the  shock  which  may  cause  ’ 

a thousandfold  pressure  increase  in  the  air . j 

s 

t 

The  X-rays  are  to  be  provided  by  radioisotopes  which  are  embedded  in  the  wall  of  ^ 

the  tunnel.  The  X-rays  are  radiated  across  the  tunnel  to  a detector  imbedded  in  the  | 

tunnel  wall.  A schematic  diagram  of  this  configuration  was  presented  in  Figure  2-2.  | 

i 

When  multiple  scattering  of  X-rays  can  be  ignored  as  will  be  true  for  relatively  | 

small  areal  densities,  and  when  both  the  emission  and  detector  solid  angles  are  minimized, 
the  time  variation  of  the  received  radiation  will  be  a measure  of  the  corresponding 
time  variation  of  the  areal  mass  density  across  the  diameter  of  the  tunnel.  The  method 
of  data  reduction  of  two  or  more  such  systems  (operating  at  different  X-ray  frequencies) 
to  extract  areal  density  and  wall  displacement  information  will  be  presented  next. 

3.2.1  Reduction  of  Radioisotope  Source  Data 

Assume  a radioisotope  source  has  an  activity,  S (disintegrations/sec)  and 
produces  an  X-ray  photon  of  specified  energy  on  each  disintegration  with  a mass  attenu- 
ation coefficient  of  6 for  air,  which  we  will  assume  to  be  the  same  as  for  the  lofted 

I -2 

> material.  The  X-ray  flux  received,  in  units  of  photons  cm  , will  then  be 

i 

i F = ^ ^ (3.1) 

: 4irr 

; I 

j where  p is  the  intervening  areal  mass  density  and  r is  the  distance  between  the  source 

i and  the  detector,  initially  about  1.0  meter. 

If  r can  be  assumed  to  be  constant,  then  monitoring  the  flux,  F,  will  produce  a 
record  of  Bp  which  can  be  solved  for  p.  If  r cannot  be  taken  as  constant,  then  other 
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measurements  have  to  be  made.  One  method  explored  here  is  to  use  another  but  dif- 
ferent radioisotope  source  with  a different  S»  and  make  simultaneous  measurements  of 
the  two  fluxes  with  two  different  detectors.  Equation  (3.1)  is  rewritten  for  two 
fluxes  as  follows: 


(3.2) 


Then  simple  algebraic  manipulation  produces 


(3.3) 


If  the  u value  computed  from  equation  (3.3)  is  inserted  into  either  of  the  equations 
of  (3.2),  then  r can  be  solved  for. 

Now,  F is  not  the  quantity  measured  directly  but  N,  the  total  number  of  photons 
received  by  the  detector  of  area.  A,  in  a time.  At,  where 

N = FA  At  (3.4) 


Rewriting  equation  (3.3)  in  terms  of  and  Nj,  it  follows  that: 


To  assess  the  feasibility  of  this  concept  a estimate  of  the  random  error  inherent 
in  the  method  must  be  established.  During  a small  time  interval.  At,  the  photon  count, 

N,  will  fluctuate  by  AN.  This  happens  because  of  the  statistical  nature  of  the  dis- 
integration and  the  absorption  process.  Due  to  such  fluctuations,  the  determination 
of  p will  be  in  some  error.  An  estimate  of  this  error  now  is  made  in  general  terms  and 
subsequently,  particularized  to  the  case  of  two  specific  radioisotopes  of  chosen  strength. 
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The  assumption  is  made  that  Poisson  statistics  applies  for  the  realization  of  a gi- 
ven value  of  N.  Under  these  circumstances,  the  proper  procedure  to  determine  the  errors 
In  u Is  to  determine  the  mean  square  fluctuation  in  u from  Equation  (3.5),  assuming 
and  N2  are  uncorrelated  and  obey  the  proper  statistics.  The  procedure  is  to  differen- 
tiate Equation  (3.5)  and  replace  derivatives  by  A,  which  signifies  fluctuation. 


AU  = - 1 ("^1  "^2 


(^1-^2)  V ""l  ""2, 


(3.6) 


Now  both  sides  of  Equation  (3.6)  are  squared  and  an  ensemble  average  performed.  Then 
the  premise  of  AN^^  and  AN2  being  uncrorrelated  is  used  to  drop  cross  terms.  The  result 


is : 


<(Ap)2> 


(6 


/<^ 

1-62)^  \ 


<(ANj^)2>  ^ <(AN2)2>) 


(3.7) 


Nr 


Here,  the  symbol  < y indicates  an  ensemble  average. 

Since  Poisson  statistics  applies,  <(AN]^)2>»  and  <(AN2)^>  “ N2.  Thus: 


(3.8) 


Setting  P = ^(Au)^  , the  fractional  error  and  employing  Equations  (3.2),  (3. A)  and 


(3.3),  it  follows  that 


1 / 4Trr2  e^l2  f 

Vi  (61-32)  \ At  y 


1 + ey(e2-3i)piAi\ 

\ S2A2  A 


(3.9) 
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It  Is  seen  that,  P,  the  fractional  error  In  determining  areal  density  varies  as 
r/i/” At.  Increasing  At  by  a factor  of  10  will  decrease  the  fractional  error  by  a factor 
of  3.2,  P as  a function  of  u will  have  a minimum  value  when  0^^  m/2  * .6  (assuming 
6-  > 0-).  When  [B.m/2]  is  small,  then  P 'V'  1 . 


Now  examine  a particular  case.  Assume  as  sources  100  Curies  each  of  Cobalt  60 
(Co^®)  and  Ytterbium  169  (Yb^^^)  [half-life  of  32  days],  the  first  of  which  produces  y 
rays  of  1 MeV  energy  and  the  latter  produces  mainly  50  KeV  X-rays . Also  assume  Aj^  “ A2  “ 
10  cm^  and  At  • 10"^  sec.  At  a distance  of  1.0  meter  in  STP  air,  these  will  produce 
Nj^  = N2  = 9.4  X 103  counts  on  an  average  during  the  10“^sec,  For  50  keV  photons, 

0j^  = .20  cm2/gm  and  for  1 MeV  photons,  ° cm^/gm.  Inserting  these  values  into 
Equation  (3.8)  and  dividing  by  p one  obtains  a fraction  error  prediction  of: 


p=  0jl0?J6  L eO.ly  ,li  + » (3.10) 

^ ’^o  M 

where  r is  the  initial  distance  assumed  to  be  1.0  meter.  At  STP,  p for  1.0  meter  of 
o 

air  is  .129  gm/cm^.  Each  atmospheric  density  increment  in  a one-meter  path  length  in- 
creases this  quantity  by  .129  gm/cm^.  The  following  table  shows  the  value  of  P (the 
fractional  error  in  p determination)  as  a funtion  of  p. 


Table  3-1. 

P,  Fractional  Error  in  Areal  Mass  Determination  as  Functions  of  Areal 
(Using  100  Curie  Sources  of  Co^^  and  Yb^^^  with  10  cm^  Detectors  and 

Resolution  of  10  ^sec.) 


Mass . 
Time 
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These  errors  are  uncomfortably  large  for  small  values  of  u,  which  occur  when  the  j 

air  Is  only  slightly  compressed.  A y of  .322  corresponds  to  a density  Increase  of  a fac-  | 

tor  of  2.5  since  y Is  0.129  for  one  meter  of  air  at  STP.  These  errors,  as  pointed  out  j 

before,  can  be  reduced  by  Increasing  At  (assuming  the  time  scales  of  the  physical  proces-  j 

ses  allow  such  reduction) . The  detector  area  A can  also  be  Increased  to  decrease  P? 
however,  the  counting  rate  for  A - 10  cm^  Is  already  10®  sec“l  which  Is  getting  close  to  , 

saturating  the  fastest  currently  available  detectors.  Errors  could  also  hypothetically 
be  reduced  at  small  y by  employing  lower  energy  X-ray  photons  than  50  keV  and  thus 
augmenting  For  example,  an  11  keV  photon  has  = 4 cm^/gm.  Using  that  value  In 

Equation  (3.4)  gives  P = 0.016  for  a 100  Curie  source  at  y » .322.  This  value  of  P Is 
approximately  1/20  that  In  Table  3-1.  Corresponding  to  this  result,  a source  of  1/4  Curie 
at  11  keV  will  give  the  same  precision  as  100  Curies  at  50  keV.  Unfortunately,  radio- 
active sources  with  strengths  of  the  magnitude  of  1 Curie  at  this  energy  are  not  available. 

The  strongest  available  source  at  this  energy  Is  a Promethlum-147  activated  X-ray  source 
with  a target  of  Selenium.  The  output  with  one  Curie  of  Pm  Is  "v  1 mllllcurle  at  11  keV  In 
a width  of  1 keV. 

j Another  example  Is  Iodlne-125  with  half  life  60  days  and  with  a .027  KeV  photon.  At  ' 

■|  this  energy  6 = .4  cm^/gm.  and  P = 0.13  for  100  Curies  at  y = .32.  This  constitutes  a j 

[ reduction  In  P by  a factor  of  2.45.  The  conclusion  Is  that  16.5  Curies  of  will  be  j 

[ equivalent  to  100  Curies  of  Yb^®^.  j 

!i  i 

1'  ! 

i|  It  would  appear  at  this  stage  of  the  analysis  that  the  proper  way  to  design  the 

experiment  Is  to  relax  the  time  requirement  and  use  100  Curies  each  of  and  Co 

. : ( 

; ' t 

The  error  In  r comes  very  simply  from  Equation  (1.1)  by  taking  a logarithmic  derlva-  j 

‘I  tlve  and  using  N Instead  of  F.  Thus,  j 

k' 

M 

' ' AN 

fi  ^ _ 2 - 6,  Ay  . (3.11) 

(i  Nj  r i j 


Combining  Equations  (3.6)  and  (3.4),  one  obtains 


Squaring  and  averaging  as  before,  one  obtains: 


or 


A<(Ar)2> 

2 


r 


3i-B2 


(3.13) 


2V((Ar)2> 1_ 

r 81-82 


(3.14) 


When  n = .322  gm/cm^,  = 9400  “ N2 , 81  = .2,  61  = *06  and 

2 1/2-1  2 
[<(Ar)  >]  ' .r  = 0.016.  When  y = 20  gm/cm  , = 171.1,  N2  = 2823  and 

2 1/2  -1 

[<(Ar)  >]  .r  = 0.021.  The  projected  measurements  thus  are  quite  accurate 
over  a whole  array  of  areal  mass  densities. 


It  appears  that  the  conceptual  basis  for  the  radioisotope  source  experiment  is 
quite  sound.  A further,  engineering  systems  study  should  be  undertaken  to  determine 
if  all  the  engineering  problems  outlined  in  the  next  section  can  be  solved  and 
Implemented. 
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4.  SUGGESTED  AREAS  FOR  FURTHER  STUDY 

For  both  the  methods  outlined  in  this  report,  there  remains  the  general  problem  of 
hardware  design  to  Implement  these  techniques.  The  proposed  equipment  must  perform 
well  while  being  subject  to  extremely  hostile  environments  in  and  adjacent  to  the  shock 
tunnel  walls.  Unless  adequately  designed,  constructed  and  protected,  the  measuring 
devices  will  be  subject  to  extreme  pressure  and  thermal  loads.  The  required  survival 
time  should  allow  for  the  complete  passage  of  the  shock  front  boundary,  (~  1 msec). 

This  section  will  identify  those  areas  that  need  further  study  to  insure  the 
success  of  both  the  gas  bag  and  radioisotope  source  methods. 

4.1  Study  Areas  Associated  with  the  Gas  Bag  and  Pyrometer  Experiment 

4.1.1  Preheating  of  Downstream  Air 

Radiation  from  the  shock  front  will  preheat  the  air  in  front  of  the  shock. 

This  problem  has  not  been  accounted  for  in  the  computer  program  work  to  date.  An 
effort  should  be  undertaken  to  incorporate  these  effects  into  the  computer  code  to 
see  how  the  Interface  matching  of  temperature  is  affected. 

4.1.2  Demonstration  of  Interface  Matching 

A test  program  should  be  outlined  which  would  Involve  the  use  of  TRW’s  shock 
tube  to  demonstrate  the  feasibility  of  interface  matching  using  gas  bags  of  cooled 
neon. 

4.1.3  Determination  of  Optimum  "Gas  Bag"  Material 

It  previously  has  been  stated  that  the  bag  material  alimost  certainly  will 
vaporize  upon  contact  with  the  propagating  shock  front.  It  is  expected  that  the  volume 
of  neon  gas  should  remain  virtually  intact  until  after  the  passage  of  the  shock  front. 
Further  study  is  required  to  determine  the  time  between  evaporation  of  the  bag  material 
and  the  passage  of  the  shock  front.  The  effects  of  procursor  radiation  on  bag  survival 
should  be  examined. 

4.1.4  Chemical  Interaction  of  Vaporized  Gas  Bag  and  Air  Shock  Front 

The  atomic,  molecular  and  ionic  constituents  resulting  from  the  evaporization 
of  the  thin  bag  material  will  certainly  participate  to  some  extent  in  the  radiative 
and  chemical  processes.  Computer  studies  should  be  conducted  to  determine  what  side 
effects  this  may  produce  in  trying  to  make  pyrometer  temperature  measurements. 
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4.2  Study  Areas  Associated  with  the  Radioisotope  Source  Experiment 
4.2.1  Determination  of  Fraction  of  Lofted  Material  Entrained  in  the  Air 

The  attenuation  coefficient  B is  virtually  the  same  for  air  and  for  Si02  at 

the  X-ray  photon  energy  of  1 MeV.  At  50  KeV,  however,  the  attenuation  coefficient  for  Si02 

169 

is  approximately  50%  greater  than  that  of  air.  Yb  emits  seven  different  specific 
energy  X-ray  photons,  ranging  from  50  to  308  KeV.  It  is  thus  conceivable  that  by 
using  two  different  detectors  to  receive  the  Yb^^^  radiation,  one  with  an  aluminum 
filter,  that  the  relative  proportions  of  air  and  of  SIO2  could  be  determined.  It  is 
suggested  that  the  feasibility  and  predicted  precision  of  such  a determination  be 
investigated.  Table  4-1,  derived  from  data  given  in  Reference  [6]  is  presented  below: 
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Table  4-1.  Photon  Emissions  From  Yb  with  Associated  Absorption 
Coefficients*  of  Si,  Nj  or  O2  and  A1 


Photon  Emissions 

6 Si 

B(N2  or  O2) 

B AJ, 

-f 

n 

KeV 

cm^/gm 

cm^/gm 

cm^/gm 

1.85 

50 

.426 

.196 

.368 

0.45 

63 

.314 

.18 

.27 

0.18 

110 

.18 

.15 

.17 

0.11 

131 

.17 

.14 

.16 

0.22 

177 

.14 

.13 

.14 

0.35 

198 

.128 

.12 

.12 

0.10 

308 

.108 

.10 

.10 

*The  absorption  coefficients  at  131  and  177  were  crudely 
estimated  from  logarithmic  plots  in  Reference  1. 

tnumber  of  photons  emitted  per  disintegration 


4.2.2  Selection  of  Detector  Study 

8 —1 

The  X-ray  and  y-ray  detectors  will  have  to  count  at  a rate  of  ~10  sec  . This 
rate  will  saturate  normal  Nal  detectors.  A plastic  solution  scintillator  seems  to 
have  short  enough  decay  times,  but  a program  to  research  and  test  various  detector/ 
scintillator  systems  should  be  initiated. 
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4.2.3  Design  Study  to  Narrow  the  Emlttance  and  Acceptance  Solid  Angles 

Areal  densities  can  be  computed  simply  only  If  the  X-ray  propagates  In  a narrow 
beam  from  the  radioisotope  source  to  the  detector.  To  assure  this  baffles  will  have 
to  be  designed  and  Incorporated  Into  the  mounting  of  the  sources  and  detectors. 

Here  a design  study  program  would  be  needed. 

4.2.4  Study  of  Source  and  Detector  Survivability 

The  radioisotope  sources  and  detectors  will  be  subject  to  radiant  and  convective 

heating  by  the  shock  Ionized  air.  Air  heated  to  12,000  °K  (1  eV)  will  radiate 
4 2 

2 to  8 • 10  cal/cm  sec  to  the  Isotope  source.  A protective  shield  of  aluminum  0.1  cm 

-3 

thick  absorbing  this  much  heat  Input  will  be  vaporized  In  less  than  10  sec.  However, 

-3 

with  a reflectivity  of  90%,  It  will  reach  a temperature  of  only  293  'C  In  1.28  x 10 
sec.  Thus,  reflective  sheets  of  aluminum  (perhaps  gold  plated)  seem  adequate  pro- 
tection against  radiative  heat  transfer  from  shock  heated  air.  Estimates  of  convective 
heat  transfer  will  also  have  to  be  made  for  survivability  computations.  The  radio- 
isotope source  and  X-ray  detector  will  have  to  be  constructed  and  mounted  so  they  are 
not  crushed  or  disintegrated  by  an  overpressure  of  1000  atmospheres.  Tests  should 

be  conducted  to  varlfy  the  design  at  these  extreme  temperature  and  pressure  levels. 

i 
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